Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,575,840 A→C 5.6% 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 5.6% 767.2 / 3.0 262pseudogene (793/3567 nt)hsdRpseudogene, premature stop codon, host restriction; endonuclease R
Reads supporting (aligned to +/- strand):  ref base A (146/94);  new base C (2/20);  total (148/114)
Fisher's exact test for biased strand distribution p-value = 3.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

AAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCC  >  CP009273/4575692‑4575988
                                                                                                                                                     |                                                                                                                                                    
aaGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGaa                                                                                                                                                      >  2:2151825/1‑149 (MQ=255)
  gCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGGTGCGCATCAATCAGGAAGc                                                                                                                                                    >  1:115556/1‑149 (MQ=255)
   cccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCg                                                                                                                                                   <  1:1502933/149‑1 (MQ=255)
   cccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCg                                                                                                                                                   <  1:4519927/149‑1 (MQ=255)
   cccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCg                                                                                                                                                   <  1:1907985/149‑1 (MQ=255)
   cccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCGTTATTGACGCCGGGTTCCGGACGTGCGCCGTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCg                                                                                                                                                   >  2:261640/1‑149 (MQ=255)
    ccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCAGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGa                                                                                                                                                  >  1:4036618/1‑149 (MQ=255)
     cTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGAc                                                                                                                                                 >  2:3048322/1‑149 (MQ=255)
        aTTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGActct                                                                                                                                              >  1:5119838/1‑149 (MQ=255)
         tttcccGTTTGATTTTTTGCGGGTCGGCCTTGGGGAATTGTGTTTTTTGCGGCTGGTTTCGGGACGGGCGCTTTTGGAGATGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTTCTTTAGGATTTTGCGATTAAATCAGGAAGCGActct                                                                                                                                              <  2:3522970/146‑1 (MQ=255)
           aCCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCtctt                                                                                                                                           >  2:4085656/1‑149 (MQ=255)
             ccGTTTCATCTTTTCCGGTGGGCCTTTGGGAATTGGCTTTTTTGACGCGGGGTTCGGGAGGTGCGCCTTTGGTGAAGCGCTGGGTTTTGCTGTCGGCCTGCCTGCCTGTTTTACGCTGTTGCGCTTCAATCTGGAAGCGACTCtcttct                                                                                                                                         <  2:221860/149‑1 (MQ=255)
              cGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCCGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTcttctt                                                                                                                                        >  2:2939699/1‑149 (MQ=255)
              cGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTcttctt                                                                                                                                        >  1:3890464/1‑149 (MQ=255)
                 ttCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGc                                                                                                                                     >  1:2391750/1‑149 (MQ=255)
                  tCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCt                                                                                                                                    >  2:4538403/1‑149 (MQ=255)
                  tCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCt                                                                                                                                    >  2:4525802/1‑149 (MQ=255)
                   cATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTa                                                                                                                                   >  1:538697/1‑149 (MQ=255)
                   cATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTa                                                                                                                                   >  1:2123497/1‑149 (MQ=255)
                    aTCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTaa                                                                                                                                  >  2:226823/1‑149 (MQ=255)
                     tCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAg                                                                                                                                 >  2:4320791/1‑149 (MQ=255)
                     tCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGGTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAg                                                                                                                                 >  2:1705316/1‑149 (MQ=255)
                       ttttCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGt                                                                                                                               >  2:4676135/1‑149 (MQ=255)
                              ggTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAgtg                                                                                                                         >  2:1421561/1‑149 (MQ=255)
                                    ccATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCt                                                                                                                   >  2:3346427/1‑149 (MQ=255)
                                       ttCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTCAGGTTGAGTGTGCGCTTGa                                                                                                                >  2:5006581/1‑149 (MQ=255)
                                       ttCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGa                                                                                                                >  1:2615650/1‑149 (MQ=255)
                                         cGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATg                                                                                                              >  2:4776999/1‑149 (MQ=255)
                                           gaaaTGGATTTTTTGACGCGGGGTTCCGGACGTGCGCCTTTGTAGATGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGTTTTACGCAGTTGCGAATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGc                                                                                                            <  1:3838301/147‑1 (MQ=255)
                                             aaTGGCTTTATTGCCGCCGGGTTCGGGCCGTGCGCTTTTGTATAAGCGCTGGTTTTTGCTGTCGGCCTGCCACCCTGCTTTACGCAGTTGCGAAACAATCTGCAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCt                                                                                                          <  1:1616422/149‑1 (MQ=255)
                                              aTGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTg                                                                                                         >  1:2549600/1‑149 (MQ=255)
                                                ggCTTTATTGCCGCCGGGTTCCGGACGGGCGCCTTTGGAGAAGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCTGTTGCGCATCAATCCGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGAt                                                                                                       <  2:1394681/149‑1 (MQ=255)
                                                ggCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGAt                                                                                                       >  2:2813455/1‑149 (MQ=255)
                                                 gattttttGACGCCGGGTTCCGGAGGTGCGCCTTTGGATAAGCGCAGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATc                                                                                                      <  1:3359565/147‑1 (MQ=255)
                                                 gCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATc                                                                                                      >  2:1289789/1‑149 (MQ=255)
                                                     tATTGCCGCGGGGTTCGGGACGGGCGCCTTTGGATATGCGCTGGTTTTTGCTGGCGGCCTGCCAGCCTGCTTAACGCTGTTGCGAATAAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTa                                                                                                  <  2:4366007/149‑1 (MQ=255)
                                                     tATTGCCGCCGGGTGCCGGCCGTGCGCCTTTGGAGAAGCGCTGGTTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCTGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTa                                                                                                  <  2:4364469/149‑1 (MQ=255)
                                                      aTTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTaa                                                                                                 >  2:4478543/1‑149 (MQ=255)
                                                      aTTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTaa                                                                                                 >  2:3959558/1‑149 (MQ=255)
                                                        tGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAAtt                                                                                               >  1:4907235/1‑149 (MQ=255)
                                                           cGCCGGGGTCCGGACGTGCGCCTTGGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTAGGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCt                                                                                            <  2:3838449/149‑1 (MQ=255)
                                                                 gTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTgg                                                                                      >  1:3169915/1‑149 (MQ=255)
                                                                   tCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCTGCCTGCTTTACGCTGTTGCGAATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTa                                                                                    <  2:1060415/149‑1 (MQ=255)
                                                                    ccGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTaa                                                                                   >  1:2457088/1‑149 (MQ=255)
                                                                      ggACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCGAAGGTTGAGTGTGCGCTTGATGGCCTGTGCGGTAATTTCTTTGTGGTAAGc                                                                                 >  2:2495228/1‑149 (MQ=255)
                                                                        aCGGGCGCCTTTGTTGAAGCGCTGGTTTTTGCTGTCGCCCGGCCAGCCTGCTTTACGCTGTTGCGAAACAACCAGAAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCt                                                                               <  2:896290/149‑1 (MQ=255)
                                                                          gTGCGCCTTTGGAGACGCGCTGGTTTTTGCTGTCGGCCTGCCCGCCTGCTTTACGATGTTGCGCATCAATCTGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTa                                                                             <  2:24591/149‑1 (MQ=255)
                                                                           tGCGCCTTTGGTGAAGCGCTGGGTTTTGGTGTCGGCCTGCCAGCCTGATTTACGCTGTTGCGAATAAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTAc                                                                            <  2:1016423/149‑1 (MQ=255)
                                                                                   tGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGt                                                                    >  1:2174267/1‑149 (MQ=255)
                                                                                     gagaAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTc                                                                  >  2:2345675/1‑149 (MQ=255)
                                                                                     gagaAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTc                                                                  >  2:4953866/1‑149 (MQ=255)
                                                                                       gaAGCGCAGGGTTTTGCGGTCGGCCTGCCAGCCTGCTTTCCGCAGTTGCGAATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCgg                                                                <  2:2540870/149‑1 (MQ=255)
                                                                                         cgcgcAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTc                                                              <  2:4319099/148‑1 (MQ=255)
                                                                                         aGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTc                                                              >  2:720336/1‑149 (MQ=255)
                                                                                          gcgcTGGTTTTTGCTGTCGGCCGCCCAGCCTGCTTTACGCTGTTGCGAAAAAACAAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCt                                                             <  1:1411027/149‑1 (MQ=255)
                                                                                            tctggtttttGGTGTGGGCCTGCCAGCGTTTTTTACGCTGTTGCGAATAACTCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGt                                                           <  1:4408098/143‑1 (MQ=255)
                                                                                             ctggtttttGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATAAATCAGGAAGCGACTCTCTTCTTCGCTTAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGtt                                                          <  2:4002276/144‑1 (MQ=255)
                                                                                             ctggtttttGCTGTCGGCCTGCAAGCCTGCTTTACGCAGTTGCGCATCAATCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGtt                                                          <  2:2075037/144‑1 (MQ=255)
                                                                                               gggTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGtttt                                                        >  2:4477782/1‑149 (MQ=255)
                                                                                               gggTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGtttt                                                        >  2:5106531/1‑149 (MQ=255)
                                                                                                   tttGCTGTCGGCCTGCCAGCCTGATTTACGCTGTTGCGAATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCa                                                    <  1:5081849/149‑1 (MQ=255)
                                                                                                    ttGGTGTCGGCCTGCCCGCCTGTTTTACGCTGGTGCGAACCAATCAGCACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAg                                                   <  1:4529003/149‑1 (MQ=255)
                                                                                                     tGCTTTCGGCCTGCCTGCCTTCTTTTGGATTTTGCTAATCAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGt                                                  <  1:4684803/149‑1 (MQ=255)
                                                                                                            ggCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTACTTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGt                                           >  1:1351806/1‑149 (MQ=255)
                                                                                                             gCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGtt                                          >  2:114720/1‑149 (MQ=255)
                                                                                                              ccTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttc                                         <  2:4409538/149‑1 (MQ=255)
                                                                                                              ccTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttc                                         <  2:1499528/149‑1 (MQ=255)
                                                                                                              ccTGCCAGCCTGATTTACGCTGTTGCGAATCAATCAGCACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttc                                         <  1:2212859/149‑1 (MQ=255)
                                                                                                               cTGCCAGCCTGATTTACGCTGTTGCGCATCAATCCGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttct                                        <  1:2728947/149‑1 (MQ=255)
                                                                                                                tGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttctt                                       >  2:592815/1‑149 (MQ=255)
                                                                                                                tGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttctt                                       >  1:294237/1‑149 (MQ=255)
                                                                                                                 gCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttcttc                                      <  2:5136841/149‑1 (MQ=255)
                                                                                                                 gCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttcttc                                      >  1:5195533/1‑149 (MQ=255)
                                                                                                                  ccAGCGTGCTTTACGCTGTTGCGCATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTtcttct                                     <  1:1983180/149‑1 (MQ=255)
                                                                                                                  ccAGCCTGCTTTACGCAGTTGCGCAAAAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTtcttct                                     <  2:422597/149‑1 (MQ=255)
                                                                                                                   cAGCGTGATTTACGCTGTTGCGAATCAATCAGGACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTtctc                                    <  1:4064569/149‑1 (MQ=255)
                                                                                                                   cAGCCTGCTTTACGCAGTTGCGCATCACTCTGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTtctc                                    <  1:1249425/149‑1 (MQ=255)
                                                                                                                    aGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCg                                   >  2:557416/1‑149 (MQ=255)
                                                                                                                     gCCTGTTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCgg                                  <  1:2904746/149‑1 (MQ=255)
                                                                                                                      ccTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGc                                 >  2:3659458/1‑149 (MQ=255)
                                                                                                                       cTGCTTTAGGCTGTGGCGCATCAATCCGCACGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCg                                <  1:4891293/149‑1 (MQ=255)
                                                                                                                       cTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCg                                >  2:473677/1‑149 (MQ=255)
                                                                                                                           tttCCGCAGTTGCGAATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCt                            <  1:864266/149‑1 (MQ=255)
                                                                                                                              aCGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCtgtg                         >  1:3574658/1‑149 (MQ=255)
                                                                                                                                      gcgcATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGt                 >  1:930160/1‑149 (MQ=255)
                                                                                                                                        gcACCAATCAGTAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGc               <  1:5006581/149‑1 (MQ=255)
                                                                                                                                         cATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGcc              >  1:4692443/1‑149 (MQ=255)
                                                                                                                                         cATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGcc              >  2:2681838/1‑149 (MQ=255)
                                                                                                                                            aaatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCa           <  1:2135737/148‑1 (MQ=255)
                                                                                                                                             aatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCaa          >  1:257244/1‑149 (MQ=255)
                                                                                                                                             aatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCaa          >  2:2590517/1‑149 (MQ=255)
                                                                                                                                             aatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCaa          >  1:854049/1‑149 (MQ=255)
                                                                                                                                               tcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTGTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAgg        >  2:1844520/1‑149 (MQ=255)
                                                                                                                                                caGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGc       <  2:639892/149‑1 (MQ=255)
                                                                                                                                                 aGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCg      <  2:2864948/149‑1 (MQ=255)
                                                                                                                                                  ggAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGa     >  1:1907165/1‑149 (MQ=255)
                                                                                                                                                   gAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAg    <  2:3390796/149‑1 (MQ=255)
                                                                                                                                                   gAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAg    >  1:715436/1‑149 (MQ=255)
                                                                                                                                                    aaGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGc   >  2:4953320/1‑149 (MQ=255)
                                                                                                                                                     aGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGcc  <  2:4026983/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: