Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,567,552 C→A 14.7% R75L (CGA→CTA)  mcrC ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,567,5520CA14.7% 873.1 / 89.4 266R75L (CGA→CTA) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base C (111/116);  new base A (21/18);  total (132/134)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 268 total aligned reads displayed.

TGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTA  >  CP009273/4567404‑4567701
                                                                                                                                                    |                                                                                                                                                     
tGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCg                                                                                                                                                      >  2:2061536/1‑150 (MQ=255)
tGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCg                                                                                                                                                      <  1:4536797/150‑1 (MQ=255)
    tttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCt                                                                                                                                                  <  1:1093586/150‑1 (MQ=255)
     ttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCtt                                                                                                                                                 <  2:3407317/150‑1 (MQ=255)
     ttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCtt                                                                                                                                                 <  1:854259/150‑1 (MQ=255)
      tAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCttt                                                                                                                                                <  2:1883311/150‑1 (MQ=255)
        aCTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGa                                                                                                                                              >  1:4545918/1‑150 (MQ=255)
        aCTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGa                                                                                                                                              >  1:4198911/1‑150 (MQ=255)
         cTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                                                             >  1:1447597/1‑150 (MQ=255)
         cTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                                                             <  2:3687219/150‑1 (MQ=255)
         cTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                                                             <  2:1754895/150‑1 (MQ=255)
           tttCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGc                                                                                                                                           >  1:4075794/1‑150 (MQ=255)
             tCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCa                                                                                                                                         >  1:103055/1‑150 (MQ=255)
              cATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAg                                                                                                                                        >  1:336047/1‑150 (MQ=255)
              cATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAg                                                                                                                                        >  1:336104/1‑150 (MQ=255)
                tGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAATCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGa                                                                                                                                      >  2:89449/1‑150 (MQ=255)
                tGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGa                                                                                                                                      >  1:2145762/1‑150 (MQ=255)
                 gCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGaa                                                                                                                                     >  2:529131/1‑150 (MQ=255)
                  cTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAAt                                                                                                                                    >  2:1375648/1‑150 (MQ=255)
                   ttaattaaTATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATg                                                                                                                                   <  1:291848/150‑1 (MQ=255)
                     aattaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAt                                                                                                                                 <  1:3761407/150‑1 (MQ=255)
                      attaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAtc                                                                                                                                >  1:2713996/1‑150 (MQ=255)
                      attaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAtc                                                                                                                                >  1:2372944/1‑150 (MQ=255)
                       ttaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAtct                                                                                                                               <  1:2773619/150‑1 (MQ=255)
                       ttaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAtct                                                                                                                               <  1:3174620/150‑1 (MQ=255)
                        taatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAtctc                                                                                                                              <  1:2533897/150‑1 (MQ=255)
                        taatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGAtctc                                                                                                                              <  1:3632888/150‑1 (MQ=255)
                         aatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCg                                                                                                                             >  2:123602/1‑150 (MQ=255)
                           tatGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGgt                                                                                                                           >  1:2169940/1‑150 (MQ=255)
                              ggCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGtt                                                                                                                        >  1:1669432/1‑150 (MQ=255)
                               gCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTa                                                                                                                       >  1:2842507/1‑150 (MQ=255)
                               gCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTa                                                                                                                       <  1:4316759/150‑1 (MQ=255)
                               gCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTa                                                                                                                       >  1:1488539/1‑150 (MQ=255)
                                cTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAg                                                                                                                      <  2:4137560/150‑1 (MQ=255)
                                cTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAg                                                                                                                      >  2:1049949/1‑150 (MQ=255)
                                    tgtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGAtt                                                                                                                  <  1:4327983/150‑1 (MQ=255)
                                       gCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTa                                                                                                               <  2:3703795/150‑1 (MQ=255)
                                        cTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCAAGATAAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGACCTTTGAAGCCAGGAATGAACTCGGTATTAGCATTGTaa                                                                                                              >  1:1941238/1‑150 (MQ=255)
                                              ataATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagc                                                                                                        >  2:2382833/1‑150 (MQ=255)
                                              ataATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagc                                                                                                        <  1:1920402/150‑1 (MQ=255)
                                               taATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagct                                                                                                       >  1:3113266/1‑150 (MQ=255)
                                                 aTTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctca                                                                                                     >  1:4471802/1‑150 (MQ=255)
                                                 aTTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctca                                                                                                     <  2:3459014/150‑1 (MQ=255)
                                                  ttCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctcaa                                                                                                    >  2:2328100/1‑150 (MQ=255)
                                                  ttCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctcaa                                                                                                    <  2:3515090/150‑1 (MQ=255)
                                                   tCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctcaag                                                                                                   >  1:1957572/1‑150 (MQ=255)
                                                     ggTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGcc                                                                                                 >  2:2650132/1‑150 (MQ=255)
                                                      gTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGccc                                                                                                >  1:2289411/1‑150 (MQ=255)
                                                        tAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTc                                                                                              <  1:1395643/150‑1 (MQ=255)
                                                          gCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCgg                                                                                            <  2:4318096/150‑1 (MQ=255)
                                                            cAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCg                                                                                          >  1:3913375/1‑150 (MQ=255)
                                                            cAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCg                                                                                          >  1:3050282/1‑150 (MQ=255)
                                                              gCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTg                                                                                        <  1:2810451/150‑1 (MQ=255)
                                                               cGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGa                                                                                       <  1:3861576/150‑1 (MQ=255)
                                                                gtgtCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGaa                                                                                      <  1:191812/150‑1 (MQ=255)
                                                                     ttCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTg                                                                                 <  2:1309518/150‑1 (MQ=255)
                                                                       cATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTa                                                                               <  2:2654706/150‑1 (MQ=255)
                                                                       cATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTa                                                                               <  2:1581545/150‑1 (MQ=255)
                                                                        aTTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaa                                                                              <  1:533046/150‑1 (MQ=255)
                                                                        aTTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaa                                                                              >  2:4070064/1‑150 (MQ=255)
                                                                          tAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaaaa                                                                            >  1:4332980/1‑150 (MQ=255)
                                                                           aaGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAAc                                                                           <  2:4170612/150‑1 (MQ=255)
                                                                                tatCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtt                                                                      <  2:1864920/150‑1 (MQ=255)
                                                                                   cAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttat                                                                   >  2:1477298/1‑150 (MQ=255)
                                                                                   cAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttat                                                                   >  2:458092/1‑150 (MQ=255)
                                                                                   cAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttat                                                                   >  2:2977334/1‑150 (MQ=255)
                                                                                   cAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttat                                                                   >  2:664306/1‑150 (MQ=255)
                                                                                    aaaaGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttatt                                                                  <  2:2544888/150‑1 (MQ=255)
                                                                                      aaGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttattta                                                                >  2:1854131/1‑150 (MQ=255)
                                                                                            tGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacat                                                          <  1:2425269/150‑1 (MQ=255)
                                                                                              aCGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacatac                                                        <  1:2223517/150‑1 (MQ=255)
                                                                                              aCGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacatac                                                        >  1:3574320/1‑150 (MQ=255)
                                                                                                ggTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATAccc                                                      >  1:3706451/1‑150 (MQ=255)
                                                                                                 gTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATAcccc                                                     <  1:3966277/150‑1 (MQ=255)
                                                                                                   tttCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAg                                                   >  2:4295454/1‑150 (MQ=255)
                                                                                                   tttCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAg                                                   >  2:2771262/1‑150 (MQ=255)
                                                                                                      cccATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGat                                                >  1:1518684/1‑150 (MQ=255)
                                                                                                          tGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCa                                            <  2:1766323/150‑1 (MQ=255)
                                                                                                          tGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCa                                            >  1:4439679/1‑150 (MQ=255)
                                                                                                           gATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaa                                           <  2:1516405/150‑1 (MQ=255)
                                                                                                           gATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaa                                           >  1:607810/1‑150 (MQ=255)
                                                                                                            aTTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaag                                          >  1:569336/1‑150 (MQ=255)
                                                                                                               aaGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCGTTATTTAATACATACCCCAGGATATCaagaag                                       >  2:2883864/1‑150 (MQ=255)
                                                                                                                  aTGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGAtt                                    <  2:2145762/150‑1 (MQ=255)
                                                                                                                    ggAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGt                                  >  2:2595896/1‑150 (MQ=255)
                                                                                                                     gAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGtt                                 <  2:4529150/150‑1 (MQ=255)
                                                                                                                      aaGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTa                                >  2:3236457/1‑150 (MQ=255)
                                                                                                                       aGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTAc                               >  1:815229/1‑150 (MQ=255)
                                                                                                                        gCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTAcc                              >  1:2919835/1‑150 (MQ=255)
                                                                                                                            cgTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGt                          <  1:3198306/150‑1 (MQ=255)
                                                                                                                              tATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTAt                        >  1:2507173/1‑150 (MQ=255)
                                                                                                                                ttGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATgg                      >  2:77473/1‑150 (MQ=255)
                                                                                                                                 tGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGc                     >  1:3246997/1‑150 (MQ=255)
                                                                                                                                   ttttAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCtt                   <  2:2249510/150‑1 (MQ=255)
                                                                                                                                   ttttAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCtt                   <  2:1178648/150‑1 (MQ=255)
                                                                                                                                      tAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCaa                <  1:3693048/150‑1 (MQ=255)
                                                                                                                                               ctATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTAcc                              >  1:4114048/1‑127 (MQ=255)
                                                                                                                                               atATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCt       <  2:3288545/149‑1 (MQ=255)
                                                                                                                                                tATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTg      <  1:1708348/150‑1 (MQ=255)
                                                                                                                                                  ttCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCt    >  1:1325577/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTA  >  CP009273/4567404‑4567701

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: