Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,281 C→G 8.1% 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.1% 584.0 / 6.7 209A221G (GCT→GGT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (106/82);  new base G (4/9);  total (110/99)
Fisher's exact test for biased strand distribution p-value = 8.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

AAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAT  >  CP009273/3779133‑3779428
                                                                                                                                                     |                                                                                                                                                   
aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                    >  1:1440529/1‑150 (MQ=255)
   tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGt                                                                                                                                                 >  2:3897245/1‑150 (MQ=255)
    ttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                                >  1:113811/1‑150 (MQ=255)
    ttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                                <  2:1013034/150‑1 (MQ=255)
     tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                               >  1:2128686/1‑150 (MQ=255)
      gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                              <  1:628317/150‑1 (MQ=255)
        cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa                                                                                                                                            <  2:916992/150‑1 (MQ=255)
            cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                        <  1:3173063/150‑1 (MQ=255)
             tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                       >  2:2529729/1‑150 (MQ=255)
              gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                      >  1:1164287/1‑150 (MQ=255)
              gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                      <  2:3195560/150‑1 (MQ=255)
                 cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                   <  1:2890853/150‑1 (MQ=255)
                  aaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                  <  1:3269349/150‑1 (MQ=255)
                  aaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGGTTGCAGGCTGCGTAt                                                                                                                                  >  1:1565636/1‑150 (MQ=255)
                    aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                                <  1:814390/150‑1 (MQ=255)
                     ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCACGCTGCGTaatcg                                                                                                                               >  1:1103163/1‑146 (MQ=255)
                     ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGCGGGTTGGCAGGCGGCGTATGcg                                                                                                                               >  1:1103127/1‑149 (MQ=255)
                           gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgac                                                                                                                         >  2:3609380/1‑150 (MQ=255)
                            tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacg                                                                                                                        <  1:826063/150‑1 (MQ=255)
                             aGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacga                                                                                                                       <  2:1779308/150‑1 (MQ=255)
                               cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaa                                                                                                                     <  2:4165606/150‑1 (MQ=255)
                                gTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaa                                                                                                                    <  2:2475424/150‑1 (MQ=255)
                                  ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatg                                                                                                                  >  1:667378/1‑150 (MQ=255)
                                   catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                                 >  2:488129/1‑150 (MQ=255)
                                     tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgatg                                                                                                               <  1:601579/150‑1 (MQ=255)
                                       atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCGTATTCTCGCGACGAAAATGATGaa                                                                                                             >  1:4053171/1‑150 (MQ=255)
                                       atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGACAATGATGaa                                                                                                             >  1:2421923/1‑150 (MQ=255)
                                       atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCAGGTTTGCATGCTGCTTATGCTCGCGACGAAAATGATTaa                                                                                                             >  2:1323431/1‑150 (MQ=255)
                                        ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAt                                                                                                            <  2:2748091/150‑1 (MQ=255)
                                          ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTc                                                                                                          <  1:1438745/150‑1 (MQ=255)
                                            tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGt                                                                                                        <  1:3007827/150‑1 (MQ=255)
                                              gTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGa                                                                                                      >  2:2882024/1‑150 (MQ=255)
                                              gTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGa                                                                                                      >  1:873828/1‑150 (MQ=255)
                                               ttACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaa                                                                                                     >  2:1478994/1‑150 (MQ=255)
                                                tactacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaaa                                                                                                    <  2:503684/150‑1 (MQ=255)
                                                  ctacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGc                                                                                                  >  1:805891/1‑150 (MQ=255)
                                                   tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCg                                                                                                 <  1:2764692/150‑1 (MQ=255)
                                                           ggACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTAt                                                                                         >  1:3412468/1‑150 (MQ=255)
                                                             aCCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATcc                                                                                       >  2:1716041/1‑150 (MQ=255)
                                                              ccGCTATAACGTTGGGTATCGCGTATAAATAGCTTATGACCTGATGACTCGGGGGCTGGGCGTGTTCCAGGCGTATACCGCGTTTGGTGGTTTGCTGGGTGCTTTTGCTCGCGACTAAATTGTTTAATTCTTGAAAGCGACCGTTATCCg                                                                                      <  2:778838/150‑1 (MQ=255)
                                                                 cgATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGc                                                                                   >  2:1904651/1‑150 (MQ=255)
                                                                  gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCt                                                                                  >  1:1332837/1‑150 (MQ=255)
                                                                   aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTg                                                                                 >  2:3558279/1‑150 (MQ=255)
                                                                     aacgtttGGGATCGCTTTGAAAAAGCTTTTGACCTGCTGACTCGGGCGCTGGGCTAGTTCCCGGCCTATCCCGCGGTTGCTGGTTTTCAGGCTGATTTTGCTCGCGACTAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGaa                                                                               <  1:1760465/145‑1 (MQ=255)
                                                                     aaCCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAATCGACCGTTATCCGTGCTGaa                                                                               >  1:2224472/1‑150 (MQ=255)
                                                                      ccggTTGGGATCGCGTATAAAAAGCTTATGCCCTGCTGACTCGGGCGCCGGGCGTGTTCCCGGCCGATCCCGCGGTTGCTGGTTTGCAGGGTGCTTATGCTCGCGACTAAATTTATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAg                                                                              <  2:2947426/147‑1 (MQ=255)
                                                                      aCCGTTGGTATCGCGTAGAAAAAGCTTATGACGTGCTGACTCTGGCGCTGGGCTAGTTCCAGGCCGATACCGCGGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAg                                                                              <  2:2556824/150‑1 (MQ=255)
                                                                         gtggggATCGCGTATAAAACGGTTTTGACCTGCGGACTCGGGCGCTGGGCGTGTTCCCGGCCGATACCGCCGTTGCGGGTTTGCTGGCTGCTTAGGCTCGCGACGAAAATTATGAATTGGTGAAAGCGACCGTTATCCGTGCTGAAGGTc                                                                           <  1:3241904/147‑1 (MQ=255)
                                                                          ttGGGATCGCTTAGAAAATGATTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCACGCGACTAAAATGATTATTTCGTTAACGCGACCGTTATCCGTGCTGAAGGTCa                                                                          <  1:2447340/150‑1 (MQ=255)
                                                                           tGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAg                                                                         <  2:225078/150‑1 (MQ=255)
                                                                           tGGGATCGCGTAGAAAAAGCTGATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCCGGCCGATCCCGCGGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGCCGAAAATGATGAATTGGTGAACGCGACCGTTATCCGTGCTGAAGGTCAg                                                                         <  2:2402119/150‑1 (MQ=255)
                                                                            gggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                        >  2:3409976/1‑150 (MQ=255)
                                                                                aTCGCTTATAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCTGGCCTATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAg                                                                     <  2:1600475/150‑1 (MQ=255)
                                                                                  cgcgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTATGCTCGCGACTAAAATGATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                                                   <  2:2237267/150‑1 (MQ=255)
                                                                                   gcgTAGAAACCGCTTAGGACCTGCTGACTCGGGCGCCGGGCGAGTTCCCGGCCTATACCGCGGTTGGTGGTTTGCAGGGTGCTTATGCTCGCGACGAAAATGATGAATTGGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcg                                                                  <  2:3292813/150‑1 (MQ=255)
                                                                                          aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTCCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATTATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                           <  2:4042071/150‑1 (MQ=255)
                                                                                             aGCTTATGACGTGCTGACTCTGGCGCTGGGCGTGTTCCCGGCCTATCCCGCGGTTGTGGTTTTGCCGGCTGCTTATGCTCGCCACGAAAATTATTAATTCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                        <  1:2700502/150‑1 (MQ=255)
                                                                                              gtttATGACCTGCGGACTCTGGCGCTGGGCGTGTTCCAGGCCGTTACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACGAAATTTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                       <  1:334224/148‑1 (MQ=255)
                                                                                                ttattaCCTGCTGTCTCTGGCGCTGGGCGATTTCCTGGCCTATACCGCTTTTGCTGGTTTGATGTTTGCTTTTGCTCGCGACTAAATTTTTGAATTCTTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                     <  1:220460/150‑1 (MQ=255)
                                                                                                ttATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCACAGGAAGa                                                     >  1:435708/1‑150 (MQ=255)
                                                                                                ttATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                     >  2:3169226/1‑150 (MQ=255)
                                                                                                     aCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg                                                >  2:1054661/1‑150 (MQ=255)
                                                                                                       gtgctgACTGGGGCGCTGGGGGTGTTCCTGGCGGATCCCGCGGTGGCTGGTTTGCAGGGTGATTTGGCTCGCGACTAAATTGTTGAATTGGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                              <  1:3889816/149‑1 (MQ=255)
                                                                                                       ctgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                              >  1:4147384/1‑150 (MQ=255)
                                                                                                           tgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTTTTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCGGAAGGTCAGCCAGCCGCGGCAATGGAAGACCGCGATGCAc                                          >  1:3031960/1‑150 (MQ=255)
                                                                                                             aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGCTGCACTg                                        >  1:1951082/1‑150 (MQ=255)
                                                                                                             aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTg                                        <  2:1055825/150‑1 (MQ=255)
                                                                                                             aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGCTGGTTTTGAGGGGGCGTATGCTCGCGCCGCAAATTATGACTTAGTGAAAGCCACCCGTATCCGTGCGGAAGGTCCGCCCCACGCGGCCAAGGAAGACCGCGATGCACTg                                        >  1:4078845/1‑150 (MQ=255)
                                                                                                              ctctGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGa                                       >  2:4352969/1‑150 (MQ=255)
                                                                                                              ctctGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACCCGGCAATGGAAGACGGCGATGCACTGa                                       >  1:494790/1‑150 (MQ=255)
                                                                                                                  ggCGCTGGGCTAGTTCCAGGCCGATACCGCCGTTGATGTTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                   <  2:517747/150‑1 (MQ=255)
                                                                                                                  ggCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                   >  2:3173346/1‑150 (MQ=255)
                                                                                                                      cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATg                               <  2:1252201/150‑1 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGGGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                              >  2:770476/1‑150 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCAGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCCTTATACGTGCTGAAGGTCAGCCCGACGCGGCAAAGGAAGAAGGCGATGCACAGATTTTAAtta                              >  1:2330736/1‑148 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGTAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACTGATTTTCATGa                              >  1:121093/1‑150 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                              >  2:118525/1‑150 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACGTATTTTAATGa                              >  1:2782883/1‑150 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAATCGTGAAAGCGACCCTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATTCACTGATTTTAATGa                              >  1:2782917/1‑150 (MQ=255)
                                                                                                                         ggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAAc                            >  2:984893/1‑150 (MQ=255)
                                                                                                                          gggAGTTCCTGGCCGTTCCCGCCGTTTGTGGTTTGCAGGCTGATTTTGGTCGCGATTAAATTGATTAATTGTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGTCGGCGATGCACTGATTTTCATGAACt                           <  1:3084651/148‑1 (MQ=255)
                                                                                                                          gCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt                           >  2:1561030/1‑150 (MQ=255)
                                                                                                                          gCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAGGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt                           <  2:1440529/150‑1 (MQ=255)
                                                                                                                           cttgTTCCAGGCCGATCCCGCGGTTGGTGTTTTGCTGGCTTCTTATGCTCGCGACGAAATTTTTGAATGCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACtt                          <  2:3312227/147‑1 (MQ=255)
                                                                                                                              gTTCCAGGCCTATACCGCGGTTGGTGGTTTGCAGGCTGCTTATGCTCGCGACTAAATTTATTATTTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg                       <  1:1637666/150‑1 (MQ=255)
                                                                                                                              gTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACGGATTTTCATGAACTTCCg                       >  2:3044108/1‑150 (MQ=255)
                                                                                                                               ttCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGt                      >  2:1590806/1‑150 (MQ=255)
                                                                                                                                 ccAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGc                    >  2:267816/1‑150 (MQ=255)
                                                                                                                                  cAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAATGTCAGCCACACGCGGCAATGGAAGACGGCGATGCACAGATTTTCATTAACATCCGTGCt                   >  1:3826332/1‑150 (MQ=255)
                                                                                                                                   tggCCTATACCGGCGTTGTTGGTTGGCTGGGTGCTTATGCTCGCGACTAAATTTTTTAATTCTTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                  <  1:2955139/149‑1 (MQ=255)
                                                                                                                                   aGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                  >  1:100871/1‑150 (MQ=255)
                                                                                                                                   aGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                  >  2:2440804/1‑150 (MQ=255)
                                                                                                                                     gCCGATCCCGCCGTTGCTGGTTTGCAGGCTGCTTAGGCTCGCGACGAAAATGATGAATGCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAc                <  2:1906797/150‑1 (MQ=255)
                                                                                                                                        gATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAATCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgc             >  2:36253/1‑150 (MQ=255)
                                                                                                                                         aTACCGCCGTTGCTGTTTTGCTGGCTGATTATGCTCGCGACGAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcg            <  2:113811/150‑1 (MQ=255)
                                                                                                                                         aTACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcg            >  2:2796010/1‑150 (MQ=255)
                                                                                                                                              gccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt       >  1:1559138/1‑150 (MQ=255)
                                                                                                                                                ggTTGCTGGTTTGCTGGCGGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGa     <  2:4129948/149‑1 (MQ=255)
                                                                                                                                                  ttGCTGGTTGGCTGGGTGGTTTTGCTCGCGACGAAATTGTTGATTTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaaa   <  2:3283810/150‑1 (MQ=255)
                                                                                                                                                   tGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAt  <  1:3609380/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAT  >  CP009273/3779133‑3779428

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: