Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,024,378 T→A 5.2% K148N (AAA→AAT rfaH ← transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,024,3780TA 5.2% 352.8 / 6.3 115K148N (AAA→AATrfaHtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (53/56);  new base A (3/3);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 129 total aligned reads displayed.

TGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTTATA  >  NC_000913/4024230‑4024527
                                                                                                                                                    |                                                                                                                                                     
acccccATATTTCGCAACGTATTGCGCACTGGCGTGGAGAAGAAGCCGCATGGCGGGCTGCCACCACGGATCCTAATGGCCACACACTGTTTGGGATTGCGTGTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCtt                                                                                                                                                      <  2:134407/148‑1 (MQ=255)
 gCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAGGGCTGGCTGCCACCACGGATGCGAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCttt                                                                                                                                                     <  2:968732/150‑1 (MQ=255)
   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                                                                                                                        <  1:1754482/150‑1 (MQ=255)
                              ggCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                                                                                        >  2:250767/1‑150 (MQ=255)
                                 gTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagcgag                                                                                                                     <  2:855448/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                                                                                                        <  2:444692/150‑1 (MQ=255)
                                              cGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCggg                                                                                                        <  1:1550036/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)
                                                   ggctggctGCCACCACGGATGCTAATGTCAAAACACTGTTTTGGATTGCGTTTTAGAGTATGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCgg                                                                                                   <  2:1249406/150‑1 (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)
                                                                       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                                                                           >  1:1460404/1‑150 (MQ=255)
                                                                           aTGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGa                                                                           >  2:1550465/1‑150 (MQ=255)
                                                                            tGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGaa                                                                          >  1:1642737/1‑150 (MQ=255)
                                                                               cAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAAcgc                                                                       >  2:1384802/1‑150 (MQ=255)
                                                                                 aaaCACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACgcgc                                                                     >  2:1387048/1‑150 (MQ=255)
                                                                                 aaaCACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACgcgc                                                                     >  2:1330393/1‑150 (MQ=255)
                                                                                        gTTTGGGCTTGCGTTGTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCGCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt                                                              <  1:1634034/150‑1 (MQ=255)
                                                                                        gTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt                                                              <  1:493168/150‑1 (MQ=255)
                                                                                        gTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt                                                              <  1:1834792/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)
                                                                                                          gagTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCAcc                                            >  1:1598385/1‑150 (MQ=255)
                                                                                                           agTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCAccc                                           >  2:2006963/1‑150 (MQ=255)
                                                                                                            gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCTCCCg                                          <  2:1157617/150‑1 (MQ=255)
                                                                                                            gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCg                                          >  1:1440878/1‑150 (MQ=255)
                                                                                                            gTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCg                                          >  1:1798537/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)
                                                                                                                                                    |                                                                                                                                                     
TGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGTAATAATCACCTTATCACCCGGATAAGGGGTTGCCGGATCGACAATGTCTTTCGGTTTATA  >  NC_000913/4024230‑4024527

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: