Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,206,219 A→G 5.4% I140V (ATA→GTA)  mcrA → putative 5‑methylcytosine/5‑hydroxymethylcytosine‑specific restriction nuclease; 5‑methylcytosine DNA binding protein; e14 prophage gene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,206,2190AG 5.4% 772.4 / 24.6 223I140V (ATA→GTA) mcrAputative 5‑methylcytosine/5‑hydroxymethylcytosine‑specific restriction nuclease; 5‑methylcytosine DNA binding protein; e14 prophage gene
Reads supporting (aligned to +/- strand):  ref base A (103/108);  new base G (7/5);  total (110/113)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 245 total aligned reads displayed.

TAATAGCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTA  >  CP009273/1206071‑1206367
                                                                                                                                                    |                                                                                                                                                    
  aTAGCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCtata                                                                                                                                                    >  1:3509808/1‑149 (MQ=255)
   tAGCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataa                                                                                                                                                   <  1:2897978/149‑1 (MQ=255)
   tAGCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataa                                                                                                                                                   <  1:2336331/149‑1 (MQ=255)
    aGCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataat                                                                                                                                                  >  2:3822333/1‑149 (MQ=255)
     gCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataata                                                                                                                                                 <  1:3208544/149‑1 (MQ=255)
     gCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataata                                                                                                                                                 >  2:37237/1‑149 (MQ=255)
     gCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataata                                                                                                                                                 >  1:2567335/1‑149 (MQ=255)
     gCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTataata                                                                                                                                                 >  1:2048674/1‑149 (MQ=255)
        cccGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGt                                                                                                                                              >  1:2295349/1‑149 (MQ=255)
           gcgATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGga                                                                                                                                           >  1:123602/1‑149 (MQ=255)
           gcgATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGga                                                                                                                                           <  2:2211339/149‑1 (MQ=255)
              atatACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGaa                                                                                                                                        <  1:4307326/149‑1 (MQ=255)
                 taCGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTa                                                                                                                                     >  2:3761904/1‑149 (MQ=255)
                    gCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTTGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCa                                                                                                                                  >  2:4431798/1‑149 (MQ=255)
                    gCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCa                                                                                                                                  <  2:232794/149‑1 (MQ=255)
                     cTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCag                                                                                                                                 >  1:4217147/1‑149 (MQ=255)
                       tGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCagag                                                                                                                               >  2:3091385/1‑149 (MQ=255)
                          aGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCtt                                                                                                                            >  1:2270644/1‑149 (MQ=255)
                          aGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCtt                                                                                                                            >  2:4614268/1‑149 (MQ=255)
                           gATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCttt                                                                                                                           <  1:152046/149‑1 (MQ=255)
                               tgtgGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCTGAGCTTTcaca                                                                                                                       <  1:2285480/149‑1 (MQ=255)
                               tgtgGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTcaca                                                                                                                       >  2:2979848/1‑149 (MQ=255)
                                 tgGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGc                                                                                                                     >  1:2386294/1‑149 (MQ=255)
                                 tgGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGc                                                                                                                     >  1:3952885/1‑149 (MQ=255)
                                    ttATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTa                                                                                                                  >  2:4031136/1‑149 (MQ=255)
                                     tATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTAc                                                                                                                 >  1:5002210/1‑149 (MQ=255)
                                      aTCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACa                                                                                                                >  1:2612835/1‑149 (MQ=255)
                                        cAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAga                                                                                                              >  2:2768818/1‑149 (MQ=255)
                                           gCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAgatga                                                                                                           <  1:4078150/149‑1 (MQ=255)
                                           gCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAgatga                                                                                                           <  2:1482189/149‑1 (MQ=255)
                                           gCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAgatga                                                                                                           <  2:1127478/149‑1 (MQ=255)
                                                ttttAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                                      >  2:740497/1‑149 (MQ=255)
                                                ttttAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                                      >  1:4842990/1‑149 (MQ=255)
                                                ttttAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAAt                                                                                                      <  1:4910789/149‑1 (MQ=255)
                                                 tttAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTGTAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATc                                                                                                     >  1:1373285/1‑149 (MQ=255)
                                                        gTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTg                                                                                              >  2:2827368/1‑149 (MQ=255)
                                                            ggggAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAata                                                                                          <  2:3424351/149‑1 (MQ=255)
                                                             gggAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAatat                                                                                         >  1:1739429/1‑149 (MQ=255)
                                                                  aaaGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAggg                                                                                    >  1:3043812/1‑149 (MQ=255)
                                                                  aaaGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAggg                                                                                    <  2:3231145/149‑1 (MQ=255)
                                                                  aaaGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAggg                                                                                    >  2:4220032/1‑149 (MQ=255)
                                                                  aaaGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAggg                                                                                    <  1:4400713/149‑1 (MQ=255)
                                                                    aGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGtt                                                                                  <  2:1727331/149‑1 (MQ=255)
                                                                     ggAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTa                                                                                 <  1:4632795/149‑1 (MQ=255)
                                                                      gAAATTCCTTCCGTCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATATTGACAGTTTTTGGGCTTCTATATTACGTGGAGAACTATCATAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGTTTAg                                                                                <  2:2491034/149‑1 (MQ=255)
                                                                      gAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAg                                                                                <  2:5070112/149‑1 (MQ=255)
                                                                          ttccttccGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTaaa                                                                            <  1:1328740/149‑1 (MQ=255)
                                                                           tccttccGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAaat                                                                           <  1:5139263/149‑1 (MQ=255)
                                                                             cttccGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAaatta                                                                         >  1:370623/1‑149 (MQ=255)
                                                                              ttccGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAaattaa                                                                        <  2:1771427/149‑1 (MQ=255)
                                                                               tccGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAaattaat                                                                       <  2:4323882/149‑1 (MQ=255)
                                                                                  ggCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAAttaattaa                                                                    <  1:2609228/149‑1 (MQ=255)
                                                                                       tAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGaaaa                                                               >  2:1043822/1‑149 (MQ=255)
                                                                                          aacgaaacgaATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGt                                                            <  2:2032350/149‑1 (MQ=255)
                                                                                             gaaacgaaTATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTGTAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGa                                                         <  1:3842034/149‑1 (MQ=255)
                                                                                              aaacgaaTATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAg                                                        >  1:1120290/1‑149 (MQ=255)
                                                                                               aacgaaTATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGt                                                       >  1:2864319/1‑149 (MQ=255)
                                                                                                  gaaTATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCaa                                                    <  2:2528063/149‑1 (MQ=255)
                                                                                                    atatTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAAcc                                                  >  1:4713821/1‑149 (MQ=255)
                                                                                                       ttGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGa                                               >  2:5173792/1‑149 (MQ=255)
                                                                                                       ttGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGa                                               <  1:3975933/149‑1 (MQ=255)
                                                                                                       ttGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGa                                               >  2:3557154/1‑149 (MQ=255)
                                                                                                        tGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAg                                              <  1:2223119/149‑1 (MQ=255)
                                                                                                         gATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAgg                                             >  1:1529154/1‑149 (MQ=255)
                                                                                                         gATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAgg                                             >  2:40137/1‑149 (MQ=255)
                                                                                                           tAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGc                                           <  1:1816142/149‑1 (MQ=255)
                                                                                                            aaaTGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCt                                          >  2:5136278/1‑149 (MQ=255)
                                                                                                              aTGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTcc                                        <  1:2960866/149‑1 (MQ=255)
                                                                                                              aTGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTcc                                        <  2:2722972/149‑1 (MQ=255)
                                                                                                                    cAGGTATTTATAGTGACAGTTTTTGGGCGTCTGTAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGaaa                                  >  1:3452176/1‑149 (MQ=255)
                                                                                                                    cAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGaaa                                  >  1:4881411/1‑149 (MQ=255)
                                                                                                                      ggTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAAcc                                >  2:4209626/1‑149 (MQ=255)
                                                                                                                      ggTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAAcc                                >  1:1027556/1‑149 (MQ=255)
                                                                                                                        tATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAg                              >  2:4528381/1‑149 (MQ=255)
                                                                                                                         aTTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGt                             <  2:4912666/149‑1 (MQ=255)
                                                                                                                          tttATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGtt                            <  1:2827368/149‑1 (MQ=255)
                                                                                                                          tttATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGtt                            >  2:3220154/1‑149 (MQ=255)
                                                                                                                          tttATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGtt                            <  2:3308122/149‑1 (MQ=255)
                                                                                                                          tttATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGtt                            >  2:621220/1‑149 (MQ=255)
                                                                                                                          tttATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGtt                            >  1:4460511/1‑149 (MQ=255)
                                                                                                                           ttATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTg                           <  1:3014620/149‑1 (MQ=255)
                                                                                                                            tataGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGa                          <  1:241382/149‑1 (MQ=255)
                                                                                                                             ataGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAg                         <  1:3749221/149‑1 (MQ=255)
                                                                                                                              taGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAgg                        <  1:16678/149‑1 (MQ=255)
                                                                                                                               aGTGACAGTTTTTGGGCGTCTGTAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGt                       <  2:1373389/149‑1 (MQ=255)
                                                                                                                                gTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTa                      >  2:2367637/1‑149 (MQ=255)
                                                                                                                                gTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTa                      >  1:1903729/1‑149 (MQ=255)
                                                                                                                                      gTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGa                >  2:191678/1‑149 (MQ=255)
                                                                                                                                      gTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGa                <  1:4437338/149‑1 (MQ=255)
                                                                                                                                      gTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGa                <  1:1841964/149‑1 (MQ=255)
                                                                                                                                       tttttGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGAc               >  1:1947942/1‑149 (MQ=255)
                                                                                                                                          ttGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTAc            <  2:3509808/149‑1 (MQ=255)
                                                                                                                                               cGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACaaaaa       <  1:1964988/149‑1 (MQ=255)
                                                                                                                                                 tCTGTAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGt     >  1:2171611/1‑149 (MQ=255)
                                                                                                                                                 tCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGt     <  2:2048674/149‑1 (MQ=255)
                                                                                                                                                 tCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGt     <  1:311851/149‑1 (MQ=255)
                                                                                                                                                   tataATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGttt   <  1:2392836/149‑1 (MQ=255)
                                                                                                                                                    ataattCGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATAGGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTa  <  2:2616162/149‑1 (MQ=255)
                                                                                                                                                    ataataCGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTa  <  2:4451121/149‑1 (MQ=255)
                                                                                                                                                    ataataCGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTa  >  1:2903845/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAATAGCCCCCGCGATATACGCTTGAAGATGTGTGGTTATCAGGCTTATTTTAGTCGTACGGGGAGAAAGGAAATTCCTTCCGGCAATAGAACGAAACGAATATTGATAAATGTTCCAGGTATTTATAGTGACAGTTTTTGGGCGTCTATAATACGTGGAGAACTATCAGAGCTTTCACAGCCTACAGATGATGAATCGCTTCTGAATATGAGGGTTAGTAAATTAATTAAGAAAACGTTGAGTCAACCCGAGGGCTCCAGGAAACCAGTTGAGGTAGAAAGACTACAAAAAGTTTA  >  CP009273/1206071‑1206367

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: