Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 219,803 A→T 13.9% T162S (ACG→TCG)  gmhB → D,D‑heptose 1,7‑bisphosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273219,8030AT13.9% 833.1 / 112.9 288T162S (ACG→TCG) gmhBD,D‑heptose 1,7‑bisphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (131/117);  new base T (21/19);  total (152/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 288 total aligned reads displayed.

CATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCC  >  CP009273/219656‑219950
                                                                                                                                                   |                                                                                                                                                   
caTCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCg                                                                                                                                                   >  1:4629117/1‑150 (MQ=255)
 aTCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGc                                                                                                                                                  >  2:4142893/1‑150 (MQ=255)
   ccGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCt                                                                                                                                                <  2:2837912/150‑1 (MQ=255)
    cGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTg                                                                                                                                               <  2:2633422/150‑1 (MQ=255)
     ggggATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATTCAGGCAGAGGTTGCGGCGAACGTGGGAACAAAAGTGATGGTGCGTACGGGTAAACCTATTACTCCTGa                                                                                                                                              >  1:2006523/1‑150 (MQ=255)
     ggggATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGa                                                                                                                                              >  2:2576188/1‑150 (MQ=255)
        gATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGTTATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGAGTACGGGTAAACCTATTACGCCTGAAGc                                                                                                                                           <  1:2083250/150‑1 (MQ=255)
          tGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAg                                                                                                                                         >  1:571270/1‑150 (MQ=255)
              tttGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGaaaa                                                                                                                                     >  1:3794746/1‑150 (MQ=255)
               ttGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAAc                                                                                                                                    <  1:2464095/150‑1 (MQ=255)
               ttGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAAc                                                                                                                                    <  1:1904296/150‑1 (MQ=255)
                     gCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACgcggcg                                                                                                                              <  2:1063870/150‑1 (MQ=255)
                      cACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGcggcgg                                                                                                                             >  2:681182/1‑150 (MQ=255)
                      cACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGcggcgg                                                                                                                             >  2:2994195/1‑150 (MQ=255)
                       aCGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGa                                                                                                                            >  2:3607045/1‑150 (MQ=255)
                       aCGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGa                                                                                                                            >  1:2147956/1‑150 (MQ=255)
                       aCGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGa                                                                                                                            >  1:3259780/1‑150 (MQ=255)
                          cgATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTg                                                                                                                         >  2:4405666/1‑150 (MQ=255)
                           gATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTgg                                                                                                                        >  1:2328761/1‑150 (MQ=255)
                           gATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTgg                                                                                                                        >  1:2488956/1‑150 (MQ=255)
                              tattTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGgtg                                                                                                                     <  1:105966/150‑1 (MQ=255)
                                tttGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGtt                                                                                                                   >  2:1619997/1‑150 (MQ=255)
                                   gCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATTCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa                                                                                                                >  1:4296206/1‑150 (MQ=255)
                                   gCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa                                                                                                                >  1:1730524/1‑150 (MQ=255)
                                    cATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAAt                                                                                                               >  2:596339/1‑150 (MQ=255)
                                         tGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCt                                                                                                          <  2:2665903/150‑1 (MQ=255)
                                          gATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTg                                                                                                         >  2:4444550/1‑150 (MQ=255)
                                             atGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa                                                                                                      >  2:3046313/1‑150 (MQ=255)
                                             atGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCa                                                                                                      <  1:1395499/150‑1 (MQ=255)
                                                 ccGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGAcc                                                                                                  >  2:353860/1‑150 (MQ=255)
                                                  cGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCt                                                                                                 >  1:3282758/1‑150 (MQ=255)
                                                  cGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACATATTAAGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGACTGGCAGACCt                                                                                                 >  1:1113426/1‑150 (MQ=255)
                                                    cttcttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGc                                                                                               >  1:773359/1‑150 (MQ=255)
                                                    cttcttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGc                                                                                               >  1:3095892/1‑150 (MQ=255)
                                                     ttcttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGcc                                                                                              <  2:4305552/150‑1 (MQ=255)
                                                       cttATATGGTGGGCGATAAATTCGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGc                                                                                            >  2:3620996/1‑150 (MQ=255)
                                                          atatGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAg                                                                                         >  1:3760111/1‑150 (MQ=255)
                                                          atatGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAg                                                                                         >  2:1525853/1‑150 (MQ=255)
                                                           tatGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGc                                                                                        >  2:2803572/1‑150 (MQ=255)
                                                             tggtggGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGa                                                                                      <  2:3128337/150‑1 (MQ=255)
                                                              ggtggGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCAGCAAGCGAt                                                                                     >  2:4110286/1‑150 (MQ=255)
                                                                 gggCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaa                                                                                  <  1:662599/150‑1 (MQ=255)
                                                                 gggCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaa                                                                                  <  1:1821496/150‑1 (MQ=255)
                                                                    cGATAAATTAGAGGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATaaaaaa                                                                               >  1:990752/1‑150 (MQ=255)
                                                                     gATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAag                                                                              >  1:4165754/1‑150 (MQ=255)
                                                                        aaaTTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAagcag                                                                           <  2:3395490/150‑1 (MQ=255)
                                                                          aTTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAgcagca                                                                         >  2:3267317/1‑150 (MQ=255)
                                                                            tAGAAGATATGCAGGCAGCGGTTGCGGCGTAGGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATTGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaaa                                                                       <  2:825214/150‑1 (MQ=255)
                                                                            tAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaaa                                                                       <  2:2131001/150‑1 (MQ=255)
                                                                             agaagaTATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGTGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaaaa                                                                      <  2:2370417/150‑1 (MQ=255)
                                                                             agaagaTATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCaaaa                                                                      >  2:603647/1‑150 (MQ=255)
                                                                                agaTATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAATACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAAcc                                                                   <  1:2790780/150‑1 (MQ=255)
                                                                                agaTATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAAcc                                                                   <  1:3649932/150‑1 (MQ=255)
                                                                                 gaTATGCAGGCAGCGGTTGCGGCGAACGTGGGATCAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCATGGCATACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCg                                                                  <  2:3663136/150‑1 (MQ=255)
                                                                                 gaTATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCg                                                                  <  2:3663120/150‑1 (MQ=255)
                                                                                  atatGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCgg                                                                 >  1:3736275/1‑150 (MQ=255)
                                                                                   tatGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGc                                                                <  2:2231128/150‑1 (MQ=255)
                                                                                       caggcagCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACaa                                                            <  1:1503583/150‑1 (MQ=255)
                                                                                         ggcagCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATg                                                          <  1:3958821/150‑1 (MQ=255)
                                                                                          gcagCGGTTGCGGCGAACGTGGGAACAAAAGAGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGa                                                         >  1:632097/1‑150 (MQ=255)
                                                                                            agCGGTTGCGGCGAACGTGGGAACAATTGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTATATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGAtt                                                       <  1:3744501/150‑1 (MQ=255)
                                                                                            agCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGAtt                                                       <  1:865986/150‑1 (MQ=255)
                                                                                            agCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGATCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGAtt                                                       <  2:3812059/150‑1 (MQ=255)
                                                                                            agCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGAtt                                                       >  2:2121493/1‑150 (MQ=255)
                                                                                               ggTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTaaa                                                    >  2:507900/1‑150 (MQ=255)
                                                                                               ggTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTaaa                                                    >  2:1807089/1‑150 (MQ=255)
                                                                                               ggTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTaaa                                                    <  1:2376782/150‑1 (MQ=255)
                                                                                                 ttGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAg                                                  >  1:2320652/1‑150 (MQ=255)
                                                                                                  tGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGa                                                 <  2:3694189/150‑1 (MQ=255)
                                                                                                   gcggcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGAt                                                <  1:1876508/150‑1 (MQ=255)
                                                                                                    cggcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATg                                               <  1:819330/150‑1 (MQ=255)
                                                                                                     ggcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGa                                              >  1:3351155/1‑150 (MQ=255)
                                                                                                     ggcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCAAGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGTCAGACCTGCCGCAAGCGGTAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGa                                              >  1:13643/1‑150 (MQ=255)
                                                                                                      gcgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAg                                             <  1:2928453/150‑1 (MQ=255)
                                                                                                       cgAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGc                                            >  1:2715530/1‑150 (MQ=255)
                                                                                                        gAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCg                                           >  1:931184/1‑150 (MQ=255)
                                                                                                         aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCgg                                          >  2:2000607/1‑150 (MQ=255)
                                                                                                         aaCGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCgg                                          <  2:3794746/150‑1 (MQ=255)
                                                                                                           cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGtt                                        >  2:3914680/1‑150 (MQ=255)
                                                                                                           cGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGtt                                        >  1:4360152/1‑150 (MQ=255)
                                                                                                            gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTg                                       >  1:2321919/1‑150 (MQ=255)
                                                                                                            gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTg                                       <  2:499185/150‑1 (MQ=255)
                                                                                                            gTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTg                                       <  2:3954806/150‑1 (MQ=255)
                                                                                                              gggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGaa                                     <  2:2447622/150‑1 (MQ=255)
                                                                                                              gggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGaa                                     >  2:2203675/1‑150 (MQ=255)
                                                                                                              gggAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGaa                                     >  2:1152460/1‑150 (MQ=255)
                                                                                                                gAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAAt                                   >  2:1130881/1‑150 (MQ=255)
                                                                                                                    aaaaGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATaaaa                               >  2:3270214/1‑150 (MQ=255)
                                                                                                                     aaaGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATaaaaa                              >  2:3894562/1‑150 (MQ=255)
                                                                                                                        gTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGc                           <  2:1366550/150‑1 (MQ=255)
                                                                                                                              gTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCAttttt                     <  1:1619974/150‑1 (MQ=255)
                                                                                                                               tGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTc                    >  1:4239215/1‑150 (MQ=255)
                                                                                                                                gCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTcc                   <  2:1387777/150‑1 (MQ=255)
                                                                                                                                  gTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGc                 <  2:3755104/150‑1 (MQ=255)
                                                                                                                                     cGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTg              >  1:2599868/1‑150 (MQ=255)
                                                                                                                                     cGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTg              <  2:792518/150‑1 (MQ=255)
                                                                                                                                          aaaCCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTc         >  1:3972616/1‑150 (MQ=255)
                                                                                                                                           aaCCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTcc        <  1:3429955/150‑1 (MQ=255)
                                                                                                                                              cTATTTCGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGa     >  1:3247377/1‑150 (MQ=255)
                                                                                                                                                 ttACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCTGCACAATGATTAAAAGATGATCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGcc  <  1:1008029/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATCCGGGGATGCTTTTGTCAGCACGCGATTATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGTTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCACAATGATTAAAAGATGAGCGGTTGAAATAAAAATGCATTTTTCCGCTTGTCTTCCTGAGCC  >  CP009273/219656‑219950

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: