Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,575,840 A→C 9.1% pseudogene (793/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,575,8400AC 9.1% 631.8 / 8.3 229pseudogene (793/3567 nt)hsdRpseudogene, premature stop codon, host restriction; endonuclease R
Reads supporting (aligned to +/- strand):  ref base A (122/79);  new base C (3/25);  total (125/104)
Fisher's exact test for biased strand distribution p-value = 6.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

AAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCC  >  CP009273/4575692‑4575988
                                                                                                                                                    |                                                                                                                                                    
aaGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGTCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTGTGGAGAATCTCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCCGGaa                                                                                                                                                      >  2:293200/1‑149 (MQ=255)
aaGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGaa                                                                                                                                                      >  2:599942/1‑149 (MQ=255)
  gCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGc                                                                                                                                                    <  1:3269718/149‑1 (MQ=255)
    ccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGa                                                                                                                                                  <  1:4205711/149‑1 (MQ=255)
     cTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGAc                                                                                                                                                 >  1:1728045/1‑149 (MQ=255)
      tGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGAct                                                                                                                                                >  1:2430745/1‑149 (MQ=255)
       gATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGActc                                                                                                                                               >  2:1754166/1‑149 (MQ=255)
       gATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGActc                                                                                                                                               >  2:4058130/1‑149 (MQ=255)
           aCCCGTTTCATCTTTTCCGGTCGGACATTCGGCAATGGCTATATAGACGCCGGGTACCGGACGAGCGCCTTTGGAGAAGCGCAGCGTTTTGCTGTCTGCCTGCCAGCCTGCTTTACGCAGTAGCGCATCAATCAGGAAGCGACTCtctt                                                                                                                                           >  2:361288/1‑149 (MQ=255)
             ccGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCtcttct                                                                                                                                         >  2:529618/1‑149 (MQ=255)
               ttttCATATTTTCCGGTCGGCCATTCGGAAATGGTTTTATTGACGCCGGTTTCCGGACGTGCGCATTTGGAGAAGCGCTGGTTTTTGCTGTCGGCCTGCAAGCCTGATTTACGATTTTGCGAATCAATCAGGAAGCGACTCTCttcttc                                                                                                                                       <  2:997630/148‑1 (MQ=255)
               gTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCttcttc                                                                                                                                       >  1:3318154/1‑149 (MQ=255)
               gTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCttcttc                                                                                                                                       >  2:2730720/1‑149 (MQ=255)
                  tCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCt                                                                                                                                    >  1:2346365/1‑149 (MQ=255)
                    aTCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTaa                                                                                                                                  >  1:530463/1‑149 (MQ=255)
                    aTCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTaa                                                                                                                                  >  1:530634/1‑149 (MQ=255)
                     tCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAg                                                                                                                                 >  1:2133769/1‑149 (MQ=255)
                     tCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAg                                                                                                                                 >  2:3800403/1‑149 (MQ=255)
                      cTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAgg                                                                                                                                >  1:1480029/1‑149 (MQ=255)
                       ttttCCGGTGGGCAATTCGGAATTGGCTTTTTTGAGGCGGGGTTCGGGACGTGCGCCTTTGGATACGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGATTTACGCTTTTGCGAATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGt                                                                                                                               <  2:875893/149‑1 (MQ=255)
                             ggTCGGCCATTCGGCATTGGCTTTATGGCCGCCGGGTTCCGGACGTGCGCCTTGGGAGAAGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGATTTACGCTGTTGCTCATAAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAgtg                                                                                                                         <  1:3358624/149‑1 (MQ=255)
                                  gCCATTCGGCAATGGCTTTTTTGACGCGGGGTTCCGGACGTGCGCCTTTGGAGACGCGCTGGTTTTGGCTGTGGGCCTGCCAGCCTGCTTTACGCTGTTGCGCATCAATCTGCACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTgcgc                                                                                                                    <  1:2531649/149‑1 (MQ=255)
                                   ccATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCt                                                                                                                   >  2:1545285/1‑149 (MQ=255)
                                    cATTCGGCATTGGTTTTATGGACGCGGGGTGCCGGACGTGCGCCTTTGGAGACGCGCTGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTCCGCTGTTGCGAATCAATCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCtt                                                                                                                  <  2:3896399/149‑1 (MQ=255)
                                    cATTCGGCATGGGATTTTTTGACGCCGGGTGCGGGCCGTGCGCGTTTGGAGACGCGCTGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTCCGCTGTTGCGCATCAATCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCtt                                                                                                                  <  2:1395214/149‑1 (MQ=255)
                                    cATTCGGCAATGGATTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGATAAGCGCAGGGTTTTGCTTTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCtt                                                                                                                  <  2:2260308/149‑1 (MQ=255)
                                     aTTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTg                                                                                                                 >  2:2491978/1‑149 (MQ=255)
                                      ttCGGCAATGGCTTTATTGCCGCCGGGTTCGGGAGGTGCGCCTTTGGAGACGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCTGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGa                                                                                                                <  2:2779715/149‑1 (MQ=255)
                                      ttCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGa                                                                                                                >  1:1210677/1‑149 (MQ=255)
                                      ttCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGGTGAGTGTGCGCTTGa                                                                                                                >  2:3033217/1‑149 (MQ=255)
                                       gcggaattGGATTTTTTGCCGCGGGGTTCCGGCCGTGCGCCTTGGTTGACGCGCCGGTTTTTGCTGTCGGCCGGCCCGCCTGCTTTACGCTGTTGCGAAACAATCAGCAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGAt                                                                                                               <  1:1069511/142‑1 (MQ=255)
                                          gCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGc                                                                                                            >  1:4111396/1‑149 (MQ=255)
                                          gCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGc                                                                                                            >  2:1164933/1‑149 (MQ=255)
                                           cAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGcc                                                                                                           >  2:660972/1‑149 (MQ=255)
                                            attGGGTTTTTTGCCGCGGGGTTCGGGCGGTGCGCGTTGGGAGACGCGCTGGGTTTGGCTGGCGGCCGGCCCGCCTGCTTTACGCTGGGGCGCACCACTCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCt                                                                                                          <  1:3710072/147‑1 (MQ=255)
                                            aaTGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCt                                                                                                          >  1:1015895/1‑149 (MQ=255)
                                             ttGGCTTTATGGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCTGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCTGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTg                                                                                                         <  2:3580059/148‑1 (MQ=255)
                                                gCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTCAGGTTGAGTGTGCGCTTGATGGCCTGATc                                                                                                      >  1:3423516/1‑149 (MQ=255)
                                                  tttATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCgg                                                                                                    >  2:226449/1‑149 (MQ=255)
                                                     aTTGCCGCGGGGTTCGGGACGTGCGCCTTTGGAGACGCGCTGGGTTTTGCTGGCGGCCTGCCAGCCTGCTTTACGCTGTGGCGCATCAAACAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTaa                                                                                                 <  2:459445/149‑1 (MQ=255)
                                                      ttGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCCGGTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAAt                                                                                                >  1:3963298/1‑149 (MQ=255)
                                                      ttGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGACCGGTAAt                                                                                                >  2:3966936/1‑149 (MQ=255)
                                                       tgccgcgGGGTTCCGGCCGTGCGCTTGTGGAGAAGCGCTGGGTTTGGGTGTCGGCCTGCCAGCCTGTTTTACGCAGTGGCGAATAAAACAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAAtt                                                                                               <  1:3927096/146‑1 (MQ=255)
                                                        gccGCCGGGTTCCGGACGTGCGCCTTGGGAGAAGCGCTGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTGGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAAttt                                                                                              <  2:4400190/147‑1 (MQ=255)
                                                         aCGCCGGGTTCCGGACGTGCGCCTTTGGAGATGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTc                                                                                             <  2:2635904/149‑1 (MQ=255)
                                                           gcggggTTCGGGCGGTGCGCCTGTGGATATGCGATGGGTTTTGCTGTGGGCCTGCCAGCGTGATTTAGGCTGTTGCTAATAAATCTGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCtt                                                                                           <  1:1572673/146‑1 (MQ=255)
                                                           gCCGGGTTCCGGACGTGCGCCTTGGGAGATGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGAATAAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCtt                                                                                           <  2:3521406/149‑1 (MQ=255)
                                                               ggTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTtgtg                                                                                       >  1:1417984/1‑149 (MQ=255)
                                                                 ttCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGt                                                                                     >  2:3039727/1‑149 (MQ=255)
                                                                      gACGTGCGCCTTTGTAGACGCGCTGGGTTTTGCGGTCGGCCTGCCCGCCTGCTTTACGCTGTTGCGAATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGcc                                                                                <  1:4025831/149‑1 (MQ=255)
                                                                      gACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGACGCGACTCTCTTCTTCGCTCAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGcc                                                                                >  2:1099959/1‑149 (MQ=255)
                                                                      gACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGcc                                                                                >  2:968755/1‑149 (MQ=255)
                                                                      gACGTGCGCATTTGTAGATGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTTTTTTACGATGTTGCGCAAAAAAACGCACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGcc                                                                                <  1:3569557/149‑1 (MQ=255)
                                                                        ggTGCGCCTTGGGAGAAGCGCTGGTTTTTGGTGTCGGCCTGCCAGCCTGCTTTCCGCTGTTGCGCATCAATCCGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCtt                                                                              <  1:1683974/148‑1 (MQ=255)
                                                                         gTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTa                                                                             >  2:3702963/1‑149 (MQ=255)
                                                                           gcgcCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACg                                                                           >  2:4015984/1‑149 (MQ=255)
                                                                             gcCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGtt                                                                         >  2:4350731/1‑149 (MQ=255)
                                                                              ccTGTGTAGACGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTCCGCTGTTGCGCATAAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTc                                                                        <  2:1189593/149‑1 (MQ=255)
                                                                                tttGGTGATGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGAATAAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCt                                                                      <  2:1219844/149‑1 (MQ=255)
                                                                                tttGGAGAAGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCTGTTGCGCATCAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCt                                                                      <  2:3439647/149‑1 (MQ=255)
                                                                                   ggAGATGCGCTGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTCCGCTGTGGCGCATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGtt                                                                   <  2:155404/149‑1 (MQ=255)
                                                                                    gttacgcgcTGGGTTTGGCTGTCGGCCTGCCAGCCTGTTTTACGCTGTGGCGCATCAAACAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTc                                                                  <  1:3185697/144‑1 (MQ=255)
                                                                                      tacgcgcTGTGTTTTGCGGGCGGCCTGCCAGCCTGGTTTACCCAGTTGCTCAAAAATCAGCAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCgg                                                                <  1:3909331/146‑1 (MQ=255)
                                                                                         gcgcAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCt                                                             >  2:4309592/1‑149 (MQ=255)
                                                                                         gcgcAGGGTTGTGCTGTCGGCCTGCCAGCCTGATTTACGCTGTTGCGCATAAAGCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCt                                                             <  2:3960785/149‑1 (MQ=255)
                                                                                             tgggTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGttt                                                         <  2:4473507/148‑1 (MQ=255)
                                                                                              gggTTTTGCTGTCGGCCTGCCAGCCGGATTTACGCTGTTGCGCATAAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGtttt                                                        <  1:272732/149‑1 (MQ=255)
                                                                                              gggTTTGGCGGCCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGtttt                                                        <  2:2556379/149‑1 (MQ=255)
                                                                                                 ttttGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTCATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGcc                                                     >  1:2052740/1‑149 (MQ=255)
                                                                                                 ttttGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGcc                                                     >  1:690096/1‑149 (MQ=255)
                                                                                                 ttttGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACATTCCTGTTCGGTCTGTTTTGcc                                                     >  1:4355805/1‑149 (MQ=255)
                                                                                                     gCTGTCGGCCTGCCAGCCTGCTTTACGCTGTGGCGCATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGtt                                                 <  1:2488983/149‑1 (MQ=255)
                                                                                                           ggCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGt                                           <  1:4565469/149‑1 (MQ=255)
                                                                                                            gCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGtt                                          <  2:1319896/149‑1 (MQ=255)
                                                                                                            gCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGtt                                          >  2:2315744/1‑149 (MQ=255)
                                                                                                            cccggCCAGCCTGCTTTCCGCTGTTGCGCAACAACCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGtt                                          <  1:954125/145‑1 (MQ=255)
                                                                                                                gCCAGCCTGCTTTCCGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttcttc                                      <  2:1654732/149‑1 (MQ=255)
                                                                                                                gCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttcttc                                      <  2:2042894/149‑1 (MQ=255)
                                                                                                                 cccgCCTGATTTACGCTGTGGCCCATAAATCCGCACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTtcttct                                     <  1:912876/146‑1 (MQ=255)
                                                                                                                   aGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCg                                   >  1:2708466/1‑149 (MQ=255)
                                                                                                                   aGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCg                                   >  2:2333867/1‑149 (MQ=255)
                                                                                                                    gCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCgg                                  >  2:3316566/1‑149 (MQ=255)
                                                                                                                        gCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAg                              >  1:2837083/1‑149 (MQ=255)
                                                                                                                        gCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAg                              >  1:3878300/1‑149 (MQ=255)
                                                                                                                         cTTTACGCCGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGc                             >  2:2024042/1‑149 (MQ=255)
                                                                                                                             aCGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCtgtg                         >  1:3551542/1‑149 (MQ=255)
                                                                                                                              cGCTGTGGCGCATCCATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGc                        <  2:2420648/149‑1 (MQ=255)
                                                                                                                              cGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGc                        <  2:475887/149‑1 (MQ=255)
                                                                                                                               gCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCt                       >  2:2353687/1‑149 (MQ=255)
                                                                                                                                  gTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCGAAGGGTGAGTGTGCGCTTGATGGCCTGATCGGTCATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTcc                    >  2:1888696/1‑149 (MQ=255)
                                                                                                                                         atcaatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCg             >  1:3494908/1‑149 (MQ=255)
                                                                                                                                         atcaatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCg             >  1:3496026/1‑149 (MQ=255)
                                                                                                                                          tcaatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGc            >  1:1906874/1‑149 (MQ=255)
                                                                                                                                          tcaatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGc            >  1:3364305/1‑149 (MQ=255)
                                                                                                                                            aatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCaa          >  2:1678655/1‑149 (MQ=255)
                                                                                                                                            aatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCaa          <  1:457408/149‑1 (MQ=255)
                                                                                                                                                 ggACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGa     <  2:53501/149‑1 (MQ=255)
                                                                                                                                                 ggAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTACTTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGa     >  2:3681379/1‑149 (MQ=255)
                                                                                                                                                 ggAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGa     >  1:818168/1‑149 (MQ=255)
                                                                                                                                                    aGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGcc  <  1:4466298/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCC  >  CP009273/4575692‑4575988

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: