Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,750 C→T 5.4% intergenic (‑133/+31) ylaB ← / ← ylaC putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7500CT 5.4% 264.6 / 4.4 93intergenic (‑133/+31)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (54/34);  new base T (3/2);  total (57/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

GTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTC  >  NC_000913/478605‑478898
                                                                                                                                                 |                                                                                                                                                    
ggcGTGTTCTCACCAGCATATCCGTATTAACACAGTGTAGGGAAAGTATATCTGATTGGGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCaga                                                                                                                                                   <  2:608415/147‑1 (MQ=255)
 tCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCagag                                                                                                                                                  >  2:676273/1‑149 (MQ=255)
  cGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTTGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                                                                                 >  1:425834/1‑149 (MQ=255)
           aCCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                                                        >  1:731666/1‑149 (MQ=255)
           aCCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAAAGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                                                                        >  1:666683/1‑149 (MQ=255)
            ccAGCTTATCCTTATTAACACAGTGTATGGAAAGTATATCTGATTGTGGATAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTGATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                                                                       <  2:592261/149‑1 (MQ=255)
            ccAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                                                                       <  2:310395/149‑1 (MQ=255)
             cAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTgg                                                                                                                                      <  1:508752/149‑1 (MQ=255)
              aGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                                                                     <  1:1033285/149‑1 (MQ=255)
               gCATATCCTTATTAACACAGTGTAGGGAAAGTATCTCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                                                                    <  2:631507/149‑1 (MQ=255)
                  tatCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGCTGGAAACGATCATTGTCATTTGACCAGTCGCAGAGCTTTACGCCTGGTATGt                                                                                                                                 >  1:624034/1‑149 (MQ=255)
                   atCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaat                                                                                                                                <  1:1026424/149‑6 (MQ=255)
                   atCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                                                                <  1:616375/149‑1 (MQ=255)
                   atCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTa                                                                                                                                <  2:465727/149‑1 (MQ=255)
                     ccTTATTAACACAGTGTAGGGAAAGTATATCTGAGTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGATACGATCTTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                                                              >  2:768307/1‑149 (MQ=255)
                       ttattaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                                                            <  1:216342/149‑1 (MQ=255)
                       ttattaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGATGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTTGAAACGATAATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                                                            >  1:659242/1‑149 (MQ=255)
                         attaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                                                          >  1:60228/1‑149 (MQ=255)
                         attaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                                                          >  2:578113/1‑149 (MQ=255)
                          ttaACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                                                                         >  2:762572/1‑149 (MQ=255)
                             acacAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCGGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACtt                                                                                                                      <  2:136687/149‑1 (MQ=255)
                              cacaGTGTCGGGAAAGTATATATCATTGAGGAAAAAGGGTGGCTGGAGATAATGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTTTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                                                                     <  1:179655/149‑1 (MQ=255)
                              cacaGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                                                                     >  2:214092/1‑149 (MQ=255)
                               acaGTGTAGGGAAAGTATATCTGATTGTGGATAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                                                                    >  1:30444/1‑149 (MQ=255)
                               acaGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                                                                    <  1:676273/149‑1 (MQ=255)
                               acaGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGAAGGGTGGAAACGATCATTGTCATTTGACCACGCGCAGAGCTATACGCCTGGTATGTAACCCTAACTTaa                                                                                                                    >  2:360393/1‑149 (MQ=255)
                                   tgtAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTACCCCTAACTTAAGATg                                                                                                                >  2:697146/1‑149 (MQ=255)
                                    gtAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGt                                                                                                               <  1:869414/149‑1 (MQ=255)
                                     tAGGGAAAGTATATCTGATTGTGGAAAAAGGGTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                                                              <  2:870328/149‑1 (MQ=255)
                                     tAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                                                              >  1:213552/1‑149 (MQ=255)
                                      aGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                                                             >  2:717660/1‑149 (MQ=255)
                                      aGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                                                             >  1:1077738/1‑149 (MQ=255)
                                         gAAAGTATATCTGATTGTGGAAAAAGCCTGGCTGGAGATCAGGAGGTAGAGCTGGTATTTAATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGCGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAtt                                                                                                          <  1:1064026/149‑1 (MQ=255)
                                          aaaGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGCGGTAGAGCTGGTATTTCATGAAGTGATGTGTGGAAACGATCATTGTCATTTGAACAGGCGCAGAGCTTTACGCATGTTATGTAACCCTAACTTAAGATTATGATTc                                                                                                         >  1:739554/1‑149 (MQ=255)
                                          aaaGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACACTAACTTAAGATGTTGATTc                                                                                                         >  2:983591/1‑149 (MQ=255)
                                            aGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCgg                                                                                                       <  1:403965/149‑1 (MQ=255)
                                                     gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGATGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCTGCgcgagcga                                                                                              >  1:319930/1‑149 (MQ=255)
                                                     gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGATACGATCATAGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgccagcgc                                                                                              >  2:861216/1‑148 (MQ=255)
                                                     gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                                                              >  1:328557/1‑149 (MQ=255)
                                                     gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                                                              >  2:197818/1‑149 (MQ=255)
                                                       ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                                                            >  1:78006/1‑149 (MQ=255)
                                                       ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                                                            <  1:214092/149‑1 (MQ=255)
                                                       ttGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTGGATTCGGCGCGAGCGAgg                                                                                            <  1:1138/149‑1 (MQ=255)
                                                            tataaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGATGTGATGGGTGGAAACGATCTTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                                                                       <  2:624034/146‑1 (MQ=255)
                                                             aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                      >  2:526670/1‑149 (MQ=255)
                                                             aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                      >  1:44355/1‑149 (MQ=255)
                                                             aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                      >  1:466987/1‑149 (MQ=255)
                                                             aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                                                                      >  2:1152281/1‑149 (MQ=255)
                                                              aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCGTAACTTAAGAGGTTGATTCGGCGCGAGCGAGGGTaaaaa                                                                                     >  2:183163/1‑149 (MQ=255)
                                                              aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                                                                     >  1:253767/1‑149 (MQ=255)
                                                                  ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                                                                 <  2:544556/149‑1 (MQ=255)
                                                                   gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGTCGCGAGCGAGGTTAAAAATATCg                                                                                >  2:327689/1‑149 (MQ=255)
                                                                           gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCGTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                                                                        >  1:613725/1‑149 (MQ=255)
                                                                            agaTAAGGACGTAGAGCTGGTATTTCATCAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCAGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                                                                       <  1:983582/149‑1 (MQ=255)
                                                                              aTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGATGATTCGGCTCGAGCGAGGGTAAAAATATCGTAAAATGACAg                                                                     >  1:1149260/1‑149 (MQ=255)
                                                                                  ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                                                                 >  2:982577/1‑149 (MQ=255)
                                                                                  ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                                                                 >  1:214651/1‑149 (MQ=255)
                                                                                  ggaggTAGAGCGGGTATTTCATGAAGTGATGTGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGTTGTTGATTCGGCGCGATCGAGGGTATAAATATCGTATAAATACAGTTCa                                                                 >  2:231454/1‑149 (MQ=255)
                                                                                       tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                                                            <  1:536308/149‑1 (MQ=255)
                                                                                         gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                                                          >  1:789539/1‑149 (MQ=255)
                                                                                            cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa                                                       >  1:148738/1‑149 (MQ=255)
                                                                                            cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa                                                       >  2:254031/1‑149 (MQ=255)
                                                                                             tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGTGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc                                                      <  2:1094427/149‑1 (MQ=255)
                                                                                                        gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt                                           >  2:355655/1‑149 (MQ=255)
                                                                                                              aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAgtt                                     <  1:1044127/149‑1 (MQ=255)
                                                                                                                  gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAACAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTc                                 <  1:143023/149‑1 (MQ=255)
                                                                                                                  gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTc                                 <  2:1173087/149‑1 (MQ=255)
                                                                                                                  gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTc                                 <  2:108273/149‑1 (MQ=255)
                                                                                                                  gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTc                                 <  2:1077738/149‑1 (MQ=255)
                                                                                                                   tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                                <  2:425834/149‑1 (MQ=255)
                                                                                                                   tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                                >  1:102642/1‑149 (MQ=255)
                                                                                                                   tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTcc                                >  2:152002/1‑149 (MQ=255)
                                                                                                                    ggAAAAGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAATACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCt                               <  2:1145054/149‑1 (MQ=255)
                                                                                                                     gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCtt                              >  2:337413/1‑149 (MQ=255)
                                                                                                                     gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCtt                              >  1:757217/1‑149 (MQ=255)
                                                                                                                      aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTa                             >  1:57846/1‑149 (MQ=255)
                                                                                                                       aaCGATCATTGTGATTTGACCAGGCGCAGAGCTTTTCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTAt                            <  2:1177704/149‑1 (MQ=255)
                                                                                                                        aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACTCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATg                           <  1:578113/149‑1 (MQ=255)
                                                                                                                         cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGa                          <  1:1152281/149‑1 (MQ=255)
                                                                                                                          gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGaa                         >  2:660457/1‑149 (MQ=255)
                                                                                                                                    aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac               >  1:84660/1‑149 (MQ=255)
                                                                                                                                    aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac               >  2:702851/1‑149 (MQ=255)
                                                                                                                                    aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGac               >  2:206235/1‑149 (MQ=255)
                                                                                                                                       ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcg            >  1:583006/3‑149 (MQ=255)
                                                                                                                                         cccAGGCGCAGCGCTTTCCGCCTGGTATGTAACCCTAACATCAGCTGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTACAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg          <  1:768307/148‑1 (MQ=255)
                                                                                                                                         aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg          >  1:467463/1‑149 (MQ=255)
                                                                                                                                         aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg          >  1:450432/1‑149 (MQ=255)
                                                                                                                                         aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACAcgcg          >  2:429906/1‑149 (MQ=255)
                                                                                                                                          ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGCCGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCgg         >  2:1019716/1‑149 (MQ=255)
                                                                                                                                          ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCgg         >  1:127504/1‑149 (MQ=255)
                                                                                                                                           cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGga        <  1:254031/149‑1 (MQ=255)
                                                                                                                                           cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGga        >  1:976475/1‑149 (MQ=255)
                                                                                                                                           cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGga        >  2:106499/1‑149 (MQ=255)
                                                                                                                                              gcgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGagag     <  2:466987/149‑1 (MQ=255)
                                                                                                                                                 cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTc  <  2:731666/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTTGCAGTTGTTCCTTATGAACATCCGCGACACGCGGAGAGTTC  >  NC_000913/478605‑478898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: