Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 576,183 G→T 10.1% intergenic (‑282/‑107) 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
*CP009273576,1830GT10.1% 569.6 / 10.2 186intergenic (‑282/‑107)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (104/53);  new base T (2/26);  total (106/80)
Fisher's exact test for biased strand distribution p-value = 2.60e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

GTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCG  >  CP009273/576035‑576331
                                                                                                                                                     |                                                                                                                                                    
gTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                      <  1:2298617/149‑1 (MQ=255)
gTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                      >  1:1070490/1‑149 (MQ=255)
gTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                      <  1:2057618/149‑1 (MQ=255)
gTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATg                                                                                                                                                      <  1:609641/149‑1 (MQ=255)
 ttGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGt                                                                                                                                                     <  2:4405924/149‑1 (MQ=255)
    tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatt                                                                                                                                                  <  1:2767433/149‑1 (MQ=255)
    tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatt                                                                                                                                                  <  1:5097031/149‑1 (MQ=255)
    tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatt                                                                                                                                                  <  2:2505599/149‑1 (MQ=255)
     gACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                 >  1:4365875/1‑149 (MQ=255)
       cATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                               <  1:1943097/149‑1 (MQ=255)
       cATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttat                                                                                                                                               <  1:2006916/149‑1 (MQ=255)
        aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                              >  2:1318764/1‑149 (MQ=255)
        aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                              >  2:2994136/1‑149 (MQ=255)
        aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                              <  1:3937660/149‑1 (MQ=255)
         tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                             >  2:2349064/1‑149 (MQ=255)
          ggaGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTAtttattta                                                                                                                                            <  1:5187795/147‑1 (MQ=255)
            aGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                          <  1:4851467/149‑1 (MQ=255)
             gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                         <  1:669430/149‑1 (MQ=255)
             gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                         >  2:1096145/1‑149 (MQ=255)
              atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                                        <  2:4746141/149‑1 (MQ=255)
              atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCGATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                                        >  2:3319201/1‑149 (MQ=255)
                gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                      >  1:1638711/1‑149 (MQ=255)
                 atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                     <  1:1312426/149‑1 (MQ=255)
                 atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                     <  1:4085933/149‑1 (MQ=255)
                 atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                     <  2:2385267/149‑1 (MQ=255)
                  tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                    <  2:1850477/149‑1 (MQ=255)
                  tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                    >  1:4572829/1‑149 (MQ=255)
                  tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                    >  2:2779906/1‑149 (MQ=255)
                    attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                  <  2:3799747/149‑1 (MQ=255)
                     ttattaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                 <  2:158756/149‑1 (MQ=255)
                     ttattaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                 <  2:2748749/149‑1 (MQ=255)
                      tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTg                                                                                                                                <  1:4394078/149‑1 (MQ=255)
                         tatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                             >  1:3027197/1‑149 (MQ=255)
                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaa                                                                                                                         >  2:3103239/1‑149 (MQ=255)
                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAaaagaa                                                                                                                         >  1:1062446/1‑149 (MQ=255)
                              ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                        >  1:1965976/1‑149 (MQ=255)
                               tttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                       <  1:1096145/149‑1 (MQ=255)
                                   tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                   <  2:4895126/149‑1 (MQ=40)
                                   tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                   <  1:1046240/149‑1 (MQ=40)
                                   tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                   <  2:747530/149‑1 (MQ=40)
                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACg                                                                                                                 >  2:5155129/1‑149 (MQ=39)
                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAAGGAAACg                                                                                                                 >  1:1816806/1‑149 (MQ=255)
                                      gTCCTTTCCGGCGATCCGACTGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                <  2:82345/149‑1 (MQ=255)
                                      gTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                >  2:2931882/1‑149 (MQ=39)
                                        ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACa                                                                                                              >  1:4723313/1‑149 (MQ=39)
                                          tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATGTATTTAAAATACCCCCTGAAAAGAAAGGACACGACAgg                                                                                                            >  2:3291181/1‑149 (MQ=255)
                                           ttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                           >  1:1220897/1‑149 (MQ=39)
                                             ccGGCGAGCCGACTGGTTCGGGGGGGGCGACCTGGGGGGTTTTCGATTTTTTTTAAAATTTTCGGTTTTATGGTGTTTCGGTTATTATTGGTTTTAATTTATTTTTTTTATTTAAAATACCCCCTTAAAAGAACGCAAACGACAGGTGc                                                                                                         <  2:2133019/149‑1 (MQ=2)
                                                gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                      >  2:1015600/1‑149 (MQ=38)
                                                gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                      >  2:3597286/1‑149 (MQ=38)
                                                gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                      >  2:1821342/1‑149 (MQ=38)
                                                    tCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                  >  1:2867010/1‑149 (MQ=38)
                                                     ccGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCg                                                                                                 >  2:4207503/1‑149 (MQ=39)
                                                       gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAg                                                                                               >  2:1140562/1‑149 (MQ=39)
                                                        aCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                              >  1:984842/1‑149 (MQ=39)
                                                        aCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                              >  2:2454863/1‑149 (MQ=39)
                                                        aCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                              >  2:5081159/1‑149 (MQ=39)
                                                           ggttcgggggCGGCGACCTCGGGGTTTTTCTCTTTTTATGAAATTTTGCCGTTTTATGGTGTTTCGTTTCTTCTTCGTCGAAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCttt                                                                                           <  2:2391495/143‑1 (MQ=14)
                                                           ggtcccGGGGCGGCGACCGCGCGGGTTTCCGATATTTATTAAAATTTCCCGGTTTATGGTGTTTCGGTTCTTCTTCGTCTAAAATTATTTTATTTATTTAAAATACCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                           <  2:5161507/144‑1 (MQ=255)
                                                            tttcgggggCGGCCCCCTCGTGTGTTTTCGATATTAATCAAATTTTTCCGTTTTAAGGTGTGTCTTTTCTTCTCCGTCTTAAATTAATTTTTTTATTAAAAAAACCCCCTAAAAACAAACCACACGACAGGTGCTGAAAGCGAGCtttt                                                                                          <  2:3487825/144‑1 (MQ=255)
                                                                   ggcggcGACGTCGGGGTTTTTCTATATTTATTAAATTTTGCGGGTTTTTGGTTTTTCGTTTCTTCTTCGTCTTAATTTATTTTATTTTTTTAAAAAACCCCCTTAAAATAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCctc                                                                                   <  2:1449888/149‑1 (MQ=255)
                                                                      ggcgACCTGGCGGTTTTTAGATTTTTTTTAAAATTTTCGGGTTTAAGGTTTTTCCGTTCTTTTTCGTGTTAACTTAATTTATTTATTTAAAATACCCCCGCAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGt                                                                                <  2:2322871/149‑1 (MQ=14)
                                                                       gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATATCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTc                                                                               >  2:4120023/1‑149 (MQ=255)
                                                                         gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGAAAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                             >  1:1029384/1‑149 (MQ=255)
                                                                         gACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTATTCTTCGTCGTAACTTAATTTATTTATTTAAAATAACCCCAGAAAAGAAAGGAAAAGACATGTGCTGAAAGCGAGCGTTTTGGCCTCTGTCGt                                                                             >  1:3494212/1‑149 (MQ=14)
                                                                          aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATAACCCATGAAAAGAAAGGAAACTACAGGTGCATAAAGCAAGATTTTTGGCCACTGTAGtg                                                                            >  2:2605254/1‑148 (MQ=14)
                                                                          aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                            >  1:3687445/1‑149 (MQ=255)
                                                                          aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                            >  2:3447342/1‑149 (MQ=255)
                                                                           ccTCGGGGTTTTTCGATATTTATTAAAATTTTCCGTTTTATGGTGTTTCCTTTATTCTTCTTCTTAAATTAATTTATTTATTAAAAATACCCCCTTAAAACAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGttt                                                                           <  2:4948207/149‑1 (MQ=255)
                                                                           accccggggTTTTCGGCTTTTTATGAAAATTTTCCGTTTTATGGTTTTTCCGTACTTCTCCGTCTTAACTTATTGTATTTATTTAAAATACCCCCTCAAAATAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                            <  2:722280/143‑1 (MQ=255)
                                                                            cgggggggTTTGCGCTTTTTTTGAAATTTGGCGGGTTTATGGTGTTGCCGTTGTTCTTCGTCGTAACTTTATTTTTTTTTTTAAAACCCCCCGTCAAACGACAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTc                                                                          <  2:4555618/144‑1 (MQ=255)
                                                                            cTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTc                                                                          >  2:4186556/1‑149 (MQ=255)
                                                                             tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTcc                                                                         >  1:1328439/1‑149 (MQ=255)
                                                                                 gggTTTTCGTTTTTTATTAATATTTTCCGGTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTTTTTATTTAAAATCCCCCCTGAAAAGAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtc                                                                     <  2:4163345/149‑1 (MQ=11)
                                                                                  ggTTTTCGCTCTTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTGTCCTTtct                                                                    >  2:3504719/1‑149 (MQ=255)
                                                                                  ggTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                    >  2:172855/1‑149 (MQ=255)
                                                                                           cTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttttt                                                            >  2:3813320/1‑149 (MQ=255)
                                                                                                 aTGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTGCTTTATCTGTTTTTTTCCgg                                                      >  1:1440538/1‑148 (MQ=255)
                                                                                                  tGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTg                                                     >  2:2082239/1‑149 (MQ=255)
                                                                                                  tGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTg                                                     >  2:679780/1‑149 (MQ=255)
                                                                                                    aaaaTTTTCGGGTTTATGTTTTTTCGTTTATTCTTCTTCTTAACTTAATTTATTTATTTAAAATCCCCCCTGAAAATAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                   <  2:4928791/149‑1 (MQ=14)
                                                                                                             cGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTCTCTGTTTTTGTCCGTGGAATTTGCAAt                                          >  2:4039005/1‑149 (MQ=255)
                                                                                                              gttttATGGTGTTTCCGTTCTTCTTCGTCTTAAATTAATTTTTTTATTTAAAATACCCCATTAAAACAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                         <  2:1991321/147‑1 (MQ=14)
                                                                                                                    aGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTc                                   >  2:867137/1‑149 (MQ=255)
                                                                                                                       tgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaac                                >  2:3173186/1‑149 (MQ=255)
                                                                                                                        gtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTGTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaaca                               >  2:4600096/1‑149 (MQ=255)
                                                                                                                             cGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTGAGTTGTATTTGGAATTGAATTCaaaaaaaaa                          >  1:4053424/1‑148 (MQ=14)
                                                                                                                             cGTTCTTCGTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTACCTGGAATGTGCAATGGAAGTCAACAAAAag                          >  2:4609330/1‑149 (MQ=255)
                                                                                                                                  tacttcGTCTTAAATTAATGATTTTATTTAAAATACCCCCTTAAAACAAACTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                     <  2:3077849/147‑1 (MQ=14)
                                                                                                                                     ttcttcTTAAATTATTTTTTTTATAAAAAAAACCCCAACAAAACAAACCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggc                  <  1:1187768/149‑1 (MQ=14)
                                                                                                                                     ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggc                  >  2:2373342/1‑149 (MQ=255)
                                                                                                                                      tcgtcgTAAATTAATTTATTTATTAAAAAAACCCCCTGAAAAGAATGTACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                 <  2:4606658/149‑1 (MQ=14)
                                                                                                                                       cgtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTTATCTGTTTTTTTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgcttggtg                >  1:5176978/1‑144 (MQ=255)
                                                                                                                                         tcttAAATTAATGTATTTATTTAAAAAACCCCCTAAAAATACAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGAc              <  2:958443/146‑1 (MQ=14)
                                                                                                                                         tataaaattAATTTATTTAAAAAAAAACCCCCATAAAAAAAAACAACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGAc              <  1:905813/142‑1 (MQ=14)
                                                                                                                                             aaCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGTAGTCAACAAAAAGCAGCTGGCTGACAttt          >  2:392830/1‑149 (MQ=255)
                                                                                                                                              aCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCCGCTGGCTGACAtttt         >  1:5178942/1‑149 (MQ=255)
                                                                                                                                              aCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCAGGCTGACCtttt         >  1:5179020/1‑149 (MQ=255)
                                                                                                                                                 tAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCgg      >  1:2497572/1‑149 (MQ=255)
                                                                                                                                                  attGTTTTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt     <  2:1965300/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCG  >  CP009273/576035‑576331

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: