Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,281 C→G 7.0% 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 7.0% 780.5 / 4.5 252A221G (GCT→GGT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (120/110);  new base G (5/17);  total (125/127)
Fisher's exact test for biased strand distribution p-value = 1.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

AAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATC  >  CP009273/3779132‑3779429
                                                                                                                                                     |                                                                                                                                                    
aaaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGc                                                                                                                                                      >  1:993179/1‑150 (MQ=255)
 aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                     <  2:3747901/150‑1 (MQ=255)
 aaaTGTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                     <  2:3825006/150‑1 (MQ=255)
  aaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                    >  1:354102/1‑150 (MQ=255)
   aTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTgg                                                                                                                                                   >  2:3314039/1‑150 (MQ=255)
     ttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                                 >  1:4456477/1‑150 (MQ=255)
      tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                >  1:1067145/1‑150 (MQ=255)
      tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                >  2:5154405/1‑150 (MQ=255)
      tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                >  1:601959/1‑150 (MQ=255)
      tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                >  1:2177364/1‑150 (MQ=255)
       gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                               <  1:2272142/150‑1 (MQ=255)
        ccGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGc                                                                                                                                              >  2:612730/1‑150 (MQ=255)
         cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa                                                                                                                                             >  1:1569194/1‑150 (MQ=255)
         cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa                                                                                                                                             <  2:3578523/150‑1 (MQ=255)
          gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAg                                                                                                                                            <  2:1615603/150‑1 (MQ=255)
           cgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAgg                                                                                                                                           >  1:2999216/1‑150 (MQ=255)
             cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                         <  2:5234516/150‑1 (MQ=255)
             cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                         >  2:3373327/1‑150 (MQ=255)
              tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                        >  2:1929201/1‑150 (MQ=255)
              tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                        <  1:4569762/150‑1 (MQ=255)
                ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct                                                                                                                                      >  1:4669468/1‑150 (MQ=255)
                  cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                    <  1:4048061/150‑1 (MQ=255)
                    aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATg                                                                                                                                  <  2:4651903/150‑1 (MQ=255)
                     aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                                 <  2:460854/150‑1 (MQ=255)
                       gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTc                                                                                                                               <  2:3988788/150‑1 (MQ=255)
                          gcgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcgcg                                                                                                                            <  2:4332213/150‑1 (MQ=255)
                            gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTAGcgac                                                                                                                          >  1:5098898/1‑150 (MQ=255)
                              aGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacga                                                                                                                        >  1:771630/1‑150 (MQ=255)
                              aGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCGGCTTATGCTCGcgacga                                                                                                                        >  1:487985/1‑150 (MQ=255)
                               gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaa                                                                                                                       <  1:1146973/150‑1 (MQ=255)
                               gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaa                                                                                                                       <  1:1284689/150‑1 (MQ=255)
                               gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaa                                                                                                                       >  2:4808426/1‑150 (MQ=255)
                                cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaa                                                                                                                      <  2:2356675/150‑1 (MQ=255)
                                cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGaaa                                                                                                                      >  1:2766559/1‑150 (MQ=255)
                                   ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatg                                                                                                                   >  1:4495112/1‑150 (MQ=255)
                                     atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgat                                                                                                                 >  2:4896584/1‑150 (MQ=255)
                                      tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgatg                                                                                                                >  1:2911919/1‑150 (MQ=255)
                                       catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                               >  2:1329610/1‑150 (MQ=255)
                                       catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                               >  2:4121198/1‑150 (MQ=255)
                                       catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCGGCGTATTCGCGCGACGCAAAtgatga                                                                                                               >  1:1041209/1‑150 (MQ=255)
                                        atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAAGATGaa                                                                                                              >  1:4040734/1‑150 (MQ=255)
                                           ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTc                                                                                                           <  1:3936150/150‑1 (MQ=255)
                                             tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTAGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGt                                                                                                         >  2:5050959/1‑150 (MQ=255)
                                               gTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGa                                                                                                       >  2:335359/1‑150 (MQ=255)
                                                ttACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaa                                                                                                      <  1:4966145/150‑1 (MQ=255)
                                                ttACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaa                                                                                                      <  1:1920974/150‑1 (MQ=255)
                                                ttACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaa                                                                                                      >  2:4723849/1‑150 (MQ=255)
                                                ttACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaa                                                                                                      <  1:150520/150‑1 (MQ=255)
                                                 tactacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGtaa                                                                                                     >  1:3142087/1‑150 (MQ=255)
                                                  actacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                    <  2:2310973/150‑1 (MQ=255)
                                                    tacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCg                                                                                                  <  1:3670366/150‑1 (MQ=255)
                                                      cgcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGGTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAc                                                                                                >  1:3276202/1‑150 (MQ=255)
                                                       gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAcc                                                                                               >  2:3635299/1‑150 (MQ=255)
                                                              aCCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATcc                                                                                        >  2:1704599/1‑150 (MQ=255)
                                                                     tAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCCGGCCGATACCGCGGTTGCTGGTTTGCTGGCTGCTTAGGCTCGCGACGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGa                                                                                 <  2:2675253/150‑1 (MQ=255)
                                                                             gggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGc                                                                         >  1:2578595/1‑150 (MQ=255)
                                                                             gggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGc                                                                         >  1:2356584/1‑150 (MQ=255)
                                                                                 tCGCTTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGTGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGATTATGCTCGCGACGAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGa                                                                     <  1:2601082/150‑1 (MQ=255)
                                                                                  cgctTATAAAAAGCTTTGTACCTGCTTACTGGGGCGCTGGGCTAGTACCCGGCCTATACCGCCGTTGGTGTTTTGCAGGCTGCTTATGCTCGCGACGAAAATTTTTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                                                    <  1:1988361/150‑1 (MQ=255)
                                                                                  cgcgTAGAAAAAGCTTTTGACCTGCTGACTCGGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                                                    <  2:4272818/150‑1 (MQ=255)
                                                                                       aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCGGAAGGTCAGCCAGACGCGGc                                                               >  2:3509078/1‑150 (MQ=255)
                                                                                        gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCa                                                              <  2:378036/150‑1 (MQ=255)
                                                                                        gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCa                                                              <  2:2999216/150‑1 (MQ=255)
                                                                                          aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                            >  2:17525/1‑150 (MQ=255)
                                                                                           aaaGCTTATGACGTGCTGACGCGGGCGCCGGGCGAGTGCCAGGCCGATCCCGCGGTTGCTGGTTTGCAGGCTGCTTAGGCGCGCTACTAAATTGATTAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                           <  1:210217/150‑1 (MQ=255)
                                                                                              gCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                        <  2:5019167/150‑1 (MQ=255)
                                                                                              gCTTAGGACCTGCTGACTCTGGCGCAGGGCGAGTTCCTGGCCGATACCGCCGTGGCTGGTTTGCAGGCTGCTTAGGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                        <  2:4595892/150‑1 (MQ=255)
                                                                                              gCTTAGGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCTATACCGCCGTTGGTGGTTTGAAGGCTGATTCTGCTCGCTACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                        <  2:1281913/150‑1 (MQ=255)
                                                                                                ttATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                      >  1:3902827/1‑150 (MQ=255)
                                                                                                 ttggACGTGCGGAGTCTGGGGCTGGGCGTGTTCCCGGCCGATACCGCCGGTGGGGTTTTGCTGGCGGCTTATGCGCGCGACGAAATTGTTGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGAc                                                     <  1:2866839/147‑1 (MQ=255)
                                                                                                    gACCTGCTGACTCTGGCGCTGGGCGTTTTCCTGGCCTATACCGCCGTTGGTGGTTTGCTGGTTGGTTTGGCTCGCTACTAAAATGATGAATTCTTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGc                                                  <  1:2051748/150‑1 (MQ=255)
                                                                                                     aCCTGCTGACTCTGGCGCTGGGCGTGTTCCTGGCCTATACCGCCGTTGCTGTTTTGCTGGCTGCTTTTGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg                                                 <  1:3028684/150‑1 (MQ=255)
                                                                                                      ccTGCTGACTGGGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGa                                                <  2:4456477/150‑1 (MQ=255)
                                                                                                      ccTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGCTCATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGa                                                <  2:3407677/150‑1 (MQ=255)
                                                                                                       ctgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                               >  2:1662482/1‑150 (MQ=255)
                                                                                                       ctgctgACTCGGGCGCTGGGCGAGTCCCTGGCCGATCCCGCCGTTGCTGGTTTGCCGGCTGCTTATGCTCGCGACGAAAATGAGGAATACGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGAt                                               <  1:4808426/150‑1 (MQ=255)
                                                                                                        tgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGTAAGCGACCGTTATACGTGCTGTAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                                              >  1:4124559/1‑150 (MQ=255)
                                                                                                         gctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGc                                             >  1:2791671/1‑150 (MQ=255)
                                                                                                          ctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCa                                            >  1:4231840/1‑150 (MQ=255)
                                                                                                                  ggCGCAGGGCGAGTTCCGGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCCTTATCCGTGCTGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACTGAtttt                                    >  1:1163413/1‑150 (MQ=255)
                                                                                                                   gcgcTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTc                                   <  2:2257680/150‑1 (MQ=255)
                                                                                                                   gcgcTGGGCGAGTGCCTGGCCGAGCCCGCCGTTGGGGGTTTGCAGGCTGCTTAGGCGCGCGACGAAAATGATGAATTCGGGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTc                                   <  2:2255812/150‑1 (MQ=255)
                                                                                                                    cgcgGGGCGTGTCCCAGGCCTATCCCGGGGTGGGTGTTTTGCTGGGTGCTTAGGCGCGCGACTAAATTTATGACTTGGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                                  <  1:4411094/150‑1 (MQ=255)
                                                                                                                      cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCAtt                                >  2:4555878/1‑149 (MQ=255)
                                                                                                                      cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATg                                <  2:361021/150‑1 (MQ=255)
                                                                                                                       aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                               >  1:425150/1‑150 (MQ=255)
                                                                                                                        gggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                              >  2:1099241/1‑150 (MQ=255)
                                                                                                                        gggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                              >  2:1099736/1‑150 (MQ=255)
                                                                                                                           ggtgTTCCAGGCCTATACCGCGGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACTAAAATGATTATTTCTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACtt                           <  1:1183696/147‑1 (MQ=255)
                                                                                                                             aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                         >  1:3486027/1‑150 (MQ=255)
                                                                                                                              gTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAGGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg                        <  2:2462178/150‑1 (MQ=255)
                                                                                                                                tCCAGGCCGATACCGCCGTTGGGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATCCGTGCTGAATGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTg                      >  1:351166/1‑150 (MQ=255)
                                                                                                                                 ccAGGCCGATACCGCCGTTGCTGGTTTGCAGGCGGCTTATGCTCGCGACGAAAATGATTAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGc                     <  2:2800263/150‑1 (MQ=255)
                                                                                                                                   aGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                   >  2:2981539/1‑150 (MQ=255)
                                                                                                                                     gCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAc                 <  2:2043117/150‑1 (MQ=255)
                                                                                                                                       cGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAAAGATGAATTCGTGAAAGCGACCGTTATACCTGCAGAAGGTCAGCCACACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcg               >  1:5059824/1‑150 (MQ=255)
                                                                                                                                           aCCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcgcg           <  2:1422157/150‑1 (MQ=255)
                                                                                                                                              gccgTTGCTGGTTTGAAGGTTGCTTTTGCTCGCGACTAAAATTATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt        <  2:2177364/150‑1 (MQ=255)
                                                                                                                                               cggTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAACTGTTGATTTCGGGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTg       <  1:796803/148‑1 (MQ=255)
                                                                                                                                                   tGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAAt   >  1:5236713/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATC  >  CP009273/3779132‑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: