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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,779,2810CG 6.5% 759.9 / 2.0 244A221G (GCT→GGT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (127/97);  new base G (4/15);  total (131/113)
Fisher's exact test for biased strand distribution p-value = 3.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

AAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCA  >  CP009273/3779132‑3779430
                                                                                                                                                     |                                                                                                                                                     
 aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                      <  2:1891093/150‑1 (MQ=255)
 aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                                                                                                                                                      >  1:2169526/1‑150 (MQ=255)
  aaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                     >  2:865366/1‑150 (MQ=255)
     ttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                                  >  1:974318/1‑150 (MQ=255)
     ttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                                  <  2:2427206/150‑1 (MQ=255)
      tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                                                                                                                                                 >  2:313854/1‑150 (MQ=255)
       gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                                <  1:1783448/150‑1 (MQ=255)
          gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAg                                                                                                                                             <  1:1487623/150‑1 (MQ=255)
             cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                          <  1:920429/150‑1 (MQ=255)
             cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                          <  1:2141928/150‑1 (MQ=255)
              tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                         <  2:4433987/150‑1 (MQ=255)
               gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                        <  2:3548551/150‑1 (MQ=255)
               gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                        >  2:257675/1‑150 (MQ=255)
                ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct                                                                                                                                       <  1:2579801/150‑1 (MQ=255)
                 gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCtt                                                                                                                                      >  2:4509454/1‑150 (MQ=255)
                 gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCtt                                                                                                                                      <  2:552185/150‑1 (MQ=255)
                   aaaGGCCGCGTTGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                    <  1:171561/150‑1 (MQ=255)
                   aaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                    >  2:163816/1‑150 (MQ=255)
                    aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATg                                                                                                                                   >  2:4756746/1‑150 (MQ=255)
                        ccGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcg                                                                                                                               <  1:3553687/150‑1 (MQ=255)
                          gcgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcgcg                                                                                                                             <  1:4716651/150‑1 (MQ=255)
                            gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgac                                                                                                                           <  2:1442282/150‑1 (MQ=255)
                            gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgac                                                                                                                           <  1:3680454/150‑1 (MQ=255)
                                cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaa                                                                                                                       >  2:2811957/1‑150 (MQ=255)
                                 gTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaa                                                                                                                      <  2:4658939/150‑1 (MQ=255)
                                  tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAat                                                                                                                     <  2:3857012/150‑1 (MQ=255)
                                   ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACCAAAatg                                                                                                                    >  1:51146/1‑150 (MQ=255)
                                   ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatg                                                                                                                    <  2:2097767/150‑1 (MQ=255)
                                   ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGACGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaaag                                                                                                                    >  2:3451013/1‑148 (MQ=255)
                                    catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                                   >  2:1805046/1‑150 (MQ=255)
                                    catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                                   <  1:3962363/150‑1 (MQ=255)
                                    catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                                   >  1:240947/1‑150 (MQ=255)
                                      tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatgatg                                                                                                                 >  1:2733717/1‑150 (MQ=255)
                                        atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGaa                                                                                                               >  1:298469/1‑150 (MQ=255)
                                                  actacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAg                                                                                                     >  2:2060496/1‑150 (MQ=255)
                                                      cgcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAc                                                                                                 >  2:4191295/1‑150 (MQ=255)
                                                      cgcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAc                                                                                                 <  1:3308769/150‑1 (MQ=255)
                                                      cgcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAc                                                                                                 >  2:2473111/1‑150 (MQ=255)
                                                      cgcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAc                                                                                                 >  2:722396/1‑150 (MQ=255)
                                                       gcgTTGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGAcc                                                                                                <  2:2234644/150‑1 (MQ=255)
                                                         gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGt                                                                                              >  2:2323743/1‑150 (MQ=255)
                                                          aTGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGtt                                                                                             >  2:4236024/1‑150 (MQ=255)
                                                           tGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTa                                                                                            >  2:304543/1‑150 (MQ=255)
                                                            ggACCGCGTTAACGTTTGGTTTCGGGTATAAAAAGTTTATGACCTGGTGACTCTGGCGCTGGGCTAGTTCCTGGCCTATACGGCCGTTGGTGTTTTGCTGGGTGCTTTTGCTCGCGAGTAAAATGTTTATTTCGTGAAAGCGACCGTTAt                                                                                           <  2:4561003/150‑1 (MQ=255)
                                                                 gcgATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTg                                                                                      >  1:642709/1‑150 (MQ=255)
                                                                   tatccccGTTGGGATCGCGTATAAAATGCTTATGACGTGCTGACTGTGGCGCTGGGCGAGTTCCAGGCCGTTACCGCCGTTGCTGGTTTGCAGGCTGCTTTTGCTCGCGACTAAATTGTTTATTTCGTTAAAGCGACCGTTATCCGTGCt                                                                                    <  2:965106/145‑1 (MQ=255)
                                                                     tAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGa                                                                                  >  2:1372688/1‑150 (MQ=255)
                                                                       aCCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAg                                                                                >  1:3904555/1‑150 (MQ=255)
                                                                        ccGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAgg                                                                               >  2:3711857/1‑150 (MQ=255)
                                                                               gATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCa                                                                        >  2:1096342/1‑150 (MQ=255)
                                                                                aTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAg                                                                       >  1:4540269/1‑150 (MQ=255)
                                                                                aTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAAATCGTGAAAGCGACCGGTATCCGTGCAGAAGGTCAGCCAg                                                                       >  1:2839399/1‑150 (MQ=255)
                                                                                aTCGCGGAGAAAAAGAATATGACCTGCTGACTCTGTCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTTCTTATGCTAGCGACCAAAATGAAGAAATAGTGAAAGCTACCGTTATACGTGCAGAAAGTCAGcccg                                                                       >  2:3634569/1‑148 (MQ=255)
                                                                                   gcgTATAAAAAGCTTTTGACCTGCTGACTCGGGCGCTGGGCGAGTCCCAGGCCGATACCGCCGTTGCTGTTTTGCTGGCTGATTATGCTCGCTACGAAAATTATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcg                                                                    <  1:1857447/150‑1 (MQ=255)
                                                                                    cgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTAGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGAcgc                                                                   >  1:3645828/1‑150 (MQ=255)
                                                                                    cgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgc                                                                   >  2:2841912/1‑150 (MQ=255)
                                                                                     gTAGAAACAGATTATGACGTGCTGACTCTGGGGCTGGGCGAGTTCCAGGCCGATACCGCGGTTGCTGGTTTGCAGGCTGGTTATGCTCGCGACTAAAATGATGAATTGTTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                  <  1:4345073/150‑1 (MQ=255)
                                                                                     gTAGAAAAAGCTTAGGACCTGCTGACTCTGGCGCTGGGCGTGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACGAAAATGTTTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg                                                                  <  1:4753310/150‑1 (MQ=255)
                                                                                      tAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                                 >  1:3801033/1‑150 (MQ=255)
                                                                                       aGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGc                                                                >  2:4310096/1‑150 (MQ=255)
                                                                                         aaaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGTTGAATTCGTGTAAGCGACCGTTATACGTGCGGAAGGTCAGCCAGACGCGGCaa                                                              >  1:891717/1‑150 (MQ=255)
                                                                                          aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGTGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                             <  1:74078/150‑1 (MQ=255)
                                                                                          aaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAt                                                             >  2:2595160/1‑150 (MQ=255)
                                                                                           aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGACTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                            >  2:3661130/1‑150 (MQ=255)
                                                                                            aaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCTATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATgg                                                           <  2:1708306/150‑1 (MQ=255)
                                                                                             tgCTTCTGACCTGCTGACTCTGGCGCTGGGCGAGTTCCCGGCCGATACCGCGGTTGCTGGTTTGCTGGGTGCTTATGCTCGCGACTAAATTTTTAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                         <  1:3438203/149‑1 (MQ=255)
                                                                                             aGCTTCTGACCTGGTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGTTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                          <  2:1817388/150‑1 (MQ=255)
                                                                                             aGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGa                                                          >  1:1585073/1‑150 (MQ=255)
                                                                                              gCTTTTGACCTGCTGACTGTGGCGCTGGGGGTGTTCCAGGCCTATACCGCGGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACTAAATTGTTTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                         <  2:75120/150‑1 (MQ=255)
                                                                                              gCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGTTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa                                                         <  2:1686790/150‑1 (MQ=255)
                                                                                                ttATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                       >  2:3974408/1‑150 (MQ=255)
                                                                                                ttATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                       >  2:2019061/1‑150 (MQ=255)
                                                                                                 tctGACCTGCTGACTCTGGCGCTGGGCGAGTTCCCGGCCGATCCCGCGGTTGGTGGTTTGCTGGGGGCTTTTGCTGGCGACGAAATTGATGATTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGAc                                                      <  2:2179300/148‑1 (MQ=255)
                                                                                                   tGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACgg                                                    >  2:1480600/1‑150 (MQ=255)
                                                                                                   tGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACgg                                                    >  1:4308905/1‑150 (MQ=255)
                                                                                                     aCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAAAGATGAAATCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAAAGGAAGACGGCg                                                  >  1:1199795/1‑150 (MQ=255)
                                                                                                        tgcttACTGGGGGGCTGGGCGTGTTCCCGGCCTATCCCGCGTTTGGTGTTTGGCGGGCTGATTTTGCTCGCGACGAATTTTATTAATTCTTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                                               <  1:1348803/150‑1 (MQ=255)
                                                                                                           tgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                            >  2:1342696/1‑150 (MQ=255)
                                                                                                            gACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCCGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCCGACCCGGCCAAGGAAGACGGCGATGCAca                                           >  1:3729513/1‑149 (MQ=255)
                                                                                                                 tGGCGCTGGGCGAGTTCCTGGCCGATACCGCGGTTGGTGTTTTGCTGGCTGATTATGCTCGCGACGAAATTGATTATTTCGTGACAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAttt                                      <  1:2129499/150‑1 (MQ=255)
                                                                                                                  ggCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt                                     >  2:801849/1‑150 (MQ=255)
                                                                                                                       tgggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGGTGCTTATGCTCGCGACTAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                                <  1:1124507/149‑1 (MQ=255)
                                                                                                                          gCGAGTTCCCGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTAGGCTCGCGCCGAAAATGATGAATTCGTGAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACt                             <  2:2480981/150‑1 (MQ=255)
                                                                                                                           cGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATTATGAAATAGTGAAAGCGACCCTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGGTGCACGGGTTTTCATGAACtt                            >  1:3019959/1‑150 (MQ=255)
                                                                                                                            gAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTc                           >  2:1869270/1‑150 (MQ=255)
                                                                                                                             cgTTCCTGGCCGATACCGCGGTGGGTGTTTTGCAGGCTGCTTTTGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                          <  2:3247764/149‑1 (MQ=255)
                                                                                                                             aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGCCGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCGGACGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                          <  2:2714557/150‑1 (MQ=255)
                                                                                                                             aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                          >  1:115487/1‑150 (MQ=255)
                                                                                                                               ttCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGt                        >  2:3981401/1‑150 (MQ=255)
                                                                                                                                  cAGGCAGATACCGCGGTTGCTGGTTTGCAGGCTGTTTATGCTCGCTAAGAAATTTATGATTTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCt                     <  2:3192845/150‑1 (MQ=255)
                                                                                                                                   aGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTg                    >  2:2649207/1‑150 (MQ=255)
                                                                                                                                      ccGATACCGCCGTTGGTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAcc                 <  2:4586173/150‑1 (MQ=255)
                                                                                                                                      ccGATACCGCCGTTGCGGGTTTGCAGGCGGCTTAGGCGCGCGACGAAAATGATGAATGCGTTAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAcc                 <  2:1229737/150‑1 (MQ=255)
                                                                                                                                       ctaaaCCGCCGTTGCTGGTTTGCAGGCTGATTATGCTCGCTACGAAAATTATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcg                <  2:869640/146‑1 (MQ=255)
                                                                                                                                         aTACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGTAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcg              >  2:2357735/1‑150 (MQ=255)
                                                                                                                                           aCCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATTATGAATTAGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCAATGGAATACGGCGATGCACTGATTTTAATGAACTTCCGTGCTGACcgcgcg            >  1:3101634/1‑150 (MQ=255)
                                                                                                                                           aCCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATTAAATCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACCCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACATACGTGCTGACcgcgcg            >  1:3044521/1‑150 (MQ=255)
                                                                                                                                                cggggCTGGTGTGCTGGCTGCTTAGGCGCGCGACGAAAATGATGAATGCGTGAACGCGACCGTTATCCGTGCTGACGGTCAGCCAGACGCGGCAAGGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGa       <  2:4428624/146‑1 (MQ=255)
                                                                                                                                                 gTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACTGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa      >  2:112817/1‑150 (MQ=255)
                                                                                                                                                 gTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa      >  2:3497640/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCA  >  CP009273/3779132‑3779430

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: