Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,362 C→A 7.5% intergenic (‑33/+38) 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,3620CA 7.5% 229.2 / 5.5 92intergenic (‑33/+38)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (43/41);  new base A (5/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAC  >  NC_000913/3471213‑3471503
                                                                                                                                                     |                                                                                                                                             
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:1178116/1‑150 (MQ=255)
                                                    tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                           >  2:1035650/1‑150 (MQ=255)
                                                    tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                           <  2:723975/150‑1 (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:879418/1‑150 (MQ=255)
                                                                          cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                                                                     >  2:290368/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:395139/1‑150 (MQ=255)
                                                                                                           cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                    >  1:340464/1‑150 (MQ=255)
                                                                                                           cTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                    >  1:1304708/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)
                                                                                                                                                     |                                                                                                                                             
GGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAC  >  NC_000913/3471213‑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: