Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,576,075 G→T 10.4% pseudogene (558/3567 nt) hsdR ← pseudogene, premature stop codon, host restriction; endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,576,0750GT10.4% 692.7 / 37.5 201pseudogene (558/3567 nt)hsdRpseudogene, premature stop codon, host restriction; endonuclease R
Reads supporting (aligned to +/- strand):  ref base G (98/82);  new base T (6/15);  total (104/97)
Fisher's exact test for biased strand distribution p-value = 3.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 203 total aligned reads displayed.

CTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGAC  >  CP009273/4575928‑4576220
                                                                                                                                                   |                                                                                                                                                 
cTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCg                                                                                                                                                   >  2:1403195/1‑148 (MQ=255)
cTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCg                                                                                                                                                   >  1:3192545/1‑148 (MQ=255)
  gTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGAc                                                                                                                                                 >  1:2081949/1‑148 (MQ=255)
  gTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGAc                                                                                                                                                 >  1:676085/1‑148 (MQ=255)
  gTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGAc                                                                                                                                                 <  1:1844513/148‑1 (MQ=255)
    tttGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACtt                                                                                                                                               <  2:3174652/148‑1 (MQ=255)
      tGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCt                                                                                                                                             >  1:1012852/1‑148 (MQ=255)
       gCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTg                                                                                                                                            <  1:1696295/148‑1 (MQ=255)
         cAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGcc                                                                                                                                          <  2:2157618/148‑1 (MQ=255)
            ttCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGa                                                                                                                                       >  1:2310896/1‑148 (MQ=255)
            ttCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGa                                                                                                                                       >  2:339896/1‑148 (MQ=255)
            ttCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGa                                                                                                                                       >  2:2638949/1‑148 (MQ=255)
              ccGCGTTCTTCTCGGCGAGCTGTGCGTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGt                                                                                                                                     <  1:505072/148‑1 (MQ=255)
               cgcgTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTc                                                                                                                                    <  1:1190009/148‑1 (MQ=255)
               cgcgTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTc                                                                                                                                    <  1:1190297/148‑1 (MQ=255)
                 cgTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTg                                                                                                                                  <  2:2310896/148‑1 (MQ=255)
                 cgTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTg                                                                                                                                  >  2:821816/1‑148 (MQ=255)
                    tcttctCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTgcgc                                                                                                                               >  2:2361028/1‑148 (MQ=255)
                    tcttctCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTgcgc                                                                                                                               >  1:2741456/1‑148 (MQ=255)
                    tcttctCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTgcgc                                                                                                                               >  1:581364/1‑148 (MQ=255)
                      ttctCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCtt                                                                                                                             <  2:164186/148‑1 (MQ=255)
                      ttctCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCtt                                                                                                                             <  1:433182/148‑1 (MQ=255)
                         tcGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCttttt                                                                                                                          <  1:1387289/148‑1 (MQ=255)
                         tcGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCttttt                                                                                                                          <  2:1066956/148‑1 (MQ=255)
                           ggCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTtct                                                                                                                        <  2:653716/148‑1 (MQ=255)
                              gAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGc                                                                                                                     <  1:531788/148‑1 (MQ=255)
                               aGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCa                                                                                                                    >  1:350444/1‑148 (MQ=255)
                                gCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCAc                                                                                                                   >  1:369096/1‑148 (MQ=255)
                                 cTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCAcc                                                                                                                  >  1:2622705/1‑148 (MQ=255)
                                 cTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCAcc                                                                                                                  <  1:1218931/148‑1 (MQ=255)
                                   gtgCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACctg                                                                                                                >  2:467868/1‑148 (MQ=255)
                                    tgCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACctgc                                                                                                               <  2:1465107/148‑1 (MQ=255)
                                           aGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCaa                                                                                                        >  2:1589114/1‑148 (MQ=255)
                                            gTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAg                                                                                                       <  1:1277614/148‑1 (MQ=255)
                                                cGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCtgt                                                                                                   <  1:3018965/148‑1 (MQ=255)
                                                   aaGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCtgttgt                                                                                                <  1:1739163/148‑1 (MQ=255)
                                                     ggCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGttt                                                                                              >  2:2621371/1‑148 (MQ=255)
                                                      gCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGtttt                                                                                             <  1:2845966/148‑1 (MQ=255)
                                                      gCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGtttt                                                                                             <  1:406843/148‑1 (MQ=255)
                                                            cTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTc                                                                                       <  1:2537845/148‑1 (MQ=255)
                                                               ggTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGc                                                                                    >  2:1740787/1‑148 (MQ=255)
                                                                 tttGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCAc                                                                                  <  1:1064303/148‑1 (MQ=255)
                                                                    gcgcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTc                                                                               <  2:147390/148‑1 (MQ=255)
                                                                    gcgcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTc                                                                               <  2:2501750/148‑1 (MQ=255)
                                                                     cgcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTcc                                                                              <  2:2985620/148‑1 (MQ=255)
                                                                     cgcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTcc                                                                              <  2:1545704/148‑1 (MQ=255)
                                                                     cgcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTcc                                                                              <  1:293450/148‑1 (MQ=255)
                                                                      gcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCt                                                                             <  1:1403195/148‑1 (MQ=255)
                                                                       cTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTg                                                                            <  1:2258588/148‑1 (MQ=255)
                                                                       cTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTg                                                                            <  2:789341/148‑1 (MQ=255)
                                                                        ttCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTgg                                                                           >  2:2852784/1‑148 (MQ=255)
                                                                        ttCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTgg                                                                           <  2:3274157/148‑1 (MQ=255)
                                                                         tCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGt                                                                          >  1:543072/1‑148 (MQ=255)
                                                                          cGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTg                                                                         <  2:1224677/148‑1 (MQ=255)
                                                                           ggTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGa                                                                        <  2:350444/148‑1 (MQ=255)
                                                                            gTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCa                                                                                                         <  2:2356786/114‑1 (MQ=255)
                                                                            gTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCa                                                                                                         >  1:2356786/1‑114 (MQ=255)
                                                                              ttCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATag                                                                     <  1:821816/148‑1 (MQ=255)
                                                                                ccTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATagag                                                                   >  1:3030467/1‑148 (MQ=255)
                                                                                  tgctgTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGt                                                                 <  1:2056795/148‑1 (MQ=255)
                                                                                  tgctgTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGt                                                                 >  2:1918474/1‑148 (MQ=255)
                                                                                   gctgTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGTGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGtt                                                                >  2:234954/1‑148 (MQ=255)
                                                                                     tgTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGtttt                                                              <  2:538227/148‑1 (MQ=255)
                                                                                       tttGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTca                                                            >  1:916974/1‑148 (MQ=255)
                                                                                       tttGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTca                                                            >  2:3081807/1‑148 (MQ=255)
                                                                                       tttGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTca                                                            >  1:2111480/1‑148 (MQ=255)
                                                                                           ccTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTcaccac                                                        >  1:1913437/1‑148 (MQ=255)
                                                                                           ccTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTcaccac                                                        >  1:2258530/1‑148 (MQ=255)
                                                                                           ccTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTcaccac                                                        >  1:1090114/1‑148 (MQ=255)
                                                                                               ccAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTc                                                    <  2:512800/148‑1 (MQ=255)
                                                                                                cAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTcc                                                   >  2:1260463/1‑148 (MQ=255)
                                                                                                  gAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCgg                                                 >  2:2035500/1‑148 (MQ=255)
                                                                                                  gAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCgg                                                 <  1:1773645/148‑1 (MQ=255)
                                                                                                   aaTGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGc                                                <  2:1564184/148‑1 (MQ=255)
                                                                                                    aTGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCa                                               >  1:1296808/1‑148 (MQ=255)
                                                                                                       gCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAaca                                            <  2:1077056/148‑1 (MQ=255)
                                                                                                          atatAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTTATAGAGGTTTTCACCACGTTCCGGCAACACaa                                         <  2:1025206/148‑1 (MQ=255)
                                                                                                           tataGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACaaa                                        >  2:751588/1‑148 (MQ=255)
                                                                                                           tataGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACaaa                                        >  1:62957/1‑148 (MQ=255)
                                                                                                             taGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAaca                                      <  1:810198/148‑1 (MQ=255)
                                                                                                              aGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAacac                                     >  1:1460453/1‑148 (MQ=255)
                                                                                                               gCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcacc                                    >  2:441884/1‑148 (MQ=255)
                                                                                                                ccGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccg                                   >  2:2668671/1‑148 (MQ=255)
                                                                                                                ccGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccg                                   >  2:2786689/1‑148 (MQ=255)
                                                                                                                 cGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccgg                                  >  1:2361269/1‑148 (MQ=255)
                                                                                                                   ttCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccggca                                <  2:2111480/148‑1 (MQ=255)
                                                                                                                   ttCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccggca                                >  1:840755/1‑148 (MQ=255)
                                                                                                                    tCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccggcac                               >  1:2076722/1‑148 (MQ=255)
                                                                                                                        ggCAACCAGCTTCTGCTGTTGCGCTTCTACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCggg                           <  1:3095308/148‑1 (MQ=255)
                                                                                                                          cAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGaa                         <  1:2361028/148‑1 (MQ=255)
                                                                                                                            aCCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGACGt                       >  1:521567/1‑148 (MQ=255)
                                                                                                                             ccAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTc                      <  2:2840421/148‑1 (MQ=255)
                                                                                                                             ccAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTc                      >  2:2768778/1‑148 (MQ=255)
                                                                                                                               aGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCAt                    >  2:1039074/1‑148 (MQ=255)
                                                                                                                                 cTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATaa                  >  1:1056255/1‑148 (MQ=255)
                                                                                                                                 cTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATaa                  >  1:2856203/1‑148 (MQ=255)
                                                                                                                                  ttCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAAt                 <  2:417502/148‑1 (MQ=255)
                                                                                                                                   tCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATc                >  2:3191325/1‑148 (MQ=255)
                                                                                                                                      gctgTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCttt             >  2:1392381/1‑148 (MQ=255)
                                                                                                                                                 tCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGAc  >  1:458987/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGAC  >  CP009273/4575928‑4576220

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