Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,567,552 C→A 10.6% 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,5520CA10.6% 929.5 / 55.0 255R75L (CGA→CTA) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base C (116/112);  new base A (12/15);  total (128/127)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 257 total aligned reads displayed.

TTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTT  >  CP009273/4567403‑4567700
                                                                                                                                                     |                                                                                                                                                    
 tGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCg                                                                                                                                                     >  1:405476/1‑150 (MQ=255)
 tGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCg                                                                                                                                                     >  1:595465/1‑150 (MQ=255)
     tttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCt                                                                                                                                                 <  1:1520948/150‑1 (MQ=255)
       tAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCttt                                                                                                                                               <  2:487601/150‑1 (MQ=255)
        aaCTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTg                                                                                                                                              >  1:2392889/1‑150 (MQ=255)
          cTTTTCATGCTTAATTAATATGGCTAATTTGATTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAATTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTTAt                                                                                                                                            >  1:1299144/1‑150 (MQ=255)
          cTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                                                            >  1:1178912/1‑150 (MQ=255)
          cTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                                                            <  1:3121676/150‑1 (MQ=255)
             ttCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGcc                                                                                                                                         >  2:1923937/1‑150 (MQ=255)
             ttCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGcc                                                                                                                                         <  1:944612/150‑1 (MQ=255)
              tCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCa                                                                                                                                        <  1:1366606/150‑1 (MQ=255)
              tCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCa                                                                                                                                        <  2:1712623/150‑1 (MQ=255)
              tCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCa                                                                                                                                        <  2:2126807/150‑1 (MQ=255)
               cATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAg                                                                                                                                       <  2:4091599/150‑1 (MQ=255)
                aTGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAgg                                                                                                                                      <  2:2231341/150‑1 (MQ=255)
                 tGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATTATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGa                                                                                                                                     >  2:2571700/1‑150 (MQ=255)
                   cTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAAt                                                                                                                                   <  1:2732967/150‑1 (MQ=255)
                     taattaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGa                                                                                                                                 <  1:4556846/150‑1 (MQ=255)
                      aattaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAt                                                                                                                                >  2:2081784/1‑150 (MQ=255)
                              tGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGgtgt                                                                                                                        >  2:1150367/1‑150 (MQ=255)
                                 cTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAg                                                                                                                     >  2:1024779/1‑150 (MQ=255)
                                 cTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAg                                                                                                                     >  1:3835579/1‑150 (MQ=255)
                                    aTGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGag                                                                                                                  >  2:3194618/1‑149 (MQ=255)
                                     tgtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGAtt                                                                                                                 >  2:2428604/1‑150 (MQ=255)
                                     tgtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGAtt                                                                                                                 <  1:3632267/150‑1 (MQ=255)
                                      gtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTg                                                                                                                >  2:3814663/1‑150 (MQ=255)
                                      gtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTg                                                                                                                >  1:3735334/1‑150 (MQ=255)
                                       tgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGt                                                                                                               >  2:2037279/1‑150 (MQ=255)
                                       tgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGt                                                                                                               >  1:3602469/1‑150 (MQ=255)
                                         cTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTaa                                                                                                             >  2:1858523/1‑150 (MQ=255)
                                          tttttATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAt                                                                                                            <  1:3240570/150‑1 (MQ=255)
                                          tttttATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAt                                                                                                            >  2:2261246/1‑150 (MQ=255)
                                              tataATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaag                                                                                                        >  1:648447/1‑150 (MQ=255)
                                                taATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagct                                                                                                      <  1:4125076/150‑1 (MQ=255)
                                                 aaTTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctc                                                                                                     >  1:3310376/1‑150 (MQ=255)
                                                   ttCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctcaa                                                                                                   <  2:359361/150‑1 (MQ=255)
                                                      ggTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGcc                                                                                                >  2:1418277/1‑150 (MQ=255)
                                                        ttAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCt                                                                                              <  1:2879159/150‑1 (MQ=255)
                                                         tAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTc                                                                                             <  1:2470564/150‑1 (MQ=255)
                                                             cAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCg                                                                                         >  1:3359774/1‑150 (MQ=255)
                                                             cAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGATCAAGCCCTCGGCg                                                                                         >  2:2169533/1‑150 (MQ=255)
                                                              aGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGt                                                                                        >  1:1396404/1‑150 (MQ=255)
                                                               gCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTg                                                                                       <  2:3953932/150‑1 (MQ=255)
                                                                   gtCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAg                                                                                   <  1:2316736/150‑1 (MQ=255)
                                                                     cTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCt                                                                                 >  2:4381315/1‑150 (MQ=255)
                                                                     cTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCt                                                                                 <  2:474142/150‑1 (MQ=255)
                                                                     cTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCt                                                                                 <  1:1897261/150‑1 (MQ=255)
                                                                     cTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCt                                                                                 <  2:2999998/150‑1 (MQ=255)
                                                                       tCATTAATCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGt                                                                               <  1:1966635/150‑1 (MQ=255)
                                                                         aTTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaa                                                                             >  1:782282/1‑150 (MQ=255)
                                                                          ttAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaaa                                                                            >  1:377547/1‑150 (MQ=255)
                                                                           tAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaaaa                                                                           >  2:412814/1‑150 (MQ=255)
                                                                           tAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaaaa                                                                           <  2:162798/150‑1 (MQ=255)
                                                                           tAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGAAATCAATCTCAAGCCATCGTCGTGAAAGATGTaaaa                                                                           >  2:1805923/1‑150 (MQ=255)
                                                                                atatCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCt                                                                      >  1:1411439/1‑150 (MQ=255)
                                                                                   tCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCttta                                                                   >  1:321452/1‑150 (MQ=255)
                                                                                    cAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttat                                                                  <  1:1007270/150‑1 (MQ=255)
                                                                                     aaaaGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttatt                                                                 >  1:2784682/1‑150 (MQ=255)
                                                                                      aaaGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttattt                                                                <  2:2742853/150‑1 (MQ=255)
                                                                                      aaaGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttattt                                                                >  1:590481/1‑150 (MQ=255)
                                                                                          tACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAata                                                            <  2:2455194/150‑1 (MQ=255)
                                                                                           aCTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatac                                                           >  1:1797753/1‑150 (MQ=255)
                                                                                            cTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAataca                                                          >  1:1943749/1‑150 (MQ=255)
                                                                                             tGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacat                                                         >  2:891537/1‑150 (MQ=255)
                                                                                             tGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGACCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacat                                                         >  2:2366291/1‑150 (MQ=255)
                                                                                              gACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacata                                                        >  2:1350442/1‑150 (MQ=255)
                                                                                               aCGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAatacatac                                                       <  2:3370413/150‑1 (MQ=255)
                                                                                                cGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATAcc                                                      <  1:3449019/150‑1 (MQ=255)
                                                                                                    tttCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAg                                                  <  1:2572057/150‑1 (MQ=255)
                                                                                                    tttCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCAGTAAAACCCCTTTATTTAATACATACCCAAg                                                  >  1:2070757/1‑150 (MQ=255)
                                                                                                     ttCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAgg                                                 >  1:3970599/1‑150 (MQ=255)
                                                                                                     ttCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAgg                                                 <  2:2583828/150‑1 (MQ=255)
                                                                                                        ccATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGata                                              <  1:795568/150‑1 (MQ=255)
                                                                                                          aTGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATc                                            <  1:4345296/150‑1 (MQ=255)
                                                                                                               tAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaagaa                                       >  2:1742528/1‑150 (MQ=255)
                                                                                                               tAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaagaa                                       <  2:1396268/150‑1 (MQ=255)
                                                                                                                aaGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaagaag                                      >  2:139266/1‑150 (MQ=255)
                                                                                                                aaGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaagaag                                      >  2:3032017/1‑150 (MQ=255)
                                                                                                                aaGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCaagaag                                      >  2:4169970/1‑150 (MQ=255)
                                                                                                                    tGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTg                                  <  2:353834/150‑1 (MQ=255)
                                                                                                                     ggAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGt                                 >  2:1939703/1‑150 (MQ=255)
                                                                                                                     ggAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGt                                 >  2:4305441/1‑150 (MQ=255)
                                                                                                                      gAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGGGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGtt                                >  1:498002/1‑150 (MQ=255)
                                                                                                                      gAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGtt                                <  2:50732/150‑1 (MQ=255)
                                                                                                                       aaGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTa                               <  2:564630/150‑1 (MQ=255)
                                                                                                                       aaGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTa                               <  1:4564853/150‑1 (MQ=255)
                                                                                                                        aGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTAc                              >  1:987810/1‑150 (MQ=255)
                                                                                                                        aGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATAGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTAc                              >  2:2807545/1‑150 (MQ=255)
                                                                                                                               tATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTAt                       <  2:648607/150‑1 (MQ=255)
                                                                                                                                aTTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATg                      <  1:1024779/150‑1 (MQ=255)
                                                                                                                                aTTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATg                      >  2:2917326/1‑150 (MQ=255)
                                                                                                                                 ttGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATgg                     <  2:3931606/150‑1 (MQ=255)
                                                                                                                                  tGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGc                    <  2:2279457/150‑1 (MQ=255)
                                                                                                                                          cAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGt            <  1:2066460/150‑1 (MQ=255)
                                                                                                                                           aaaCTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGtt           <  1:441076/150‑1 (MQ=255)
                                                                                                                                            aaCTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGttt          <  1:877809/150‑1 (MQ=255)
                                                                                                                                            aaCTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGttt          <  1:1042747/150‑1 (MQ=255)
                                                                                                                                            aaCTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGttt          <  2:675314/150‑1 (MQ=255)
                                                                                                                                                 tATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTg     >  1:960889/1‑150 (MQ=255)
                                                                                                                                                  aTTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGc    >  2:114627/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTT  >  CP009273/4567403‑4567700

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: