Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,281 C→G 5.8% 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 5.8% 839.2 / 2.3 270A221G (GCT→GGT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (126/125);  new base G (4/14);  total (130/140)
Fisher's exact test for biased strand distribution p-value = 2.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

AAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATC  >  CP009273/3779133‑3779429
                                                                                                                                                    |                                                                                                                                                    
aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                     >  1:31358/1‑150 (MQ=255)
aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                     <  1:5129506/150‑1 (MQ=255)
 aaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                    >  2:1818520/1‑150 (MQ=255)
   tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGt                                                                                                                                                  >  2:2114244/1‑150 (MQ=255)
    ttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGtt                                                                                                                                                 >  2:3526381/1‑150 (MQ=255)
     tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGTTTACCGTTTGGATCGCGTAGAAAAATCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                <  1:3334153/150‑1 (MQ=255)
     tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                <  2:2442373/150‑1 (MQ=255)
       ccGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGc                                                                                                                                              <  2:5012625/150‑1 (MQ=255)
       ccGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGc                                                                                                                                              <  1:521818/150‑1 (MQ=255)
        cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa                                                                                                                                             <  2:1137271/150‑1 (MQ=255)
           gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc                                                                                                                                          >  2:4226203/1‑150 (MQ=255)
            cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGATTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                         <  1:2921594/150‑1 (MQ=255)
            cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                         <  1:2223765/150‑1 (MQ=255)
            cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                         >  2:5293056/1‑150 (MQ=255)
             tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                        >  1:5189929/1‑150 (MQ=255)
             tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                        >  1:1775636/1‑150 (MQ=255)
             tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                        >  2:3113764/1‑150 (MQ=255)
              gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGgctgc                                                                                                                                       >  1:4924662/1‑150 (MQ=255)
               ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct                                                                                                                                      <  2:293299/150‑1 (MQ=255)
                 cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                    >  2:1089090/1‑150 (MQ=255)
                 cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                    <  2:1328667/150‑1 (MQ=255)
                  aaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                   >  1:1171040/1‑150 (MQ=255)
                  aaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                   <  2:2158483/150‑1 (MQ=255)
                   aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATg                                                                                                                                  >  2:4558371/1‑150 (MQ=255)
                     ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCt                                                                                                                                <  2:1301404/150‑1 (MQ=255)
                      gCCGCGTAGCGTCCATCAGTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAACAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTTTGCTc                                                                                                                               <  2:4007519/150‑1 (MQ=255)
                        cgcgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcgc                                                                                                                             <  1:742051/150‑1 (MQ=255)
                         gcgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcgcg                                                                                                                            <  1:376735/150‑1 (MQ=255)
                          cgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcga                                                                                                                           <  2:2846008/150‑1 (MQ=255)
                           gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgac                                                                                                                          >  2:545020/1‑150 (MQ=255)
                            tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacg                                                                                                                         <  2:3402763/150‑1 (MQ=255)
                                    atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgat                                                                                                                 <  1:348803/150‑1 (MQ=255)
                                     tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgatg                                                                                                                <  2:4551240/150‑1 (MQ=255)
                                         tGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAtt                                                                                                            <  1:3622936/150‑1 (MQ=255)
                                         tGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAtt                                                                                                            >  2:2728809/1‑150 (MQ=255)
                                         tGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAtt                                                                                                            >  2:4577579/1‑150 (MQ=255)
                                          ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTc                                                                                                           <  1:2115189/150‑1 (MQ=255)
                                          ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTc                                                                                                           >  2:2279525/1‑150 (MQ=255)
                                                    acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                                 <  1:2925485/150‑1 (MQ=255)
                                                    acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGa                                                                                                 <  1:1896941/150‑1 (MQ=255)
                                                      gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAcc                                                                                               <  1:4427560/150‑1 (MQ=255)
                                                       cgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAATCGACCg                                                                                              >  1:467846/1‑150 (MQ=255)
                                                        gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCCACGAAAATGATTAATTCGTGAAAGCGAccct                                                                                             >  1:65920/1‑148 (MQ=255)
                                                        gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGt                                                                                             >  2:3909565/1‑150 (MQ=255)
                                                        gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTAGCGACGAAAATGATGAATTCGTGAAACCGACCGt                                                                                             >  1:1608432/1‑150 (MQ=255)
                                                           ggACCGCGATAACGGTTGGTATCGCGTAGAAAATGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCTGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACTAAATTGATTAATTCGTGAAAGCGACCGTTAt                                                                                          <  2:4165369/150‑1 (MQ=255)
                                                            gACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGGTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATc                                                                                         <  2:3478488/150‑1 (MQ=255)
                                                                 cgATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGc                                                                                    >  2:5088988/1‑150 (MQ=255)
                                                                  gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCt                                                                                   >  2:2785714/1‑150 (MQ=255)
                                                                   aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTg                                                                                  >  2:3560483/1‑150 (MQ=255)
                                                                   aTAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTg                                                                                  >  1:1271161/1‑150 (MQ=255)
                                                                    tAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCATATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGa                                                                                 >  1:1030634/1‑150 (MQ=255)
                                                                       ccGTTGGGATCGCGTATAAAATGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAgg                                                                              <  2:1611863/150‑1 (MQ=255)
                                                                         gTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCATGCTGCTTATGCTCGCGACGAAAATGATGAATTATTTAAAGCGACCGTTATCCGTGCTGAAGGTc                                                                            >  2:252932/1‑150 (MQ=255)
                                                                         gTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCGGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTc                                                                            <  2:4892340/150‑1 (MQ=255)
                                                                          ttGGTATCGCTTATAAAAAGCTTATTACTTGCTTAATCTGGCGCAGGGCTTGTTCCTGGCCGATACCGCGTTTGCTGGTTTGCTGGTTTATTATGCTCGCGACGAAATTTTTGAATTGGTGAAAGCGACCGTTATCCGTGCTGAAGGTCa                                                                           <  1:4387542/150‑1 (MQ=255)
                                                                          ttGGTATCGCGTATAAAAAGCTTTTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTATGCTCGCGACTAAAATTTTTAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCa                                                                           <  1:3934144/150‑1 (MQ=255)
                                                                             ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGcc                                                                        >  1:2656091/1‑150 (MQ=255)
                                                                              gATCGCGTATAAAAAGCTTATTACCTGCTGACTCTGGCGCTGGGCGTGTTCCCGGCCTATCCCGCGGTTGCTGGTTTGCTGGCTGTTTATGCTCGCGACTAAATTTTTTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                       <  1:205263/150‑1 (MQ=255)
                                                                              gATCGCGTATAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGTGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATGATTAATTCTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCa                                                                       <  2:1358777/150‑1 (MQ=255)
                                                                               aTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAg                                                                      >  1:2703304/1‑150 (MQ=255)
                                                                                tCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGa                                                                     >  1:1062320/1‑150 (MQ=255)
                                                                                 cgcgTATAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGATTTCCAGGCCTATACCGCCTTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                                                    <  2:2022665/150‑1 (MQ=255)
                                                                                      aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGc                                                               >  2:4194765/1‑150 (MQ=255)
                                                                                       gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCa                                                              >  1:884152/1‑150 (MQ=255)
                                                                                        aaaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                             >  2:4319757/1‑150 (MQ=255)
                                                                                         aaaaGTTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCTATACCGCGGTTGCTGTTTTGCTGGCTGATTATGCTCGCTACTAAAATTTTTAATTCGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                            <  1:3708241/150‑1 (MQ=255)
                                                                                         aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGATGCTTATGCTCGCTACTAAAATGATTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                            <  2:4739942/150‑1 (MQ=255)
                                                                                         aaaaGCTTATGACCTGCTGACTCGGGCGCTGGGCTTGTTCCTGGCCGATACCGCCGTTGGTGTTTTGAAGGCTGCTTATGCTCGCGACTAAAATGATGACTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                            <  1:4743489/150‑1 (MQ=255)
                                                                                          aaaGTTTATTACCTGCTGAATCTGGCGCTGGGCTAGTTCCTGGCCTATACCGCCTTTGTTGTTTTGCTGGGTGATTTTGCTGGCGACTAAATTTTTTATTTCTTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                           <  1:2638024/150‑1 (MQ=255)
                                                                                          aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                           >  1:931728/1‑150 (MQ=255)
                                                                                           aaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGGGGTGTGCAGGCTTCTTTTGATCGCGCCGAAAAAGAATAATTCGTGAAACCGACCGGTATCCGTTCTGAAGGTAAGCCACCCGCGGCAATgg                                                          >  1:1856374/1‑150 (MQ=255)
                                                                                            aGCTTATGACCTGCTGACTCGGGCGCAGGGCGAGTTCCAGGCCGATCCCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                         <  2:4431613/150‑1 (MQ=255)
                                                                                              cTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                       >  1:2476266/1‑150 (MQ=255)
                                                                                               ttATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGTAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                      >  1:4527292/1‑150 (MQ=255)
                                                                                                tATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGAc                                                     <  1:3013835/150‑1 (MQ=255)
                                                                                                 aTGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg                                                    <  1:3023640/150‑1 (MQ=255)
                                                                                                  tGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACgg                                                   >  2:5316498/1‑150 (MQ=255)
                                                                                                   gACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGTATTAGTGAAAGCGACCGGTATACGTGCTGAAGGTCAGCCAGACGCGGCAAAGGAAGACGGc                                                  >  1:2742501/1‑150 (MQ=255)
                                                                                                   gACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGc                                                  >  2:1474563/1‑150 (MQ=255)
                                                                                                      gtgctgACTCGGGCGCTGGGGGTGTCCCCGGCCGATCCGGCTGTTGGTGTTTGGCTGGCTTCTTTTGCTCGCTACTAAATTTATTATTTCTTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                               <  1:1433598/149‑1 (MQ=255)
                                                                                                            aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAATAGTGAAAGCGACCGGTATCCGTGCGGAAGGTCAGCCCGACGCGGCCATGGAAGACGGCGATGCACTg                                         >  1:19822/1‑150 (MQ=255)
                                                                                                                 ggCGCTGGGCGTGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTTTGCTCGCGACGAAAATGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                    <  2:3971255/150‑1 (MQ=255)
                                                                                                                 ggCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                    <  2:378198/150‑1 (MQ=255)
                                                                                                                      aGGGCGAGTTCAAGGCCGATACCGCCGTTTGTGGTTTGCAGGCAGATTATGCTAGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATACGTGCTGTAGGTCAGACAGACGCGGCAATTGAAGACGGAGATGGACTGATTTTAATGa                               >  2:3158182/1‑150 (MQ=255)
                                                                                                                           gtgTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATCCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTc                          <  1:3257991/148‑1 (MQ=255)
                                                                                                                           gAGTTCCCGCCCGATCCCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACCAAAATTATGAATCCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTc                          <  1:1699560/150‑1 (MQ=255)
                                                                                                                            aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                         <  2:384317/150‑1 (MQ=255)
                                                                                                                            aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                         <  1:3538044/150‑1 (MQ=255)
                                                                                                                            aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                         <  2:40738/150‑1 (MQ=255)
                                                                                                                            aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                         >  2:2619854/1‑150 (MQ=255)
                                                                                                                                  aGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAAATAGTGAAAGCGACCGTTATACCTGCTGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACTGATTTTAAAGAAATTACGTGCTg                   >  1:5028466/1‑150 (MQ=255)
                                                                                                                                  aGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                   >  2:5045527/1‑150 (MQ=255)
                                                                                                                                    gCCGATACCGCGGTTGGTGGTTTGCTGGCTGATTTTGCTCGCTACTAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAc                 <  1:4577579/150‑1 (MQ=255)
                                                                                                                                    gCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATGCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAc                 <  2:695022/150‑1 (MQ=255)
                                                                                                                                     cctaTACCGCCGTTGCTGTTTTGCAGGCTGATTATGCTCGCGACAAAAATTATGAATTCTTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAcc                <  1:1978698/147‑1 (MQ=255)
                                                                                                                                       gATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgc              >  2:1716235/1‑150 (MQ=255)
                                                                                                                                               cgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGa      >  2:4742709/1‑150 (MQ=255)
                                                                                                                                               cgTTGCGGGTTTGCCGGCGGCTTATGCTCGCGACTAAAATTATGAATCCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGa      <  1:545001/150‑1 (MQ=255)
                                                                                                                                                gTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa     >  1:4044654/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATC  >  CP009273/3779133‑3779429

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