Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,357 C→T 7.6% intergenic (‑28/+43) tufA ← / ← fusA translation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,471,3570CT 7.6% 267.9 / 8.4 93intergenic (‑28/+43)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base T (4/3);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

CCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAC  >  NC_000913/3471208‑3471503
                                                                                                                                                     |                                                                                                                                                  
ccgTAGTTTTTATCCAGTTCGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTATACGTGGCCGATAGTACGAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTTTATTCGTTACTATAGTGCTCTcc                                                                                                                                                    <  2:940272/150‑1 (MQ=18)
 cgTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTccc                                                                                                                                                   >  2:832492/1‑150 (MQ=255)
  gTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcccc                                                                                                                                                  <  1:541076/150‑1 (MQ=255)
     ggTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCATATTTTTCTTTAGACACGGCTTTATTCCTTACTATAGTGCTCTCCCCGTc                                                                                                                                               >  2:1209589/1‑150 (MQ=255)
     ggTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCACCTTc                                                                                                                                               >  1:705390/1‑150 (MQ=255)
      gTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCa                                                                                                                                              >  2:54666/1‑150 (MQ=255)
       ttttAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACTTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGATCTCCCCTTCAg                                                                                                                                             <  1:148725/150‑1 (MQ=255)
       ttttAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                                                                             <  1:378289/150‑1 (MQ=255)
        tttAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAgg                                                                                                                                            <  1:66786/150‑1 (MQ=255)
          tAGCCAGTACGGTGGTGATTGCGGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                                                                          >  1:1206494/1‑150 (MQ=255)
             ccAGTACGGTGGTGATTGCAGCGGTCAGTGTAGTTTTACCGTGGGCAACGTGGCCGATAGTACCAACGTTTACGTGAGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                                                                       <  2:558884/150‑1 (MQ=255)
             ccAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                                                                       <  1:1047649/150‑1 (MQ=255)
             ccAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGATCAAATTTTTCTTTAGACACGCATAAATTCATTAATATAGTCATCTCCCCTTCATgagaga                                                                                                                                       >  1:167709/1‑150 (MQ=18)
              cAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGagagag                                                                                                                                      <  2:833709/150‑1 (MQ=255)
               attACGGTGGTTATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTTGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGc                                                                                                                                     <  2:561534/148‑1 (MQ=21)
                gTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                                                                    >  1:186904/1‑150 (MQ=255)
                ctACGGGGGTGAATTCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATCTTCCTTACAATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                                                                    <  2:746365/149‑1 (MQ=25)
                 tACGGTGGTGATTGCAGCGGTGAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                                                                   <  1:12665/150‑1 (MQ=21)
                   cGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCGGGAGAGAGCACgg                                                                                                                                 <  2:23034/150‑1 (MQ=255)
                   cGGTGGTGATTGCAGCGGTCAGAGTAGGTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                                                                 <  1:667894/150‑1 (MQ=255)
                   cGGGGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACgg                                                                                                                                 <  1:1113679/150‑1 (MQ=255)
                    ggtggtGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                                                                >  2:470963/1‑150 (MQ=255)
                      gggtgATTGCAGCGGCCAGAGTAGTTTTACCGTGGTCAACGTGGAATATATTACCTACGTTAACGTGATGTTTTGTACGTTTAAATTTTTCTTTAGACTCGGCTATATTCATTACTATAGTGCTCTATCCTGAAGGGGAGAGCActatag                                                                                                                              <  2:61384/149‑6 (MQ=255)
                       ggtgATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGCCt                                                                                                                             >  2:1033425/1‑150 (MQ=255)
                       ggtgATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACt                                                                                                                             >  1:458997/1‑150 (MQ=255)
                        gtgATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAAGGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTAAAATTTTTCTTTAGACACGGCTATATTCCTACTTTAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                                                                           <  2:558702/150‑1 (MQ=255)
                        gtgATTGCAGCGGTACGATGAGTTTAACCGTTGGCAACGTGGCCGATAGTACCAACGTTACCGTGCGGTTTTGTACGGGCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACtt                                                                                                                            <  1:246284/150‑1 (MQ=255)
                             tGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                                                                       >  1:370623/1‑150 (MQ=255)
                               cAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAtt                                                                                                                     >  2:778185/1‑150 (MQ=255)
                                aGCGGTAAGAGTAGTTTTACCGTGGTCAACGTGGCAGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGAGACGGCTTTATTACTGACTATAGTGCTCTCCCATTCAGGAGAGAGAACGGGACTTTGGAATTa                                                                                                                    >  2:194394/1‑150 (MQ=25)
                                   ggTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                 <  2:534837/150‑1 (MQ=255)
                                   ggTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                 <  1:1246097/150‑1 (MQ=255)
                                       agagTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTa                                                                                                             >  1:650861/1‑150 (MQ=255)
                                        gagTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAg                                                                                                            <  2:751130/150‑1 (MQ=255)
                                           tagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCt                                                                                                         <  1:557677/150‑1 (MQ=255)
                                            agtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCtt                                                                                                        >  1:865468/1‑150 (MQ=255)
                                            agtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGAACGGGACTTTGGTATTAACCCTTAGGCtt                                                                                                        >  1:908034/1‑150 (MQ=255)
                                                ttACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAttt                                                                                                    <  1:1211489/150‑1 (MQ=255)
                                                       ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                                                                             >  2:468663/1‑150 (MQ=255)
                                                        gtaaaCGTGGCCTATAGTACCAACTTTAAGGTGGGGTTTTTTATGTTCAAATTTTTATTTAGACACGGTTATATTCCTTATTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                                                                            <  2:1174661/147‑1 (MQ=255)
                                                        gTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                                                                            >  1:1164752/1‑150 (MQ=255)
                                                         tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                           <  2:723975/150‑1 (MQ=255)
                                                         tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                           >  2:1178116/1‑150 (MQ=255)
                                                         tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                           >  2:1035650/1‑150 (MQ=255)
                                                          cAACGTGTCCGATAGTACCAACTTTAACGGGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCt                                                                                          <  1:102687/150‑1 (MQ=255)
                                                            aCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATAATCCTTACTATAGTGCTCTCCCCTTCAGTTGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTc                                                                                        <  2:1127744/150‑1 (MQ=255)
                                                             cGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCa                                                                                       <  2:357418/150‑1 (MQ=255)
                                                               tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGATCTCCCCTTCAATAGAGAGCACGGGACTTTGGTATTAAACATTAGGCTTATTTACCACGTGCTTCaaa                                                                                     >  1:249621/1‑149 (MQ=37)
                                                                 gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTACCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTa                                                                                   >  2:989227/1‑150 (MQ=255)
                                                                 gCCGAAAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGAGCTCTCCACTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTa                                                                                   >  1:97715/1‑150 (MQ=255)
                                                                  ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTAc                                                                                  >  1:1096772/1‑150 (MQ=255)
                                                                  ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTAc                                                                                  >  1:562946/1‑150 (MQ=255)
                                                                     aTAGTACCAACGTTAACGGGCGGTTTTGTCCGTTCAAATTTTTCTTTAGTCACGGCTATATTCCTTACTATAGTGCGCTCCCCATCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                                                               <  1:505691/150‑1 (MQ=255)
                                                                      tattACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGAGACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGcc                                                                              <  1:931768/147‑1 (MQ=255)
                                                                       aGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCt                                                                             >  2:289081/1‑150 (MQ=255)
                                                                          aCCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAg                                                                          >  2:1210561/1‑150 (MQ=255)
                                                                            cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCa                                                                        >  2:272449/1‑150 (MQ=255)
                                                                            cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCa                                                                        >  1:493685/1‑150 (MQ=255)
                                                                             aaCGTAAACGTGGGGTTTTGTACGTTCAAATTTTTCTTTATACACGGTTAAATTCCTTACTATTGTGCGCTCCCCTTCCGTAGTGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCaa                                                                       <  2:966771/150‑1 (MQ=255)
                                                                              aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAAc                                                                      >  1:1201361/1‑150 (MQ=255)
                                                                               cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                                                                     >  2:290368/1‑150 (MQ=255)
                                                                               cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                                                                     >  2:879418/1‑150 (MQ=255)
                                                                               cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGAACGGGACTTTTGTATTAACCCTCAGGCTTATTTACAACGGGCTTCAATTACGGCCTGAGCAACg                                                                     >  1:878477/1‑150 (MQ=255)
                                                                                    aCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTAATATAGTGATCTCCCCTTAAGGAGAGATCTCGGGAATTTGCTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACt                                                                >  2:719055/1‑150 (MQ=255)
                                                                                     cGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCATGAGCAACGTTACTc                                                               >  1:696323/1‑150 (MQ=255)
                                                                                           gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt                                                         >  1:565156/1‑150 (MQ=255)
                                                                                            ggttGTACGTTAAAATTTTATTTTATACACGTATATATTCATTACTATACTGCTATGCACTTCAGTAGAGAGCACGTGACTTTGGTATTACCCCTTAGGCTTATTTACCAGGGGCTTAATTTAGTGCCTGAGCAACGTTTATCTGCGCtt                                                        <  2:159044/148‑1 (MQ=255)
                                                                                              ttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca                                                      <  1:854598/150‑1 (MQ=255)
                                                                                              ttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca                                                      <  1:1141847/150‑1 (MQ=255)
                                                                                               tGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCACCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat                                                     >  1:1272528/1‑150 (MQ=255)
                                                                                                   cGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTTTTTACAACGGGCTTCAATTACGGCCTGAGAAACGTTACTAGGCGCTTCATCATa                                                 >  1:1287717/1‑150 (MQ=255)
                                                                                                         aaTTTTTCTTTATACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAg                                           <  1:282749/150‑1 (MQ=255)
                                                                                                            ttttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGaa                                        <  2:370623/150‑1 (MQ=255)
                                                                                                               tCTTTAGACACGGCTATATTCCTTACTATAGTGCTCACCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTc                                     <  1:301672/150‑1 (MQ=255)
                                                                                                                cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                    <  2:186904/150‑1 (MQ=255)
                                                                                                                cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                    >  1:1304708/1‑150 (MQ=255)
                                                                                                                cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                    >  1:340464/1‑150 (MQ=255)
                                                                                                                cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                    >  1:395139/1‑150 (MQ=255)
                                                                                                                   tAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATa                                 >  1:553740/1‑150 (MQ=255)
                                                                                                                    aGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAg                                >  2:1081090/1‑150 (MQ=255)
                                                                                                                      acacGGATATATTCCTTACTATCGGGCTCTCCCCATCTGGAGAGAGAACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAgtg                              <  1:1144934/150‑1 (MQ=255)
                                                                                                                       cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGCTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTGACTCGGCGCTTCATCAGACGTCAGGAATTCCATTgtgt                             >  1:70413/1‑150 (MQ=255)
                                                                                                                         cGGCTATATTCCTTCCTATAGTGGTCTCCCCTTCTGGAGCGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTat                           <  2:1128388/150‑1 (MQ=255)
                                                                                                                         cGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACACTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTat                           >  1:162318/1‑150 (MQ=255)
                                                                                                                            cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTTCCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTGAGGAATTCCATAGTGTatgat                        >  1:471606/1‑150 (MQ=255)
                                                                                                                             tataTTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGTGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTatgatg                       <  2:167709/150‑1 (MQ=255)
                                                                                                                               tatTCCTTACTATAGTGCTCTCCCGTTCATGAGAGAGCACGGGACTTTGGTATTAACCCTTATGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATTATGCa                     <  1:736782/150‑1 (MQ=255)
                                                                                                                                         tcccgTGCTCTCTCCTGATGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATGATACGTCAGGAATTCCATAGTGTATGATGCGCGACCTTggg           >  1:69216/5‑150 (MQ=255)
                                                                                                                                          ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGt          <  1:190310/150‑1 (MQ=255)
                                                                                                                                           ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTc         >  2:485114/3‑150 (MQ=255)
                                                                                                                                            cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGGGTATGATGCACGACCTTTGGTCa        >  2:310787/2‑150 (MQ=255)
                                                                                                                                              tGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCaga      >  1:12039/1‑150 (MQ=255)
                                                                                                                                              ggcgccccccTTCCGTAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCaga      <  1:961356/144‑1 (MQ=255)
                                                                                                                                              ggCTCTCCCCGTCCGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCaga      <  2:190214/149‑1 (MQ=255)
                                                                                                                                               gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCagag     >  2:2515/1‑150 (MQ=255)
                                                                                                                                                ctctCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAgaga    <  2:370945/150‑1 (MQ=255)
                                                                                                                                                  ctcCCCTTCACGATAGAGCACGGGACTTTTGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAc  <  2:1061602/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAC  >  NC_000913/3471208‑3471503

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: