Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,172 T→A 6.4% intergenic (+44/‑186) tufB → / → secE translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1720TA 6.4% 437.6 / 13.4 139intergenic (+44/‑186)tufB/secEtranslation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (61/69);  new base A (5/4);  total (66/73)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 141 total aligned reads displayed.

TTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGG  >  NC_000913/4177025‑4177319
                                                                                                                                                   |                                                                                                                                                   
  aCCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                                  >  1:143857/1‑149 (MQ=255)
    ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                                                                   <  7:221667/146‑1 (MQ=255)
    ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                >  4:179633/1‑149 (MQ=255)
      tGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                                                                              >  6:25558/1‑149 (MQ=255)
       gATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                                                                             >  3:386897/1‑149 (MQ=255)
       gATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGccc                                                                                                                                             >  7:369619/1‑149 (MQ=255)
           cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtttt                                                                                                                                         <  8:33290/149‑1 (MQ=255)
              ccGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttttagt                                                                                                                                      <  1:32768/149‑4 (MQ=255)
                gctCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATTCCCTTTTTGCAc                                                                                                                                    <  6:323591/147‑1 (MQ=255)
                gATCGCCATGGACGACGGTCGGCGTTTCGCAAGCCGGGACGGCGGCCGTACCGTTGGCGCGGGCGTTGGAGCAAAAGTTCTGAGCTAATTGCCGACAACATTTGACGCAAGGCGCACTAACAGGGCATCATTTGATGCCCTTTTTGCAc                                                                                                                                    <  2:343172/149‑1 (MQ=25)
                 aTCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACg                                                                                                                                   <  6:193446/149‑1 (MQ=255)
                    gcgATGGACGACGGTCTGCGGTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCtt                                                                                                                                <  5:134177/149‑1 (MQ=255)
                      gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                                                                              <  3:73857/149‑1 (MQ=255)
                      gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                                                                              <  8:205089/149‑1 (MQ=255)
                      gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGGTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                                                                              <  1:339335/149‑1 (MQ=255)
                          gacgacGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttt                                                                                                                                          >  3:193939/1‑133 (MQ=255)
                             gacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAg                                                                                                                       >  8:367346/1‑149 (MQ=255)
                              acgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAg                                                                                                                       >  7:271265/1‑148 (MQ=255)
                                ggTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAc                                                                                                                    <  3:242612/149‑1 (MQ=255)
                                  tCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCt                                                                                                                  <  5:170558/149‑1 (MQ=255)
                                  tCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCt                                                                                                                  <  8:271265/149‑1 (MQ=255)
                                   cTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                                                                                 <  3:262811/149‑1 (MQ=255)
                                    tGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                                                                                >  1:1737/1‑149 (MQ=255)
                                    tGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                                                                                >  5:83455/1‑149 (MQ=255)
                                       gTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtc                                                                                                             <  4:258475/149‑1 (MQ=255)
                                       gTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGAGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTTCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATCGAATTTTTCAACGCTTTAGTAACAGaaaatggat                                                                                                              >  2:189501/1‑141 (MQ=21)
                                           cGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatc                                                                                                          <  2:286281/148‑1 (MQ=255)
                                            gCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                        <  2:44895/149‑1 (MQ=255)
                                             cAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACAAAAAGGGCATCAAATGATGCCCttttagt                                                                                                                                      >  6:41638/1‑114 (MQ=37)
                                               aTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGGAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCGTTTTGCAAGCTTTCGTACCATACCCTGGGTCATCAGGGa                                                                                                     >  2:295883/1‑149 (MQ=38)
                                               aTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                                                                     <  7:319010/149‑1 (MQ=255)
                                                tCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAt                                                                                                    <  8:254192/149‑1 (MQ=255)
                                                  cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAttt                                                                                                  >  4:230884/1‑149 (MQ=255)
                                                   gTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtttt                                                                                                 >  7:291272/1‑149 (MQ=255)
                                                    tGACGGCGGCCGTACCGTTGGCGCGGGCGGTGTAGCAAAAGTGCTGAGCTAATAGCCGCTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                                                                                <  2:374154/149‑1 (MQ=39)
                                                    tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                                                                                <  5:76603/149‑1 (MQ=255)
                                                       aGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                                                                             <  8:267104/149‑1 (MQ=255)
                                                           ggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                                                                         >  6:194671/1‑149 (MQ=255)
                                                            gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                                                                   <  8:251329/90‑1 (MQ=255)
                                                              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATaa                                                                                      <  7:350928/149‑1 (MQ=255)
                                                              cGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATaa                                                                                      >  5:238449/1‑149 (MQ=255)
                                                                    gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTg                                                                                <  5:194671/149‑1 (MQ=255)
                                                                       ggCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGGTCATCAGTGATTTTCTTTGTCATAATCATTGCt                                                                              <  1:234185/148‑1 (MQ=255)
                                                                        gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                                                                         >  8:359438/1‑136 (MQ=255)
                                                                        gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCa                                                                                         <  7:359438/136‑1 (MQ=255)
                                                                           cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGATTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgaga                                                                          >  6:58068/1‑148 (MQ=255)
                                                                             ggCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAg                                                                       >  5:26994/1‑149 (MQ=255)
                                                                              gCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                                                                      <  5:410015/149‑1 (MQ=255)
                                                                              gCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                                                                      <  3:179633/149‑1 (MQ=255)
                                                                                gTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                                                                    >  1:64574/1‑149 (MQ=255)
                                                                                 ttGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATAATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACATGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                                                                   <  6:217214/149‑1 (MQ=255)
                                                                                       cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                             <  8:291272/149‑1 (MQ=255)
                                                                                        aaaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg                                                            >  2:251028/1‑149 (MQ=255)
                                                                                           aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc                                                         <  7:271528/149‑1 (MQ=255)
                                                                                           aGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc                                                         <  3:322889/149‑1 (MQ=255)
                                                                                             ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCg                                                        >  8:313215/1‑148 (MQ=255)
                                                                                                  aGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtata                                                    <  3:95853/147‑1 (MQ=255)
                                                                                                  aGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAAt                                                  <  4:87827/149‑1 (MQ=255)
                                                                                                  aGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttta                                                                                                                                        >  8:51465/1‑62 (MQ=37)
                                                                                                  aGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttta                                                                                                                                        <  7:51465/63‑2 (MQ=37)
                                                                                                      aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCg                                               >  3:349005/1‑148 (MQ=255)
                                                                                                      aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCg                                               >  4:14398/1‑148 (MQ=255)
                                                                                                       aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                                                                      >  4:196618/1‑124 (MQ=255)
                                                                                                       aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAgg                                                                      <  3:196618/124‑1 (MQ=255)
                                                                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGc                                                                     <  5:411740/122‑1 (MQ=255)
                                                                                                          gCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaaa                                           <  6:26994/148‑1 (MQ=255)
                                                                                                           ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGc                                         >  6:289685/1‑149 (MQ=255)
                                                                                                             gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCt                                        >  6:279207/1‑148 (MQ=255)
                                                                                                             gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTGTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCt                                        >  4:222381/1‑148 (MQ=255)
                                                                                                             gATAACATTTGACGCAATGCGCAATAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCt                                        >  1:51637/1‑148 (MQ=255)
                                                                                                              aTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                                      >  8:156128/1‑149 (MQ=255)
                                                                                                              aTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACAAGAACCTGTATCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGATAAGCTaa                                      <  8:265556/149‑1 (MQ=255)
                                                                                                               tAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                                      >  3:265228/1‑148 (MQ=255)
                                                                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGCTGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCATTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcgc                                 <  2:81110/149‑1 (MQ=255)
                                                                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcgc                                 >  5:257827/1‑149 (MQ=255)
                                                                                                                   aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcgc                                 >  4:370759/1‑149 (MQ=255)
                                                                                                                        aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTc                            >  2:220045/1‑149 (MQ=255)
                                                                                                                          gCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGa                          >  1:140247/1‑149 (MQ=255)
                                                                                                                            aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaaa                                           <  1:356149/130‑1 (MQ=255)
                                                                                                                            aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAtt                        >  2:23786/1‑149 (MQ=255)
                                                                                                                             aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAttt                       <  1:220045/149‑1 (MQ=255)
                                                                                                                                  cACTAAAAGGGCATCATTTGCTGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCCTTTCGATTTGGttt                  <  6:40528/149‑1 (MQ=255)
                                                                                                                                  cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGttt                  <  3:14398/149‑1 (MQ=255)
                                                                                                                                   aCTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtat                                                     >  5:52589/1‑113 (MQ=255)
                                                                                                                                   aCTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGtat                                                     <  6:52589/113‑1 (MQ=255)
                                                                                                                                       aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTc             <  5:279207/149‑1 (MQ=255)
                                                                                                                                         aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGa                          <  7:254673/134‑1 (MQ=255)
                                                                                                                                         aGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcgc           >  3:283632/1‑149 (MQ=255)
                                                                                                                                         aGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcgc           >  2:222036/1‑149 (MQ=255)
                                                                                                                                         aGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAGTTGGTTTGCCTcgc           >  8:159334/1‑149 (MQ=255)
                                                                                                                                            gCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAAt                                     <  7:409905/120‑1 (MQ=255)
                                                                                                                                            gCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAAt                                     >  8:409905/1‑120 (MQ=255)
                                                                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaaa                                           <  7:286629/113‑1 (MQ=255)
                                                                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaaa                                           >  8:286629/1‑113 (MQ=255)
                                                                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATc       <  7:382404/149‑1 (MQ=255)
                                                                                                                                             catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATc       >  4:25335/1‑149 (MQ=255)
                                                                                                                                              atcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCAGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATcg      >  4:51465/1‑149 (MQ=255)
                                                                                                                                                catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATcgcg    >  8:304528/1‑149 (MQ=255)
                                                                                                                                                 atTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCgg   <  3:25335/149‑1 (MQ=255)
                                                                                                                                                  tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCggg  <  7:14662/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGG  >  NC_000913/4177025‑4177319

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