Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,833,730 C→A 100% R3S (CGT→AGT)  ynjE → molybdopterin synthase sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,833,7300CA97.1% 1010.1 / 15.3 292R3S (CGT→AGT) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base A (168/115);  total (173/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 294 total aligned reads displayed.

CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833583‑1833876
                                                                                                                                                   |                                                                                                                                                  
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                  >  2:4342786/1‑150 (MQ=255)
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                  >  2:1247707/1‑150 (MQ=255)
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                  >  2:1735878/1‑150 (MQ=255)
 tcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtg                                                                                                                                                 >  1:3280511/1‑150 (MQ=255)
 tcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtg                                                                                                                                                 >  1:904298/1‑150 (MQ=255)
     ggATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtc                                                                                                                                             <  1:1334141/150‑1 (MQ=255)
        tGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCa                                                                                                                                          <  1:438854/150‑1 (MQ=255)
             ctgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGa                                                                                                                                     >  1:2925853/1‑150 (MQ=255)
             ctgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGa                                                                                                                                     >  1:326062/1‑150 (MQ=255)
               gctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGAcc                                                                                                                                   >  1:4681863/1‑150 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                  >  2:2906613/1‑150 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                  <  1:576986/150‑1 (MQ=255)
                 tgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcgc                                                                                                                                 >  2:3743433/1‑150 (MQ=255)
                          tgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAAt                                                                                                                        >  2:2356278/1‑150 (MQ=255)
                           gttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATg                                                                                                                       <  1:4458885/150‑1 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  <  2:1960241/150‑1 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  <  2:1530982/150‑1 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  >  2:3417631/1‑150 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  >  1:3883244/1‑150 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  <  2:3824740/150‑1 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  <  2:972004/150‑1 (MQ=255)
                                 gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTa                                                                                                                 >  2:2072611/1‑150 (MQ=255)
                                   tatGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAgg                                                                                                               >  1:2895491/1‑150 (MQ=255)
                                    atGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAggg                                                                                                              <  2:1550276/150‑1 (MQ=255)
                                      ggCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggct                                                                                                            >  2:3139255/1‑150 (MQ=255)
                                        cgcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctgg                                                                                                          >  2:3550356/1‑150 (MQ=255)
                                         gcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggc                                                                                                         >  1:4432833/1‑150 (MQ=255)
                                          cAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggct                                                                                                        <  1:1273148/150‑1 (MQ=255)
                                             tGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgc                                                                                                     >  1:1972731/1‑150 (MQ=255)
                                               gTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgcgc                                                                                                   >  1:4199130/1‑150 (MQ=255)
                                                tCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGct                                                                                                  <  1:307327/150‑1 (MQ=255)
                                                  aGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttc                                                                                                >  2:685349/1‑150 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              >  2:573066/1‑150 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              <  1:2174970/150‑1 (MQ=255)
                                                    aaaaTGACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              <  2:4588848/150‑1 (MQ=255)
                                                      aaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCACTGGCTTTAGGGCTGGCTTGCGCTTCTTCg                                                                                            >  1:2299204/1‑150 (MQ=255)
                                                       aTTACAACACACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGt                                                                                           <  2:1972055/150‑1 (MQ=255)
                                                        ttACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTg                                                                                          >  1:2122525/1‑150 (MQ=255)
                                                        ttACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTg                                                                                          >  1:2447803/1‑150 (MQ=255)
                                                         tACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                         <  2:4084592/150‑1 (MQ=255)
                                                         tACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                         <  1:2745754/150‑1 (MQ=255)
                                                           caacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGc                                                                                       >  1:943301/1‑150 (MQ=255)
                                                           caacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGc                                                                                       <  1:2798827/150‑1 (MQ=255)
                                                            aacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCACTGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGcc                                                                                      >  1:1483102/1‑150 (MQ=255)
                                                            aacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGcc                                                                                      <  2:2686287/150‑1 (MQ=255)
                                                             aatatcTGCGATAACGCAAAGTTTTTCTCAATGAGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATTGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCg                                                                                     <  1:4313951/145‑1 (MQ=255)
                                                              caaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGc                                                                                    <  1:4056990/150‑1 (MQ=255)
                                                               aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCt                                                                                   >  2:2731306/1‑150 (MQ=255)
                                                               aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCt                                                                                   >  1:3420738/1‑150 (MQ=255)
                                                               aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCt                                                                                   >  2:1659103/1‑150 (MQ=255)
                                                                aaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTg                                                                                  >  1:222789/1‑150 (MQ=255)
                                                                aaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTg                                                                                  >  2:2453251/1‑150 (MQ=255)
                                                                      gATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTgg                                                                            >  2:4099435/1‑150 (MQ=255)
                                                                       aTAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGc                                                                           >  2:645654/1‑150 (MQ=255)
                                                                         aaCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGa                                                                         <  2:1825519/150‑1 (MQ=255)
                                                                         aaCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGa                                                                         <  2:4209742/150‑1 (MQ=255)
                                                                          aCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGaa                                                                        <  1:1913682/150‑1 (MQ=255)
                                                                          aCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGaa                                                                        >  1:1119282/1‑150 (MQ=255)
                                                                           cGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAg                                                                       <  1:3227923/150‑1 (MQ=255)
                                                                             cAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGcc                                                                     <  1:1144796/150‑1 (MQ=255)
                                                                             cAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGcc                                                                     <  1:3537048/150‑1 (MQ=255)
                                                                                 gTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCtt                                                                 >  1:1596909/1‑150 (MQ=255)
                                                                                    tttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTaca                                                              >  2:4072551/1‑150 (MQ=255)
                                                                                      tctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACt                                                            >  1:489583/1‑150 (MQ=255)
                                                                                       ctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                           <  2:2847558/150‑1 (MQ=255)
                                                                                       ctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                           >  1:4017515/1‑150 (MQ=255)
                                                                                         cAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGa                                                         >  2:4011840/1‑150 (MQ=255)
                                                                                           aTGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGAcc                                                       >  1:1361401/1‑150 (MQ=255)
                                                                                            tGCGTCAGTTCAGAATGGCGCACTAAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAACTGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGATGAACTGGCGAAGCCTCTTACACTTGACCa                                                      <  1:2846280/150‑1 (MQ=255)
                                                                                             gCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAg                                                     <  2:2277669/150‑1 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                   >  2:2018346/1‑150 (MQ=255)
                                                                                                tCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCtt                                                  <  1:1551586/150‑1 (MQ=255)
                                                                                                  aGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTca                                                >  1:3492852/1‑150 (MQ=255)
                                                                                                     tCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaaca                                             >  1:3811633/1‑150 (MQ=255)
                                                                                                      cAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaa                                            >  1:4005636/1‑150 (MQ=255)
                                                                                                        gAATGGCGCACTCAAAACTACATTGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGTCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaaca                                          <  2:841010/150‑1 (MQ=255)
                                                                                                        gAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaaca                                          <  1:4751789/150‑1 (MQ=255)
                                                                                                         aaTGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaacaa                                         >  1:2522453/1‑150 (MQ=255)
                                                                                                         aaTGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaacaa                                         >  1:723954/1‑150 (MQ=255)
                                                                                                           tGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaaa                                       <  1:3743443/150‑1 (MQ=255)
                                                                                                               gcACTCAAAACGACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGc                                   <  2:4117083/150‑1 (MQ=255)
                                                                                                                  cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaaa                                >  2:2633981/1‑150 (MQ=255)
                                                                                                                  cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaaa                                >  2:941523/1‑150 (MQ=255)
                                                                                                                   acAAAACTACAATGTCGTGATTTTCGATTAAAAGTGTTTCTCAAATGACCTCGCTGGCAATGGCTTTATGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaaat                               <  2:3881706/149‑2 (MQ=255)
                                                                                                                    cAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGc                              <  2:2381851/150‑1 (MQ=255)
                                                                                                                       aaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgat                           <  2:1392880/150‑1 (MQ=255)
                                                                                                                        aCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgata                          >  2:1353778/1‑150 (MQ=255)
                                                                                                                         cTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatag                         >  1:4243052/1‑150 (MQ=255)
                                                                                                                         cTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatag                         >  1:2402875/1‑150 (MQ=255)
                                                                                                                             aaTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATAc                     >  1:1704878/1‑150 (MQ=255)
                                                                                                                               tGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAATCGATAGATACTc                   <  2:3266462/150‑1 (MQ=255)
                                                                                                                               tGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTc                   >  2:2716223/1‑150 (MQ=255)
                                                                                                                                gTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCg                  <  1:1983379/150‑1 (MQ=255)
                                                                                                                                     gATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCa             <  2:4205159/150‑1 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    >  1:2628904/1‑150 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    <  1:2985523/150‑1 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    >  1:4534371/1‑150 (MQ=255)
                                                                                                                                               gAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTta   <  2:3624294/150‑1 (MQ=255)
                                                                                                                                               gAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTta   <  2:3499397/150‑1 (MQ=255)
                                                                                                                                                aaaaGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTtat  >  1:3679516/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833583‑1833876

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: