Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,409,935 G→T 100% A570A (GCG→GCT pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,409,9350GT98.6% 471.8 / ‑1.7 147A570A (GCG→GCTptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (61/84);  total (63/84)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 147 total aligned reads displayed.

GTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCGATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTG  >  CP009273/2409790‑2410077
                                                                                                                                                 |                                                                                                                                              
 tccgCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTAtt                                                                                                                                             >  2:2172865/1‑148 (MQ=255)
    gccgcTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAg                                                                                                                                          <  2:1745271/148‑1 (MQ=255)
     ccgcTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGt                                                                                                                                         >  1:2297933/1‑148 (MQ=255)
          gCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCt                                                                                                                                    >  1:1146195/1‑148 (MQ=255)
           cAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTg                                                                                                                                   >  1:2070008/1‑148 (MQ=255)
              cTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGAtt                                                                                                                                >  1:2541454/1‑148 (MQ=255)
                  tCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGc                                                                                                                            >  1:1264395/1‑148 (MQ=255)
                    aaaaCTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTg                                                                                                                          >  1:2645056/1‑148 (MQ=255)
                      aaCTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCg                                                                                                                        >  2:647555/1‑148 (MQ=255)
                      aaCTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCg                                                                                                                        >  2:2408918/1‑148 (MQ=255)
                         tGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGcc                                                                                                                     <  2:2377043/148‑1 (MQ=255)
                           cACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCtt                                                                                                                   <  1:3418952/148‑1 (MQ=255)
                            aCCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTc                                                                                                                  >  2:3477533/1‑148 (MQ=255)
                               gggCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGt                                                                                                               <  2:2578850/148‑1 (MQ=255)
                               gggCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGt                                                                                                               <  1:764454/148‑1 (MQ=255)
                                 gCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTAt                                                                                                             >  2:1393875/1‑148 (MQ=255)
                                    gCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGa                                                                                                          >  2:207440/1‑148 (MQ=255)
                                    gCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGa                                                                                                          <  1:2126788/148‑1 (MQ=255)
                                    gCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGa                                                                                                          <  1:2126797/148‑1 (MQ=255)
                                    gCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGa                                                                                                          >  2:208185/1‑148 (MQ=255)
                                     cTCCCTGGGATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGaa                                                                                                         <  1:3144501/148‑1 (MQ=255)
                                       cccTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAAcc                                                                                                       >  1:2713412/1‑148 (MQ=255)
                                          tGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACcgcg                                                                                                    >  2:2722790/1‑148 (MQ=255)
                                               tCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTg                                                                                               >  1:2810714/1‑148 (MQ=255)
                                               tCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTg                                                                                               <  2:1637966/148‑1 (MQ=255)
                                               tCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTg                                                                                               <  2:1636053/148‑1 (MQ=255)
                                                cTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGc                                                                                              >  2:1940890/1‑148 (MQ=255)
                                                cTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGc                                                                                              >  1:958501/1‑148 (MQ=255)
                                                cTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCGATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGc                                                                                              >  2:1036745/1‑148 (MQ=255)
                                                  tCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTa                                                                                            <  1:1465868/148‑1 (MQ=255)
                                                  tCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTa                                                                                            <  1:2395521/148‑1 (MQ=255)
                                                  tCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTa                                                                                            <  2:843996/148‑1 (MQ=255)
                                                  tCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCCGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTa                                                                                            <  1:2395576/148‑1 (MQ=255)
                                                    cGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATg                                                                                          <  2:927980/148‑1 (MQ=255)
                                                    cGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATg                                                                                          <  1:2314348/148‑1 (MQ=255)
                                                    cGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATg                                                                                          <  1:3260586/148‑1 (MQ=255)
                                                      tgtTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGct                                                                                        >  2:3324160/1‑148 (MQ=255)
                                                       gtTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGctc                                                                                       >  2:1771229/1‑148 (MQ=255)
                                                        ttcttcATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGctct                                                                                      <  1:267206/148‑1 (MQ=255)
                                                        ttcttcATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGctct                                                                                      <  1:2801882/148‑1 (MQ=255)
                                                        ttcttcATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGctct                                                                                      <  2:1000241/148‑1 (MQ=255)
                                                          cttcATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTcc                                                                                    <  1:2660454/148‑1 (MQ=255)
                                                          cttcATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTcc                                                                                    <  1:342258/148‑1 (MQ=255)
                                                          cttcATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCAGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTcc                                                                                    <  1:723619/148‑1 (MQ=255)
                                                               tgctgctgcCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTc                                                                               >  1:1002936/1‑148 (MQ=255)
                                                               tgctgctgcCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTc                                                                               <  1:2949017/148‑1 (MQ=255)
                                                                gctgctgcCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCGATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTcc                                                                              >  2:1919004/1‑148 (MQ=255)
                                                                    ctgcCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCg                                                                          >  1:2627936/1‑148 (MQ=255)
                                                                        cGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTAc                                                                      <  1:1703645/148‑1 (MQ=255)
                                                                        cGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTAc                                                                      <  2:1253084/148‑1 (MQ=255)
                                                                           aaCAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTc                                                                   >  2:1819213/1‑148 (MQ=255)
                                                                                gtttacgtttacGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTg                                                              >  2:2797049/1‑148 (MQ=255)
                                                                                gtttacgtttacGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTg                                                              >  2:2797434/1‑148 (MQ=255)
                                                                                gtttacgtttacGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTg                                                              <  2:1411123/148‑1 (MQ=255)
                                                                                  ttacgtttacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCa                                                            >  1:1938185/1‑148 (MQ=255)
                                                                                  ttacgtttacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCa                                                            <  1:2093201/148‑1 (MQ=255)
                                                                                  ttacgtttacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCa                                                            <  2:1901846/148‑1 (MQ=255)
                                                                                  ttacgtttacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCa                                                            <  1:2172865/148‑1 (MQ=255)
                                                                                    acgtttacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAgg                                                          >  1:913848/1‑148 (MQ=255)
                                                                                       tttacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAg                                                       <  1:729207/148‑1 (MQ=255)
                                                                                         tacgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCg                                                     <  2:250839/148‑1 (MQ=255)
                                                                                          acgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGa                                                    <  2:1726509/148‑1 (MQ=255)
                                                                                          acgGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGa                                                    <  1:919703/148‑1 (MQ=255)
                                                                                               gACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAg                                               <  2:315200/148‑1 (MQ=255)
                                                                                                aCTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGa                                              >  1:782110/1‑148 (MQ=255)
                                                                                                        aTCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTc                                      >  2:25420/1‑148 (MQ=255)
                                                                                                         tCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTcg                                     <  1:3146755/148‑1 (MQ=255)
                                                                                                         tCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTcg                                     >  1:49349/1‑148 (MQ=255)
                                                                                                            aCCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGa                                  <  1:1995431/148‑1 (MQ=255)
                                                                                                               cGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGc                               >  2:2661024/1‑148 (MQ=255)
                                                                                                               cGGATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGc                               >  2:3392500/1‑148 (MQ=255)
                                                                                                                ggATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCa                              <  1:2412305/148‑1 (MQ=255)
                                                                                                                ggATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCa                              <  1:2412237/148‑1 (MQ=255)
                                                                                                                 gATCCGACCGCTGAACAGCTGGCAGAACTCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCaa                             <  2:1197882/148‑1 (MQ=255)
                                                                                                                 gATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCTGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCGGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCaa                             <  2:2707600/148‑1 (MQ=255)
                                                                                                                 gATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCaa                             <  2:931296/148‑1 (MQ=255)
                                                                                                                 gATCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCaa                             <  1:528283/148‑1 (MQ=255)
                                                                                                                   tCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACt                           >  1:257438/1‑148 (MQ=255)
                                                                                                                   tCCGACCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACt                           <  1:1388463/148‑1 (MQ=255)
                                                                                                                       aCCGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTc                       >  2:1288981/1‑148 (MQ=255)
                                                                                                                        ccGCTGAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCt                      <  2:392165/148‑1 (MQ=255)
                                                                                                                             gAACAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACGCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGgcgc                 <  1:2089049/148‑1 (MQ=255)
                                                                                                                              aaCAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCa                <  2:2628819/148‑1 (MQ=255)
                                                                                                                                cAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAgg              >  1:3002112/1‑148 (MQ=255)
                                                                                                                                cAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAgg              >  1:1040877/1‑148 (MQ=255)
                                                                                                                                cAGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAgg              >  1:1644524/1‑148 (MQ=255)
                                                                                                                                 aGCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGa             >  1:2784050/1‑148 (MQ=255)
                                                                                                                                  gCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGaa            >  2:2062815/1‑148 (MQ=255)
                                                                                                                                  gCTGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGaa            >  2:2062869/1‑148 (MQ=255)
                                                                                                                                    tGGCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGaaaa          <  1:792319/148‑1 (MQ=255)
                                                                                                                                     ggCAGAAATCGCTATTCCGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGaaaaa         <  2:324355/148‑1 (MQ=255)
                                                                                                                                     ggCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGaaaaa         <  1:2761297/148‑1 (MQ=255)
                                                                                                                                     ggCAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGaaaaa         <  2:324422/148‑1 (MQ=255)
                                                                                                                                       cAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACg       <  1:1545746/148‑1 (MQ=255)
                                                                                                                                       cAGAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACg       >  1:3304877/1‑148 (MQ=255)
                                                                                                                                         gAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTCCTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTc     <  1:207440/148‑1 (MQ=255)
                                                                                                                                         gAAATCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTc     <  2:2295006/148‑1 (MQ=255)
                                                                                                                                         gAAATCGCTATTCAGTCCGCTGAGTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTc     <  2:2294264/148‑1 (MQ=255)
                                                                                                                                          aaaTCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTcc    >  2:2955526/1‑148 (MQ=255)
                                                                                                                                            aTCGCTATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTcctg  <  2:2151240/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCGGATCCGACCGCTGAACAGCTGGCAGAAATCGCGATTCAGTCCGCTGATTCCGCTGCGGCCTTCGGTATCGAACCGCGCGTTGCTATGCTCTCCTACTCCACCGGTACTTCTGGTGCAGGTAGCGACGTAGAAAAAGTTCGCGAAGCAACTCGTCTGGCGCAGGAAAAACGTCCTG  >  CP009273/2409790‑2410077

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