Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,024,382 T→A 5.2% N147I (AAT→ATT)  rfaH ← transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,024,3820TA 5.2% 378.4 / 6.3 125N147I (AAT→ATT) rfaHtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (59/59);  new base A (3/4);  total (62/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

CCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTTATATAC  >  NC_000913/4024233‑4024530
                                                                                                                                                     |                                                                                                                                                    
cccATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGAGTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTttat                                                                                                                                                      >  2:1342284/1‑150 (MQ=255)
   atatTTTGCAACGTATTGCGCAATGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGTGATTGCTTTTTAGAGTTTGCTGAACTCGGTATTCTTCACACTGTGCTTAATCTCTttattta                                                                                                                                                   <  1:363071/150‑3 (MQ=255)
        ttGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGa                                                                                                                                              <  2:1536252/150‑1 (MQ=255)
         tGCAACGTATTGCGCACTGGCGTGGAGAAGATGACGCATGGCGGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGCGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGAt                                                                                                                                             <  2:300161/150‑1 (MQ=255)
           cAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGAttt                                                                                                                                           >  1:272508/1‑149 (MQ=255)
           cAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTATGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTa                                                                                                                                           <  2:1929812/150‑1 (MQ=255)
              cGTATTCCGCCCTGGCGCGGAGACGATGCCGCATGGCTGGCTGCCCCCACGGATGCTAATCTCACAACACAGATTCGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCATCACACTGTGCTTAATCTTATTAATAAAGAGATTAAGc                                                                                                                                        <  2:1386657/150‑1 (MQ=255)
               gTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCa                                                                                                                                       >  1:1144892/1‑150 (MQ=255)
                 aTTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAAt                                                                                                                                     <  1:1188080/150‑1 (MQ=255)
                   tGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAATCACTGGGTGGGATTGCGTTTTAGAGGTTGCGGAACTCGGTATGCTTCACACTGTGCTGAATCTCTTTATTAATAAGATTAAGCAATAg                                                                                                                                   >  2:1387728/1‑150 (MQ=255)
                   tGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAg                                                                                                                                   >  2:1077625/1‑150 (MQ=255)
                    gcgcACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTTGTTAATCTCTTTATTAATAAGATTAAGCAATAga                                                                                                                                  >  2:313818/1‑149 (MQ=255)
                     cgcACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTTATTAATAAAGAGATTAAGCacagtgt                                                                                                                                 <  2:1074430/150‑7 (MQ=39)
                     cgcACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTTATTAATAAAGAGATTAAGCacagtgt                                                                                                                                 <  2:662356/150‑7 (MQ=39)
                           ggCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAAGACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                                                                                           >  1:102964/1‑150 (MQ=255)
                           ggCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                                                                                           >  2:250767/1‑150 (MQ=255)
                           ggCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                                                                                           <  1:1754482/150‑1 (MQ=255)
                                    aagaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGTGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCAc                                                                                                                  <  2:1799668/150‑1 (MQ=255)
                                    aagaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCAc                                                                                                                  <  1:321545/150‑1 (MQ=255)
                                     agaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCAcc                                                                                                                 >  2:1400268/1‑150 (MQ=255)
                                        tGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATc                                                                                                              >  1:1860495/1‑150 (MQ=255)
                                          ccGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCgg                                                                                                            >  1:915254/1‑150 (MQ=255)
                                          ccGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCgg                                                                                                            <  1:1194503/150‑1 (MQ=255)
                                           cGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCggg                                                                                                           <  1:1550036/150‑1 (MQ=255)
                                           cGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCggg                                                                                                           <  2:444692/150‑1 (MQ=255)
                                           cGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACAGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCggg                                                                                                           >  1:1340756/1‑150 (MQ=255)
                                            gCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGt                                                                                                          <  1:906180/150‑1 (MQ=255)
                                               tggctggcTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCg                                                                                                       >  1:551016/1‑150 (MQ=255)
                                                  ctggctaCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCTATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTg                                                                                                    <  1:306035/150‑1 (MQ=255)
                                                     gctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTTTTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGaaa                                                                                                 <  2:1478952/150‑1 (MQ=255)
                                                      ctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGaaaa                                                                                                <  2:1221403/150‑1 (MQ=255)
                                                      ctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCATTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGaaaa                                                                                                <  2:567569/150‑1 (MQ=255)
                                                        gCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATg                                                                                              <  2:1210577/150‑1 (MQ=255)
                                                          caccacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATCGc                                                                                            >  2:268912/1‑150 (MQ=255)
                                                           tccacGGATGCTAATGTCAAAACACTGGTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGcc                                                                                           >  1:1382074/2‑150 (MQ=255)
                                                           accacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACAATCGGGTTCGGTGAAAATGGca                                                                                           >  1:1951734/1‑149 (MQ=255)
                                                                   tGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGcc                                                                                   >  2:1252977/1‑150 (MQ=255)
                                                                    gCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCt                                                                                  >  2:2042531/1‑150 (MQ=255)
                                                                     cTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCtt                                                                                 >  2:1477397/1‑150 (MQ=255)
                                                                      tAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTc                                                                                >  1:376728/1‑150 (MQ=255)
                                                                       aaTGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCg                                                                               <  1:2025856/150‑1 (MQ=255)
                                                                        aTGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGa                                                                              >  2:1550465/1‑150 (MQ=255)
                                                                        aTGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGa                                                                              >  1:1460404/1‑150 (MQ=255)
                                                                         tGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGaa                                                                             >  1:1642737/1‑150 (MQ=255)
                                                                            cAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAAcgc                                                                          >  2:1384802/1‑150 (MQ=255)
                                                                              aaaCACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACgcgc                                                                        >  2:1330393/1‑150 (MQ=255)
                                                                              aaaCACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACgcgc                                                                        >  2:1387048/1‑150 (MQ=255)
                                                                                     gTTTGGGCTTGCGTTGTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCGCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt                                                                 <  1:1634034/150‑1 (MQ=255)
                                                                                     gTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt                                                                 <  1:1834792/150‑1 (MQ=255)
                                                                                     gTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt                                                                 <  1:493168/150‑1 (MQ=255)
                                                                                                ggTTTTCGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCt                                                      <  1:492732/149‑1 (MQ=255)
                                                                                                ggTTTTAGATTTTGCTGAACTCTGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGATCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCt                                                      <  2:1981025/149‑1 (MQ=255)
                                                                                                  ttttAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTa                                                    <  2:915253/150‑1 (MQ=255)
                                                                                                    ttAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATc                                                  >  2:86170/1‑150 (MQ=255)
                                                                                                       gagTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCAcc                                               >  1:1598385/1‑150 (MQ=255)
                                                                                                        agTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCAccc                                              >  2:2006963/1‑150 (MQ=255)
                                                                                                         gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCTCCCg                                             <  2:1157617/150‑1 (MQ=255)
                                                                                                         gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCg                                             >  1:1798537/1‑150 (MQ=255)
                                                                                                         gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCg                                             >  1:1440878/1‑150 (MQ=255)
                                                                                                         gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCATCGGTAATAATCACCTTATCACCCg                                             >  1:1349067/1‑150 (MQ=255)
                                                                                                         gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCAAGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCAGAAACCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCg                                             >  2:170039/1‑150 (MQ=255)
                                                                                                          tttGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCgg                                            <  2:452587/150‑1 (MQ=255)
                                                                                                           ttGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCCTCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGa                                           >  2:1551927/1‑150 (MQ=255)
                                                                                                            tGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGAt                                          >  2:514128/1‑150 (MQ=255)
                                                                                                            tGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGAt                                          >  1:1915584/1‑150 (MQ=255)
                                                                                                                  aCTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAgggg                                    >  1:428196/1‑150 (MQ=255)
                                                                                                                     cGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTg                                 <  2:767840/150‑1 (MQ=255)
                                                                                                                      ggTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGc                                >  1:894706/1‑150 (MQ=255)
                                                                                                                        caTTCCTCACACTGTGCTTAATCTATTAAATATACAGATTAAGAAATAGCATGGAGCGAGCCTCACAATCGGGTTCGTTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCg                              <  2:1928331/149‑1 (MQ=255)
                                                                                                                            cttcACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATc                          >  2:882015/1‑150 (MQ=255)
                                                                                                                             ttcACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCg                         >  1:1735604/1‑150 (MQ=255)
                                                                                                                               cacaCTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGAc                       >  2:1922527/1‑150 (MQ=255)
                                                                                                                                acacTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACa                      <  1:1978841/150‑1 (MQ=255)
                                                                                                                                   cTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATg                   <  2:1835002/150‑1 (MQ=255)
                                                                                                                                   cTGTGCTTAATCTCTTTATTAATAAGAATAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATg                   <  1:587359/150‑1 (MQ=255)
                                                                                                                                    tgtgCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGt                  >  1:653376/1‑150 (MQ=255)
                                                                                                                                    attgCTTAATCTTATTAATAAAGAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGt                  >  2:709932/3‑150 (MQ=255)
                                                                                                                                       gCTTAATCTTATTAATAAAGAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTTCACCCGGATAAGGGGTTGCCGGATCGACAATGTCttt              >  2:629534/1‑150 (MQ=255)
                                                                                                                                       gCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCtt               >  2:1668639/1‑150 (MQ=255)
                                                                                                                                         ttAATCTTATTAATAAAGAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTc             >  2:245801/1‑150 (MQ=255)
                                                                                                                                         ttAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTc             <  1:1893046/150‑1 (MQ=255)
                                                                                                                                         ttAATCTCTTTATAAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCAACAATGTCTTTc             <  2:272508/150‑1 (MQ=255)
                                                                                                                                          tAATCTCTTTATTAATAAGATTAAGCAATAGCCTGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGTCCTGAAAGCCTTCGAACGCGCCTTCTGTAATAATCACCTTTTCACCCGGATAATGTGTTGCCGGATCGACAATGTCTTTCg            <  1:1441331/150‑1 (MQ=255)
                                                                                                                                          tAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCg            >  1:1992796/1‑150 (MQ=255)
                                                                                                                                          tAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCg            >  1:132125/1‑150 (MQ=255)
                                                                                                                                           aaTCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCgg           <  1:2027398/150‑1 (MQ=255)
                                                                                                                                           aaTCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCgg           <  2:779910/150‑1 (MQ=255)
                                                                                                                                             tctcTTTATTAATATGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTTAAAATGGCCTGAAAGCCTTCGAACGGGCCTTCGTTAATAATCACCTTATCACCCGTATAAGGGGTTGCCGGATCGACAATGTCTTTCGGtt         <  1:837953/150‑1 (MQ=255)
                                                                                                                                             tctcTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGtt         <  1:1036924/150‑1 (MQ=255)
                                                                                                                                              ctctTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGttt        >  2:386812/1‑150 (MQ=255)
                                                                                                                                              ctctTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGttt        <  1:1184881/150‑1 (MQ=255)
                                                                                                                                               tctTTATTAATAAGATTAATCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGACTTCGAACGCGCGTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTta       <  2:302375/150‑1 (MQ=255)
                                                                                                                                               tctTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTta       <  2:102964/150‑1 (MQ=255)
                                                                                                                                               tctTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCATATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGATta       >  1:1122846/1‑150 (MQ=255)
                                                                                                                                                ctTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGCATCGACAATGTCTTTCGGTTtat      >  1:64804/1‑150 (MQ=255)
                                                                                                                                                ctTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCGGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTtat      <  1:981930/150‑1 (MQ=255)
                                                                                                                                                 tttATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTtata     <  1:823239/150‑1 (MQ=255)
                                                                                                                                                 tttATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTtata     >  1:825024/1‑150 (MQ=255)
                                                                                                                                                    attaATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTTATATAc  >  2:747330/1‑150 (MQ=255)
                                                                                                                                                    attaATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGAGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAAAGTCTTTCGGTTTATATAc  >  2:1171116/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTTATATAC  >  NC_000913/4024233‑4024530

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: