Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,759 G→C 5.6% intergenic (‑314/‑75) ylcI ← / → nohD DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7590GC 5.6% 934.9 / 2.9 302intergenic (‑314/‑75)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (149/120);  new base C (2/17);  total (152/150)
Fisher's exact test for biased strand distribution p-value = 1.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

TTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAG  >  NC_000913/580610‑580908
                                                                                                                                                     |                                                                                                                                                     
tttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGG                                                                                                                                                       >  2:3706562‑M1/1‑150 (MQ=255)
   gtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAA                                                                                                                                                    <  1:4788038‑M1/150‑1 (MQ=255)
    tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAAC                                                                                                                                                   <  1:244344‑M1/150‑1 (MQ=255)
      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGA                                                                                                                                                 >  1:2184985‑M1/1‑150 (MQ=255)
      ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGGTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGA                                                                                                                                                 >  2:3970919‑M1/1‑150 (MQ=255)
        tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACA                                                                                                                                               >  1:793452‑M1/1‑150 (MQ=255)
        tCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACA                                                                                                                                               >  2:167479‑M1/1‑150 (MQ=255)
           tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                                                            >  1:4049915/1‑150 (MQ=39)
                ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                                                       >  2:4597597/1‑150 (MQ=39)
                ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                                                       >  1:3062926/1‑150 (MQ=39)
                 gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCTTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaa                                                                                                                                      <  2:4950775/150‑1 (MQ=38)
                   gATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAg                                                                                                                                    >  2:1525210/1‑150 (MQ=38)
                      ccGACTGGTTAGGGGGCGGCGACCTCGCGGGTTTTCGCTTTTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTATTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGa                                                                                                                                 <  2:2390525/150‑1 (MQ=255)
                       cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                                                                >  2:1205250/1‑150 (MQ=39)
                       cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                                                                <  2:4288874/150‑1 (MQ=39)
                        gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                                                               >  2:5665968/1‑150 (MQ=39)
                         aCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCt                                                                                                                              >  2:3280634/1‑150 (MQ=39)
                             gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                                                          >  1:2042810/1‑150 (MQ=40)
                              ttACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                                                                         >  1:1806202/1‑150 (MQ=40)
                               tccGGGGCGGCGACCTGGCGGTTTTTCGCTATTTTTGAAATTTTTCGGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTgg                                                                                                                        <  2:5081826/148‑1 (MQ=14)
                                  ggggCGGCGCCCGGGGGGGTTTGGGATATTTCTGAAATTTGTCGGGTTTCCGGTGTTTCGGTTCTTGGGGGTGGTACGTTATTGTTTTTATTTAAAATACCCCCTGACAAGAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                                                                     <  2:2615225/150‑1 (MQ=2)
                                  ggggCGGCGACCTCGGGGGTTTCCGATATTTATTAAAATTTTCCGGTTTTTGGTGTTCCGGTTCTTCTTCGTCGTAACTTATTTTTTTTATTTAAAATACCCCCTCAAAAGAAAGGAAACGACAGaggccaaaaagctcgctttcagcac                                                                                                                     <  1:2884705‑M1/150‑26 (MQ=255)
                                     gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctg                                                                                                                  <  2:3136732‑M2/150‑4 (MQ=255)
                                      cggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgt                                                                                                                 <  2:2515149‑M2/150‑5 (MQ=255)
                                       ggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATAACCCCAGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtc                                                                                                                <  1:114765‑M2/150‑6 (MQ=255)
                                        gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTACGGTTTAAGGTTTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAAAACCCCAGAAAATAAAGGAAAAGAAAGGTGCTGAAAGCGAGCATTTTGGCCTctgtcg                                                                                                               <  1:1328364‑M2/150‑7 (MQ=255)
                                             cTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttcc                                                                                                          <  1:3048421‑M2/150‑12 (MQ=255)
                                               cgcgGGTTTTCGCTATTTATGACAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGaggccaaaaagctcgctttcagcacctgtcgtttcctt                                                                                                        <  2:4247641‑M1/150‑39 (MQ=255)
                                                gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGaggccaaaaagctcgctttcagcacctgtcgtttccttt                                                                                                       <  2:4636045‑M1/150‑40 (MQ=255)
                                                 cgGGATTTCGATTTTAATCAATATTTCCCGTTTTATGGTGTATCCTTTATTATGCGCCGTAATTAAATTTATTTTTTAAAAAAACCCCCTGAAAATAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtc                                                                                                      <  1:3882972/150‑1 (MQ=9)
                                                   ggTTTTCGCTATTTATGAAAATTTTCCGGTTGAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtctc                                                                                                    <  2:3294791/150‑1 (MQ=255)
                                                    gTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                                                   <  2:3392439/150‑1 (MQ=255)
                                                    gTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                                                   >  2:4213086/1‑150 (MQ=255)
                                                       ttCGCTATTTATGAAATTTTTCGGGTTTACGGTGTTTCCGTTCTTCTTCGTCGTAAGTTAATGTTTTTATTTAAAATACCCCCTGAAAAGAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtt                                                                                                <  2:2268471/150‑1 (MQ=11)
                                                        tCGCTATTTATTAAATTTTTCCGTTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAAATTAATTTATTTATTTAAAATACCCCCTTAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                                                               <  2:3024599/150‑1 (MQ=14)
                                                          gttttttATTAAATTTTTCGGTTTTATGTTTTTTCGTTTCTTCTGCTTATTAAATTTTTTTTTTTATTTAAAATACCCCCTCAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                                                             <  1:4435800/146‑1 (MQ=2)
                                                          gatATTTATGAAAATTTTCCGGTTTACGGTGTTTCCGTTCTTCTGCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                                                             <  2:139558/148‑1 (MQ=255)
                                                          gCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                                                             >  2:897038/1‑150 (MQ=255)
                                                          gCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                                                             >  1:735047/1‑150 (MQ=255)
                                                          gCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                                                             <  2:5246835/150‑1 (MQ=255)
                                                             aTTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTc                                                                                          >  1:3679996/1‑150 (MQ=255)
                                                                   gAAAATTTTCGGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                                                    >  2:1967254/1‑150 (MQ=255)
                                                                   gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                                                    >  1:37154/1‑150 (MQ=255)
                                                                     aatttttCCGTTTTCAGGTGTTTCGGTTCTTCTTCGTCGCAAATTAATGTATTTATTTAAAATACCCCCTTAAAATAAAGTACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAt                                                                                  <  1:4752223/147‑1 (MQ=14)
                                                                           tCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCGTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                                                            >  1:4251869/1‑150 (MQ=255)
                                                                             cGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                                                          >  2:5569067/1‑150 (MQ=255)
                                                                             cGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                                                          >  2:1277701/1‑150 (MQ=255)
                                                                               gTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                                                        >  2:4081946/1‑150 (MQ=255)
                                                                               gTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                                                        <  2:2843619/150‑1 (MQ=255)
                                                                               gTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                                                        >  2:801432/1‑150 (MQ=255)
                                                                                tttATGGTGTTCCTGTTTTTCTTCTTCTTAAATTATTTTATTTTTTTAAAATACCCCCTTAAATCAAACGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                                                       <  1:4912488/150‑1 (MQ=17)
                                                                                tttAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                                                       >  1:915346/1‑150 (MQ=255)
                                                                                 ttATGGTTTTTCCTTTCTTCTTCTTCTTAAATTATTTTATTTATTTAAAATACCCCCTTAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAg                                                                      <  2:1016668/150‑1 (MQ=14)
                                                                                  tAAGGTGTTTCGGTTCTTCTTCGTCGTAAATTAATTTATTTATTTAAAATACCCCCTGAAAACAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt                                                                     <  2:3460219/150‑1 (MQ=255)
                                                                                  tAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTGGTTTGTCTGTTTTTGTGGGTGGTATGTGGAATGGAAGt                                                                     >  1:5125076/1‑150 (MQ=17)
                                                                                      gtgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                                                 >  2:3247645/1‑150 (MQ=255)
                                                                                        gtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaacaa                                                               >  2:2506445/1‑150 (MQ=255)
                                                                                         tttCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTTCAATGGAAGTCAACaaa                                                              >  1:3613396/1‑150 (MQ=255)
                                                                                          ttCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                                                             >  1:3208512/1‑150 (MQ=255)
                                                                                          tacgtttcttcttcGTCTAAACTTATTGTTTTTTTTTAAAATCCCCCCTTAAAAGAATGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                                                             <  1:1634534/145‑1 (MQ=14)
                                                                                             cGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagc                                                          >  2:117213/1‑150 (MQ=255)
                                                                                              gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                                                         >  1:5438374/1‑150 (MQ=255)
                                                                                               ttcttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcag                                                        >  2:5461946/1‑150 (MQ=255)
                                                                                               ttcttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcag                                                        >  1:230152/1‑150 (MQ=255)
                                                                                                 cttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                                                      <  2:5805373/150‑1 (MQ=255)
                                                                                                  ttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctg                                                     >  2:624280/1‑150 (MQ=255)
                                                                                                   tcttcGTCGTAACTTAATGTATATATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCAGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctgg                                                    <  2:4834559/150‑1 (MQ=255)
                                                                                                     ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                                                  <  2:3200787/150‑1 (MQ=255)
                                                                                                         tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa                                              >  2:1589472/1‑150 (MQ=255)
                                                                                                           gtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGGTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTGCGTGGTATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtt                                            >  2:1778313/1‑150 (MQ=255)
                                                                                                               cTTAATTTATTTATTTAAAATACCCCCTTAAAACAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCg                                        <  2:3054040/150‑1 (MQ=255)
                                                                                                                ttAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGTATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGTCATTTTCgg                                       >  1:3966767/1‑150 (MQ=255)
                                                                                                                   ttttttttttttaATACACCCCCTTAAAATAATGTAAACTACTGGTGCTGAATGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc                                    <  1:2575449/145‑1 (MQ=14)
                                                                                                                   aTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc                                    <  2:4252713/150‑1 (MQ=255)
                                                                                                                       atttatttAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGt                                >  1:4507140/1‑150 (MQ=255)
                                                                                                                        tttatttaAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTAGGTGCGAGTa                               >  1:187308/1‑150 (MQ=255)
                                                                                                                           atttaAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATcc                            <  2:1887251/150‑1 (MQ=255)
                                                                                                                            tttaAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTGCGTGGCATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTAGGTGCGAGTATCCg                           >  1:105810/1‑150 (MQ=255)
                                                                                                                            tttaAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGCGTATCCg                           >  2:740438/1‑150 (MQ=255)
                                                                                                                              taAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTa                         >  2:4359727/1‑150 (MQ=255)
                                                                                                                              gcccagaCCCCCTGACACGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGCAGGTGCAATGGAAGTCAACAAACAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTa                         <  2:4321289/144‑1 (MQ=14)
                                                                                                                                    aCCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTc                   <  1:3384107/150‑1 (MQ=255)
                                                                                                                                    aCCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTc                   <  1:4343136/150‑1 (MQ=255)
                                                                                                                                      ccccTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAg                 >  2:2927022/1‑150 (MQ=255)
                                                                                                                                        ccTGAAAAGAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGaa               <  1:3469510/150‑1 (MQ=255)
                                                                                                                                         cTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAAc              >  2:3669996/1‑150 (MQ=255)
                                                                                                                                         cTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCATTCTCTGTTTTTTTCCGTGGAATGTGCAATGGCAGTCAACAAAAAACAGCTGGCGGACATTTTCGGTGCGAGTATCCGTACCAATCAGAAc              >  1:1523339/1‑150 (MQ=14)
                                                                                                                                          tGAAAAGACAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACt             <  2:4134832/150‑1 (MQ=255)
                                                                                                                                          tGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACt             <  1:3231000/150‑1 (MQ=255)
                                                                                                                                           gAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTg            <  2:4172970/150‑1 (MQ=255)
                                                                                                                                           gAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCGTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTg            >  2:1990951/1‑150 (MQ=255)
                                                                                                                                            aaacGAATGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCTATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTgg           <  1:2515013/150‑1 (MQ=255)
                                                                                                                                            aaaaGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTgg           <  1:210326/150‑1 (MQ=255)
                                                                                                                                              gataaagcaaACGACAGaggccaaaaagctcgctttcagcacctgtcgtttcctttcttttcagggtttattttaaataaatacattaagttacgacgaagaagaacggaaacaccttaaaccggaaaattttcataaatagcgaaaacc         <  1:1042136‑M1/147‑134 (MQ=255)
                                                                                                                                                gaaagGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAgg       <  1:1970266/150‑1 (MQ=255)
                                                                                                                                                 aaagGAAACGACAGaggccaaaaagctcgctttcagcacctgtcgtttcctttcttttcagggggtattttaaataaatacattaagttacgacgaagaagaacggaaacaccttaaaccggaaaattttcataaatagcgaaaacccgc      <  2:1970266‑M1/150‑137 (MQ=255)
                                                                                                                                                 aaagGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGa      >  2:3853929/1‑150 (MQ=255)
                                                                                                                                                  aagGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTctgtcgtttcctttcttttcagggggtattttaaataaatacattaagttacgacgaagaagaacggaaacaccttaaaccggaaaattttcataaatagcgaaaacccgcg     <  2:4247465‑M2/150‑113 (MQ=255)
                                                                                                                                                  aagGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGCGTATCCGTACCATTCAGAACAGGCAGGaa     >  2:2745460/1‑150 (MQ=255)
                                                                                                                                                   agGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAAc    >  2:2058722/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAG  >  NC_000913/580610‑580908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: