Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,281 C→G 8.5% A221G (GCT→GGT)  gpmM → phosphoglycero mutase III, cofactor‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,779,2810CG 8.5% 548.9 / 5.5 194A221G (GCT→GGT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (94/80);  new base G (3/17);  total (97/97)
Fisher's exact test for biased strand distribution p-value = 1.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

ATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAA  >  CP009273/3779135‑3779427
                                                                                                                                                  |                                                                                                                                                  
aTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGcg                                                                                                                                                   >  1:610678/1‑147 (MQ=255)
 tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                  >  2:2571079/1‑148 (MQ=255)
 tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                  >  2:115160/1‑148 (MQ=255)
    gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                               <  1:2775030/148‑1 (MQ=255)
      cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                             <  1:75131/148‑1 (MQ=255)
        cgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa                                                                                                                                           >  1:814537/1‑148 (MQ=255)
          cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAgg                                                                                                                                         <  1:1122908/148‑1 (MQ=255)
          cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAgg                                                                                                                                         >  2:1005033/1‑148 (MQ=255)
             ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCTAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                      <  2:3068931/148‑1 (MQ=255)
             ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                      >  2:225450/1‑148 (MQ=255)
              gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                     <  1:409623/148‑1 (MQ=255)
              gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                     <  1:3505496/148‑1 (MQ=255)
               cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct                                                                                                                                    <  2:1711515/148‑1 (MQ=255)
               cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct                                                                                                                                    <  2:2888973/148‑1 (MQ=255)
               cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct                                                                                                                                    >  1:3413823/1‑148 (MQ=255)
                 aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                  <  1:1253200/148‑1 (MQ=255)
                 aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                  <  2:2222558/148‑1 (MQ=255)
                  aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                 >  2:1836544/1‑148 (MQ=255)
                  aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                 <  1:1000776/148‑1 (MQ=255)
                  aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                 <  1:1000754/148‑1 (MQ=255)
                  aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                 >  2:494166/1‑148 (MQ=255)
                   ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATg                                                                                                                                >  1:1088792/1‑148 (MQ=255)
                    gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                               <  2:3254889/148‑1 (MQ=255)
                     ccGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCt                                                                                                                              >  2:3136017/1‑148 (MQ=255)
                         gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcgcg                                                                                                                          <  1:3121095/148‑1 (MQ=255)
                          tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcga                                                                                                                         <  2:1423494/148‑1 (MQ=255)
                             cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTNCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTTGAGGCTGCGTATTCTAGcgacga                                                                                                                      >  1:1801437/1‑148 (MQ=255)
                                ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaa                                                                                                                   <  1:1481430/148‑1 (MQ=255)
                                ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaa                                                                                                                   <  1:1983472/148‑1 (MQ=255)
                                ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaa                                                                                                                   >  2:2331372/1‑148 (MQ=255)
                                   tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                                <  1:3410187/148‑1 (MQ=255)
                                      ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                             >  2:1041975/1‑148 (MQ=255)
                                      ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                             <  2:657914/148‑1 (MQ=255)
                                      ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                             <  2:657973/148‑1 (MQ=255)
                                         gTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAtt                                                                                                          <  1:1760865/148‑1 (MQ=255)
                                          tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGGTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGACTTc                                                                                                         >  2:3285959/1‑148 (MQ=255)
                                           cGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCg                                                                                                        >  2:1747963/1‑148 (MQ=255)
                                             ttACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATCCGTg                                                                                                      >  2:850483/1‑148 (MQ=255)
                                                 tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTTTGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                  >  2:1203035/1‑148 (MQ=255)
                                                 tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                  >  1:1980685/1‑148 (MQ=255)
                                                  acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGc                                                                                                 <  1:3296589/148‑1 (MQ=255)
                                                    gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                               <  2:3106184/148‑1 (MQ=255)
                                                    gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                               <  2:1793695/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)
                                                         ggACCGCTATACCCTTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTGTGGCGCAGGGCTAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGtt                                                                                          <  2:3400680/148‑1 (MQ=255)
                                                         ggACCGCTATAACGGTTGGGATCGCGTATAAAAAGATTATGACCTGGTGACTCTGGCGCTGGGCGTGTTCCAGGCCGATACCGCCGTTGGTGTTTTGCTGGCTGCTTTTGCTCGCGACTAAAATGTTGAATTGGTGAAAGCGACCGtt                                                                                          <  1:2338814/148‑1 (MQ=255)
                                                              gcgATACCCGTTGGGATCGCGTAGAAAACGCTTATGACGTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCTATACCGCGGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTTAAAGCGACCGTTATCCg                                                                                     <  1:2965663/148‑1 (MQ=255)
                                                                taTAACCGTTGGGATGGCTTAGAAAAAGCTTATGACCTGCTGACTGTGGCGCAGGGCGTGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTAGGCGCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTg                                                                                   <  2:3358548/147‑1 (MQ=255)
                                                                 aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAATCGACCGTTATCCGTGc                                                                                  >  1:1010917/1‑148 (MQ=255)
                                                                  tAACGGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGGGTTCCCGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGGTTATGCTCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCt                                                                                 <  2:2521376/148‑1 (MQ=255)
                                                                    aCCGTTGGTATCGCGTAGAAAAAGCTTAGGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGa                                                                               <  2:2199546/148‑1 (MQ=255)
                                                                      cGTTGGGTTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCCGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAg                                                                             <  2:2793560/148‑1 (MQ=255)
                                                                        ttGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGt                                                                           >  2:61568/1‑148 (MQ=255)
                                                                         tGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTc                                                                          >  1:1356553/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:509133/1‑148 (MQ=255)
                                                                              tCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                     >  2:2360407/1‑148 (MQ=255)
                                                                                    aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                               >  1:1940290/1‑148 (MQ=255)
                                                                                     gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTTCTGGTTTTGAGGCTGCTTATGCTCGCGACGAAAATGATTCAATAGTGAAATCGACCCTTATAAGTGCTGAATGTAAGCCAGACGCgg                                                              >  1:3124104/1‑148 (MQ=255)
                                                                                     gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                              >  2:1299188/1‑148 (MQ=255)
                                                                                     gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                              <  2:3168891/148‑1 (MQ=255)
                                                                                        aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                           >  1:2936483/1‑148 (MQ=255)
                                                                                        aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                           >  2:113971/1‑148 (MQ=255)
                                                                                        aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                           >  1:2492039/1‑148 (MQ=255)
                                                                                         aaGCTTATGACCGGCTGCCGCGGGCGCGGGGCGAGTGCCCGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTAGGCTCGCGACTAAAATTATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                          <  2:758612/148‑1 (MQ=255)
                                                                                          aGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCACTg                                                         >  2:2958191/1‑148 (MQ=255)
                                                                                            cTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                       <  2:814537/148‑1 (MQ=255)
                                                                                            cTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                       >  1:1909011/1‑148 (MQ=255)
                                                                                             ttATGACCGGCTCACTCGGGCGCCGGGCGATTTCCCGCCCTATACGGCGGTTGCTGTTTTGCTGGCTGTTTTGGCTCGCTACTAAAATTTTTATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                      <  1:115160/148‑1 (MQ=255)
                                                                                              taggACCTGCTGACTCGGGCGCCGGGCGAGTTCCCGGCCTATACCGCGGTTGCTGGTTTGCAGGGTGATTATGCTCGCGACTAAATTGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                     <  1:219739/145‑1 (MQ=255)
                                                                                               aggACCTGCTGCATCTGGGGCAGGGCGTGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCTGGGTGCTTATGATCGGGACTAAAATGATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                    <  1:1183928/146‑1 (MQ=255)
                                                                                               aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                    >  2:3295035/1‑148 (MQ=255)
                                                                                               aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATCCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                    <  1:850483/148‑1 (MQ=255)
                                                                                                 gACCTGCTGACTGGGGCGCCGGGCGTGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTAGGCTCGCGACTAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                  <  2:1252141/148‑1 (MQ=255)
                                                                                                 gACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTCTGCTCGCTACTAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                  <  2:3024062/148‑1 (MQ=255)
                                                                                                       ctgACTCGGGCGCTGGGGGTGTTCCCGCCCTATACCGCCTTGGGTGTTTTGCTGGCTGATTTTGCTCGCCACTAATTTTTTTATTTGGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                                            <  1:359192/148‑1 (MQ=255)
                                                                                                          aCTCTGGCGCCGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGGTGCTTTTGCTCGCGACGAAAATGATGAATTCGTGAAACCGACCGGTATACGTGCTGAAGGTCAGCCATACGCGGCAATTGAATAAGGCGATGCAc                                         >  1:2078620/1‑148 (MQ=255)
                                                                                                                gcgcAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATACGTGCGGAAGGTCAGCCAGACGCGGCCATGGAAGAAGGCGATGCACTGAttt                                   >  1:1340294/1‑148 (MQ=255)
                                                                                                                 cgcTGGGCGAGTCCCCGGCCGATCCCGCGGTTGGTGTTTTGCTGGCTGCTTTTGCGCGCGACCAAATTTATTATTTCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                  <  1:1577533/148‑1 (MQ=255)
                                                                                                                  tctgggCGAGTCCCTGCCCTTAACCGCGGTTGCTGTTTTGCTGTGTGATTTGGCTCGCTAGGAAATTGTTTAAGCCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTc                                 <  1:2814234/145‑1 (MQ=255)
                                                                                                                   cAGGGCGAGTTCCCGGCCGATACCGCGGTTGGTGTTTTGCCGGCTGCTTATGCTCGCGACGAAAATTTTGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                <  1:1432485/148‑1 (MQ=255)
                                                                                                                    aGGGCGAGTTCCCGCCCTTTACCGCGGTTGGTGGTTTGCTGGGTGCTTTTGCTCGCGACTAAAATGTTGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCAt                               <  1:2585525/148‑1 (MQ=255)
                                                                                                                     gggCGAGTTCCAGGCCGATACCGCCGTTGGGGGTTTTCAGGCTGCGTATTCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCGGAAGGTCAGCCCGACGCGGCAAAGGAAGACGGCGTTTCACTGATTTTAATg                              >  1:1604363/1‑148 (MQ=255)
                                                                                                                      ggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCGTTTTCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGGACTGATTTTCATGa                             >  1:977431/1‑148 (MQ=255)
                                                                                                                       gCGAGTTCCCGGCCTATACCGCGTTTGCGGTTTTGCTGGCTGTTTATGCTCGCGACGAAATTTTTTAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                            <  1:3024895/148‑1 (MQ=255)
                                                                                                                       gCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAATATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACTGATTTTCATGaa                            >  1:1649784/1‑148 (MQ=255)
                                                                                                                        cGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTTGAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAAc                           >  1:3318614/1‑148 (MQ=255)
                                                                                                                          aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACtt                         <  1:3136017/148‑1 (MQ=255)
                                                                                                                                      atcccgccgTTGGTGGTTGGCTGGCTGATTATGCTCGCGACTAAAATGTTGAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcg             <  1:709743/145‑1 (MQ=255)
                                                                                                                                             ggTTGGTGGTTTGCTGGCTGCTTATGCTCGCGACGAAATTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt      <  1:3484677/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)
                                                                                                                                                tGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa   <  2:2595792/148‑1 (MQ=255)
                                                                                                                                                tGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa   >  2:1414955/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAA  >  CP009273/3779135‑3779427

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: