Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,316 G→T 10.5% D233Y (GAT→TAT)  gpmM → phosphoglycero mutase III, cofactor‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,779,3160GT10.5% 419.4 / 6.9 180D233Y (GAT→TAT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base G (92/61);  new base T (3/24);  total (95/85)
Fisher's exact test for biased strand distribution p-value = 1.60e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

TCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCT  >  CP009273/3779170‑3779463
                                                                                                                                                  |                                                                                                                                                   
   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                                                                 >  2:1041975/1‑148 (MQ=255)
   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                                                                 <  2:657944/148‑1 (MQ=255)
       tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGGTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGACTTc                                                                                                                                             >  2:3285959/1‑148 (MQ=255)
        cGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCg                                                                                                                                            >  2:1747963/1‑148 (MQ=255)
              tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTTTGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                                                      >  2:1203035/1‑148 (MQ=255)
              tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                                                      >  2:821512/1‑148 (MQ=255)
              tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                                                      >  2:2603363/1‑148 (MQ=255)
               acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGc                                                                                                                                     >  1:3021102/1‑148 (MQ=255)
               acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGc                                                                                                                                     <  1:3296589/148‑1 (MQ=255)
                cgcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCg                                                                                                                                    <  2:3418168/148‑1 (MQ=255)
                 gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                                                                   <  2:1793695/148‑1 (MQ=255)
                 gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                                                                   >  1:1951599/1‑148 (MQ=255)
                 gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                                                                   <  2:3106184/148‑1 (MQ=255)
                    aTGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTAGTGAAATCGACCg                                                                                                                                >  1:1280209/1‑148 (MQ=255)
                     tGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGt                                                                                                                               >  2:469646/1‑148 (MQ=255)
                     tGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTAGGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGt                                                                                                                               <  2:3232381/148‑1 (MQ=255)
                             taTAACCGTTGGGATCGCGTAGAAAAAGCTTAGGACCTGCTGACTCTGGCGCAGGGCGAGTTCCCGGCCGATACCGCGGTTGGTGGTTTGCAGGCTGCTTTTGCTCGCGACTAAAATTTTTATTTCGTGAAAGCGACCGTTATCCGTg                                                                                                                       <  2:313529/147‑1 (MQ=255)
                               tacccGTGGGGACCGCGTAGACAACGCTTTGGCCGTGCTGACGCGGGCGCCGGGCGAGTGCCCGGCCGATACCGCCGTTGCTGGTTGCCAGGCTGCTTTGGCGCGCGCCGAAAATGATGAATTGGTGAAAGCGACCGTTATCCGTGCt                                                                                                                     <  2:2521341/145‑1 (MQ=255)
                                 aCCGTTGGTATCGCGTAGAAAAAGCTTAGGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGa                                                                                                                   <  2:2199546/148‑1 (MQ=255)
                                   ggtgggtaTCGCTTACAAAAAGCTTATCACCGGCTGACTCTGGCGCTGGGCGATTTCCAGGCCGATCCCGCCGTTGCTGGTTTGCCGGCTGATTATGCTCGCGACTAAAATGATGACTTCGTCAAAGCGACCGTTATCCGTGCTGAAg                                                                                                                 <  2:610678/141‑1 (MQ=255)
                                   cGTTGGGTTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCCGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAg                                                                                                                 <  2:2793560/148‑1 (MQ=255)
                                    gTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAgg                                                                                                                >  2:1146331/1‑148 (MQ=255)
                                     ttGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGCGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGt                                                                                                               <  2:785403/148‑1 (MQ=255)
                                     ttGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAATCGACCGTTATCCGTGCTGAAGGt                                                                                                               >  2:61573/1‑148 (MQ=255)
                                     ttGGGATCGCGTAGAAAAAGCTTATGACCGGCTGACTCGGGCGCAGGGCGAGTGCCAGGCCGATCCCGCCGTTGCGGGTTTGCCGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGt                                                                                                               <  2:2565810/148‑1 (MQ=255)
                                      tGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCCTTATACGTGCTGAAGGTc                                                                                                              >  1:3126170/1‑148 (MQ=255)
                                        ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAg                                                                                                            >  1:487790/1‑148 (MQ=255)
                                         gATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                           >  2:1170293/1‑148 (MQ=255)
                                         gATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                                                           >  1:2951929/1‑148 (MQ=255)
                                          aTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGcc                                                                                                          >  2:3511502/1‑148 (MQ=255)
                                           tCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                                                         >  2:2360407/1‑148 (MQ=255)
                                           tCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                                                         >  2:508229/1‑148 (MQ=255)
                                           tCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                                                         >  2:509133/1‑148 (MQ=255)
                                              cgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGAc                                                                                                      >  1:672401/1‑148 (MQ=255)
                                               ttataaaatGCTTATCACGGGCTGACTCGGGCGCTGGGCTTGTGCCCGGCCTATCCCGCGGTGGGTGTTTTGCTGGCTGCTTTTGCTCGCGACTAAATTTTTGATTTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcg                                                                                                     <  1:1331216/144‑1 (MQ=255)
                                                 aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                                                   >  1:1940233/1‑148 (MQ=255)
                                                 aGAAAAAGATTATGACCTGCTGACTGTGGCGCTGGGCGTGTGCCCGGCCGATACCGCGGTTGGTGTTTTGCAGGCTGATTATGCTCGCGACCAAATTGATTATTTCTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                                                   <  1:1343207/148‑1 (MQ=255)
                                                  gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                                                                  >  2:1299188/1‑148 (MQ=255)
                                                  gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                                                                  <  2:3168891/148‑1 (MQ=255)
                                                    aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGGAGGCTGCTTATGCTCGCGACGAAAATTGTGAATTCGTGGAAGCGACCGGTATCCGTGGGGAAGGTCAGCCAGACGCGGcc                                                                                                >  1:1634653/1‑147 (MQ=255)
                                                     aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                                                               >  1:2936483/1‑148 (MQ=255)
                                                     aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                                                               >  1:2492039/1‑148 (MQ=255)
                                                       aGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCACTg                                                                                             >  2:2958191/1‑148 (MQ=255)
                                                       aGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                                                             >  1:737231/1‑148 (MQ=255)
                                                         cTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                                                           >  1:899607/1‑148 (MQ=255)
                                                           taggACCTGCTGACTCGGGCGCCGGGCGAGTTCCCGGCCTATACCGCGGTTGCTGGTTTGCAGGGTGATTATGCTCGCGACTAAATTGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                                                         <  1:219739/145‑1 (MQ=255)
                                                           tATGACGTGCTGACTCGGGCGCTGGGCGTGTTCCCGGCCTATACCGCGGTTGCTGGTTTGCAGGGTGGTTATGCTCGCGACGAAATTTTTTAATTGTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                                                         <  1:1374421/148‑1 (MQ=255)
                                                           tATGACCTGCGTACTCTGGCGCTGGGCGTGTTCCGGGCCTATCCCGCCGTTGGTGGTTTGCTGGCTGCTTATGCTCGCGACTAAATTGTTGAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                                                         <  1:219710/148‑1 (MQ=255)
                                                            aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                                                        <  2:1951599/148‑1 (MQ=255)
                                                            aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                                                        >  2:3295035/1‑148 (MQ=255)
                                                              gACCTGCTGACTGGGGCGCCGGGCGTGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTAGGCTCGCGACTAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                                                      <  2:1252141/148‑1 (MQ=255)
                                                              gACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTCTGCTCGCTACTAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                                                      <  2:3024062/148‑1 (MQ=255)
                                                               ccgtgctgAGTCTGGCTCAGGGCGAGTTCCTGGCCGATACCGCCGTTGCTGGTTTGCTTGCTGCTTATGCACGCGTCTAAAATGATGAATTCGTGATAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACgg                                                                                     <  2:2578419/145‑1 (MQ=255)
                                                               aCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACgg                                                                                     >  2:2237515/1‑148 (MQ=255)
                                                                    ctgACTCGGGCGCTGGGGGTGTTCCCGCCCTATACCGCCTTGGGTGTTTTGCTGGCTGATTTTGCTCGCCACTAATTTTTTTATTTGGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                                                                                <  1:359192/148‑1 (MQ=255)
                                                                       aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCCGCCAGACGCGGCCATGGAAGACGGCGATGCAc                                                                             >  1:3041216/1‑148 (MQ=255)
                                                                           gggCGCTGGGCTTGTTCCCGCCCGATCCCGCGTTGGGGGGTTTGCGGGGTGCTTCTGCTCGCGACCAAAATGTTTATTTCTTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAt                                                                         <  1:1660971/147‑1 (MQ=255)
                                                                              cgcTGGGCGAGTCCCCGGCCGATCCCGCGGTTGGTGTTTTGCTGGCTGCTTTTGCGCGCGACCAAATTTATTATTTCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                                                      <  1:1577533/148‑1 (MQ=255)
                                                                               tctgggCGAGTCCCTGCCCTTAACCGCGGTTGCTGTTTTGCTGTGTGATTTGGCTCGCTAGGAAATTGTTTAAGCCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTc                                                                     <  1:2814234/145‑1 (MQ=255)
                                                                                cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                                                    >  2:3398161/1‑148 (MQ=255)
                                                                                 aGGGCGAGTTCCCGCCCTTTACCGCGGTTGGTGGTTTGCTGGGTGCTTTTGCTCGCGACTAAAATGTTGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCAt                                                                   <  1:2585525/148‑1 (MQ=255)
                                                                                 aGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCCGACGCGGCAAAGGAAGACGGCGATGCACAGATTTTAAt                                                                   >  1:1868762/1‑148 (MQ=255)
                                                                                   ggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCGTTTTCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGGACTGATTTTCATGa                                                                 >  1:977431/1‑148 (MQ=255)
                                                                                   ggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAAGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                                                                 >  1:3413113/1‑148 (MQ=255)
                                                                                     cGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTTGAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAAc                                                               >  1:3318614/1‑148 (MQ=255)
                                                                                               gccgtgcccgcggTGGCTGGTTTGCGGGGTGATTATGCTCGCGACTAAATTTATGACTGGGTCACAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                     <  1:1003905/141‑1 (MQ=255)
                                                                                               gCCGAGACCGCGGTTGGGGTTTTGCTGGCTGCTTTTGCGCGCGACGAATTTTATGATTTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                                                     <  1:1041975/148‑1 (MQ=255)
                                                                                                  gATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCGGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAcc                                                  >  2:899047/1‑148 (MQ=255)
                                                                                                          ggTTGCTGGTTTGCCGGCTGCTTTTGCTCGCGACGAAAATGATGAATGGGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt                                          <  1:821512/147‑1 (MQ=255)
                                                                                                          cgTTGCTGGTTTGCTGGCTGCTTATGCGCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt                                          <  1:2603363/148‑1 (MQ=255)
                                                                                                          cgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt                                          >  2:2836934/1‑148 (MQ=255)
                                                                                                           gTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTg                                         >  2:2765987/1‑148 (MQ=255)
                                                                                                                 ggTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCa                                   >  1:1471605/1‑148 (MQ=255)
                                                                                                                  gTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCAc                                  >  1:1443327/1‑148 (MQ=255)
                                                                                                                   tttGCAGGCTGCTTATGCTCGGGACGAAATTGATGATTTGGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACt                                 <  2:1234990/148‑1 (MQ=255)
                                                                                                                    ttGCAGGCTGCTTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTc                                <  2:1801437/148‑1 (MQ=255)
                                                                                                                       cAGGCTGATTATGCTCGCGACTAAAATGTTTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTg                             <  1:225450/148‑1 (MQ=255)
                                                                                                                        aGGCTGCTTATGCTCGCGACTAAATTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGc                            <  1:3169732/148‑1 (MQ=255)
                                                                                                                            tgctTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTc                        <  1:2052534/148‑1 (MQ=255)
                                                                                                                            tgctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTc                        <  1:1160886/148‑1 (MQ=255)
                                                                                                                             tctTATGCTCGCTAATAAAATGATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCg                       <  1:3511502/147‑1 (MQ=255)
                                                                                                                               ttATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTg                     >  1:411538/1‑148 (MQ=255)
                                                                                                                                tATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGa                    <  1:176790/148‑1 (MQ=255)
                                                                                                                                 aTGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGaa                   >  1:1333050/1‑148 (MQ=255)
                                                                                                                                   gCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACg                 >  1:303384/1‑148 (MQ=255)
                                                                                                                                     tCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCt               <  1:3308832/148‑1 (MQ=255)
                                                                                                                                      cgcgACTAAAATTATGAATGCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTg              <  1:3285959/148‑1 (MQ=255)
                                                                                                                                      cgcgACGAAAATTATGAATTGGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTg              <  1:1010428/148‑1 (MQ=255)
                                                                                                                                        cgacgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAt            >  2:2572637/1‑148 (MQ=255)
                                                                                                                                        cgacgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAt            <  2:885670/148‑1 (MQ=255)
                                                                                                                                          acgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGAttt          >  1:747493/1‑148 (MQ=255)
                                                                                                                                             aaaaTTATCAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa       <  1:390730/148‑1 (MQ=255)
                                                                                                                                             aaaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa       >  2:1637447/1‑148 (MQ=255)
                                                                                                                                             aaaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGa       <  1:1167004/148‑1 (MQ=255)
                                                                                                                                               aaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATg     >  2:3126946/1‑148 (MQ=255)
                                                                                                                                               aaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATg     >  2:1356899/1‑148 (MQ=255)
                                                                                                                                                atgatgAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATgg    >  2:2508340/1‑148 (MQ=255)
                                                                                                                                                atgatgAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATgg    >  2:727638/1‑148 (MQ=255)
                                                                                                                                                 tgatgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGc   <  2:1776917/148‑1 (MQ=255)
                                                                                                                                                 tgaagaATGCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGc   <  2:1102995/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCT  >  CP009273/3779170‑3779463

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: