Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,779,281 C→G 7.2% 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.2% 504.5 / 2.2 177A221G (GCT→GGT) gpmMphosphoglycero mutase III, cofactor‑independent
Reads supporting (aligned to +/- strand):  ref base C (79/82);  new base G (6/9);  total (85/92)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

AATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAA  >  CP009273/3779134‑3779426
                                                                                                                                                   |                                                                                                                                                 
aaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGc                                                                                                                                                   <  2:3189886/148‑1 (MQ=255)
aaTGTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGc                                                                                                                                                   <  2:3189883/148‑1 (MQ=255)
  tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                 <  1:2956155/148‑1 (MQ=255)
  tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                                                                                                                                                 <  2:127056/148‑1 (MQ=255)
     gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                                                                                                                                              >  1:1758320/1‑148 (MQ=255)
       cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                            <  1:340731/148‑1 (MQ=255)
       cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                            <  1:340752/148‑1 (MQ=255)
       cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                            <  2:1729240/148‑1 (MQ=255)
       cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg                                                                                                                                            <  1:705833/148‑1 (MQ=255)
           cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAgg                                                                                                                                        <  1:2260384/148‑1 (MQ=255)
            tGGGCAAAGGCCGCGTAGCGTCCATCATTTGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc                                                                                                                                       <  2:2512048/148‑1 (MQ=255)
             gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct                                                                                                                                      >  1:267206/1‑148 (MQ=255)
              ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                     >  1:3150299/1‑148 (MQ=255)
              ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg                                                                                                                                     <  1:45789/148‑1 (MQ=255)
               gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc                                                                                                                                    <  1:3003257/148‑1 (MQ=255)
                  aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa                                                                                                                                 <  2:2616260/148‑1 (MQ=255)
                   aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt                                                                                                                                <  2:859673/148‑1 (MQ=255)
                    ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATg                                                                                                                               >  2:2980031/1‑148 (MQ=255)
                     gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                              >  2:656483/1‑148 (MQ=255)
                     gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                              >  2:2054398/1‑148 (MQ=255)
                     gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                              >  2:2054430/1‑148 (MQ=255)
                     gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGc                                                                                                                              >  2:2219655/1‑148 (MQ=255)
                      ccGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCt                                                                                                                             <  1:999500/148‑1 (MQ=255)
                        gcgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcg                                                                                                                           <  1:573575/148‑1 (MQ=255)
                          gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTcgcg                                                                                                                         <  2:2249020/148‑1 (MQ=255)
                           tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcga                                                                                                                        <  2:252341/148‑1 (MQ=255)
                           tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcga                                                                                                                        >  2:2909608/1‑148 (MQ=255)
                             gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacg                                                                                                                      >  1:770800/1‑148 (MQ=255)
                             gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacg                                                                                                                      >  1:2166867/1‑148 (MQ=255)
                              cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGcgacga                                                                                                                     <  2:1998040/148‑1 (MQ=255)
                              cGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTGATGCTCGcgacga                                                                                                                     >  2:2639419/1‑148 (MQ=255)
                                tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaa                                                                                                                   <  1:524416/148‑1 (MQ=255)
                                tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaa                                                                                                                   <  2:704876/148‑1 (MQ=255)
                                 ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGaaaa                                                                                                                  <  2:200539/148‑1 (MQ=255)
                                  catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGAACTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAat                                                                                                                 <  1:1064996/148‑1 (MQ=255)
                                    tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGGAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                               >  1:455420/1‑148 (MQ=255)
                                    tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                               <  1:2541230/148‑1 (MQ=255)
                                    tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAatga                                                                                                               <  1:2259618/148‑1 (MQ=255)
                                       ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                            <  2:2983700/148‑1 (MQ=255)
                                       ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAAtgatga                                                                                                            >  2:1165807/1‑148 (MQ=255)
                                        tGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGTGGTTTGCAGGCTGCTTATGCTCGCGCCGAAAATGATGaa                                                                                                           >  1:1535663/1‑148 (MQ=255)
                                        tGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGaa                                                                                                           <  2:765680/148‑1 (MQ=255)
                                         ggTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCGTATTCTCGCCCCGAAAATGATGAAt                                                                                                          >  1:897536/1‑148 (MQ=255)
                                           tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTAGCGACGAAAAATATGAATTc                                                                                                        >  1:2542957/1‑148 (MQ=255)
                                           tCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGCTGGTTTTAACGCTGCTTATTCTCGCGACGAAAAAGATTaaata                                                                                                        >  1:3085273/1‑145 (MQ=255)
                                               tactacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCGTATGCTCGCGACGAAAATGATGAATTCGTGa                                                                                                    >  1:3270643/1‑148 (MQ=255)
                                                 ctacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaaa                                                                                                  <  2:1311500/148‑1 (MQ=255)
                                                 ctacGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGaaa                                                                                                  <  1:2011119/148‑1 (MQ=255)
                                                   acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGc                                                                                                <  1:2678033/148‑1 (MQ=255)
                                                   acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGc                                                                                                >  2:724115/1‑148 (MQ=255)
                                                   acGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGAGGAATTCGTGAAAGc                                                                                                >  2:2847985/1‑148 (MQ=255)
                                                         tGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGt                                                                                          >  1:2285906/1‑148 (MQ=255)
                                                           gACCGCGATACCCGTTGGGATCGCGTAGAAAAAGCTTTTGCCCTGCTGACTGGGGCGCTGGGCGCGTTCCAGGCCTATACCGCCGTTGCTGGTTTGCTGGCTGCTTATGCTCGCGACTAAAATGATGAATGCGTGAAAGCGACCGTTa                                                                                        <  2:2255105/148‑1 (MQ=255)
                                                                     aCCGTTGGGATCGCGTATAAAAAGCTTATGACCTGCTGACTCTGGCGCTGGGCGAGTTCCAGGCCGATCCCGCCGTTGCTGTTTTGCAGGCTGCTTATGCTCGCGACGAAATTGATGATTTCGTGAACGCGACCGTTATCCGTGCTGa                                                                              <  2:3072141/148‑1 (MQ=255)
                                                                     aCCGTGGGGAGCGCGTATAAAAAGCTTATCACCTGCTGACTCTGGCGCTGGGCGAGTTCCCGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGATTATGCTCGCGACGAAATTGATGAATTCGTCAAAGCGACCGTTATCCGTGCTGa                                                                              <  2:2793403/148‑1 (MQ=255)
                                                                       cGTTGGGATCGCTTAGAAAAAGCTTATCACCTGCGGACTCTGGCGCAGGGCGTGTTCCAGCCCGATACCGCCGTTGCTGGTTTGCAGGCTGATTATGCTCGCGACGAAAATGATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAg                                                                            <  2:2328963/148‑1 (MQ=255)
                                                                         ttGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGt                                                                          >  2:2620173/1‑148 (MQ=255)
                                                                            ggATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAg                                                                       >  2:101757/1‑148 (MQ=255)
                                                                                cgcgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCACGCGGCTTATGCTCGCGACGAAAAAGATGAATTCGTGAACTCGACCCTTATCACTGCTGAAGGTCATCCAg                                                                   >  1:113790/1‑148 (MQ=255)
                                                                                  cgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAc                                                                 >  2:3247208/1‑148 (MQ=255)
                                                                                      gAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCgg                                                             >  2:2364077/1‑148 (MQ=255)
                                                                                       aaaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTTGTGGGTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCACACGCGGc                                                            >  1:1511877/1‑148 (MQ=255)
                                                                                       aaaaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGGTATCCGTGCTGAAGGTCAGCCAGACGCGGc                                                            >  1:3231268/1‑148 (MQ=255)
                                                                                         aaaGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTGATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCaa                                                          >  1:1599372/1‑148 (MQ=255)
                                                                                           aGCTTCTGACCTGCTGACTCGGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAATTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                        <  2:1624989/148‑1 (MQ=255)
                                                                                           aGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATg                                                        <  2:2432739/148‑1 (MQ=255)
                                                                                            gCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATgg                                                       >  2:554845/1‑148 (MQ=255)
                                                                                               tATGCCGTGCTGACTGTGGCGCTGGGCGTGTCCCCGGCGGATACGGCGGTTGGTGTTTTCCCGGCTGATTAGGCACGCGACTAAATTTTTTAATTCGTCAACGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg                                                    <  1:1165807/148‑1 (MQ=255)
                                                                                                aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                   <  1:2182816/148‑1 (MQ=255)
                                                                                                aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGa                                                   <  2:1373526/148‑1 (MQ=255)
                                                                                                   aCCTGCTGACTCTGGCGCCGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACgg                                                <  2:943288/148‑1 (MQ=255)
                                                                                                     ctgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg                                              >  1:1042564/1‑148 (MQ=255)
                                                                                                      tgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGGGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCCGACGCGGCAATGGAAGACGGCGa                                             >  1:1255474/1‑148 (MQ=255)
                                                                                                      tgctgACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTTCACGCTGCTTATTCTAGCGACGAAAATGATTAATTCGTGAAATCGACCCTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAAAGGAAGACCGCGa                                             >  1:215224/1‑148 (MQ=255)
                                                                                                        ctgACTCGGGCGCCGGGCTATTCCCCGGCCTATACCGCGTTGGGTGGTTTGCTGGGTGCTTATGCTCGCTACTAAATTTTTGATTTCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg                                           <  1:1757129/148‑1 (MQ=255)
                                                                                                          gACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTGGGGGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGTAAGCGACCGGTATCCGTGCGGTAGGTCAGCCAGACGCCGCCATGGAAGACGGGGATGCa                                         >  1:288540/1‑148 (MQ=255)
                                                                                                           aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                        >  2:1574025/1‑148 (MQ=255)
                                                                                                           aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                        >  2:1753783/1‑148 (MQ=255)
                                                                                                           aCTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc                                        >  1:382180/1‑148 (MQ=255)
                                                                                                              ctGGCGCTGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATTATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGa                                     <  2:1978049/148‑1 (MQ=255)
                                                                                                                    cAGGGCGAGTTCCAGGCCGATACCGCCTTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACTAAAATGATGAATTGGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                               <  2:2825426/148‑1 (MQ=255)
                                                                                                                    cAGGGCGAGTTCCAGGCCGATACCGCCGTTGCGGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAATCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCa                               >  2:965804/1‑148 (MQ=255)
                                                                                                                       ggCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGa                            >  2:2878376/1‑148 (MQ=255)
                                                                                                                        gCGTGTTCCAGGCCGATCCGGCCGTTGGTGGTTTGCTGGCTGGTTTTGCTCGCGACTAAAATTATGAATTCGTTAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                           <  2:3085273/148‑1 (MQ=255)
                                                                                                                        gCGTGTTCCAGGCCGATACCGCCGTTGTTGTTTTGCAGGCTGCTTAGGCTCGCGACTAAATTTATTAATTCGTCAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGaa                           <  1:145665/148‑1 (MQ=255)
                                                                                                                           aGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATTCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCGGAAGGTCAGCCAGACGCGGCCATGGAAGACGGCGATGCACTGATTTTCATGAACtt                        >  1:2669787/1‑148 (MQ=255)
                                                                                                                             ttCCAGGCCGATACCGCCTTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                      <  2:542277/148‑1 (MQ=255)
                                                                                                                             ttCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTcc                      <  2:820669/148‑1 (MQ=255)
                                                                                                                              tCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg                     >  2:2726914/1‑148 (MQ=255)
                                                                                                                              tCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg                     >  2:33058/1‑148 (MQ=255)
                                                                                                                              tCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTAATGAACTTCCg                     >  2:1572808/1‑148 (MQ=255)
                                                                                                                              gccgggCCGATCCGGCGGTTGGTGGTTGGCCGGCTGGTTATGCTCGCGCCGAAAATGATTAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg                     <  1:2897788/144‑1 (MQ=255)
                                                                                                                                   gCCGATACCGCCGTTGCTGGTTTTCAGGCTGCTTATTCTCGCGACGAAAAATATGAATTAGTGAAAGCCACCCTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTAATGAACTTCCGTGCTg                >  1:1283255/1‑148 (MQ=255)
                                                                                                                                   gCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAAATAGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACAGATTTTCATGAAATTCCGTGCTg                >  1:2929842/1‑148 (MQ=255)
                                                                                                                                       aTACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTAGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcg            >  2:1196855/1‑148 (MQ=255)
                                                                                                                                       aTACCGCCGTTGCGGGTTTTCAGGCTGCTTATGCTCGCGACGAAAATTATTAAATAGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcg            >  1:2438361/1‑148 (MQ=255)
                                                                                                                                         aCCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTGTCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACATCCGTGCTGACcgcg          >  2:359005/1‑148 (MQ=255)
                                                                                                                                           cgccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATACGTGCTGAAGGTCAGCCAGACGCGGCAAAGGAAGACGGCGATGCACTGATTTTCATTAACATCCGTGCTGACcgcgcg        >  1:718607/1‑148 (MQ=255)
                                                                                                                                                 tGCTGTTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa  <  2:422753/148‑1 (MQ=255)
                                                                                                                                                 ggCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGaa  <  2:2560971/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAA  >  CP009273/3779134‑3779426

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: