Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 108,039 C→A 88.3% intergenic (+119/+97) mutT → / ← yacG nucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP/DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273108,0390CA88.3% 512.7 / 64.5 212intergenic (+119/+97)mutT/yacGnucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP/DNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (7/18);  new base A (89/98);  total (96/116)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

AATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGC  >  CP009273/107896‑108182
                                                                                                                                               |                                                                                                                                               
aaTTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGAt                                                                                                                                             >  2:2893608/1‑148 (MQ=255)
  ttGCGAAGCTTACACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATaa                                                                                                                                           <  1:574069/148‑1 (MQ=255)
  ttGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATaa                                                                                                                                           <  1:2395380/148‑1 (MQ=255)
  ttGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATaa                                                                                                                                           >  2:2726615/1‑148 (MQ=255)
    gCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAgg                                                                                                                                         <  1:1170359/148‑1 (MQ=255)
      gAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCg                                                                                                                                       <  2:3091497/148‑1 (MQ=255)
         gCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTc                                                                                                                                    <  1:2260768/148‑1 (MQ=255)
         gCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTc                                                                                                                                    >  1:646923/1‑148 (MQ=255)
          cTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCa                                                                                                                                   >  1:2073274/1‑148 (MQ=255)
          cTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCa                                                                                                                                   >  1:685314/1‑148 (MQ=255)
           ttAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAc                                                                                                                                  >  2:2136000/1‑148 (MQ=255)
            tAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATATATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcg                                                                                                                                 <  1:186270/148‑1 (MQ=255)
            tAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcg                                                                                                                                 >  1:1449235/1‑148 (MQ=255)
            tAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcg                                                                                                                                 >  2:1110408/1‑148 (MQ=255)
            tAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcg                                                                                                                                 >  1:2303140/1‑148 (MQ=255)
             aaaCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgc                                                                                                                                >  1:1982111/1‑148 (MQ=255)
               aCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgccg                                                                                                                              >  1:2634239/1‑148 (MQ=255)
                 gTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCa                                                                                                                            <  2:2270369/148‑1 (MQ=255)
                  tCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCAt                                                                                                                           <  2:85606/148‑1 (MQ=255)
                  tCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCAt                                                                                                                           <  1:1002846/148‑1 (MQ=255)
                         gTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACAt                                                                                                                    >  1:2344333/1‑148 (MQ=255)
                         gTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACAt                                                                                                                    <  2:2806651/148‑1 (MQ=255)
                         gTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACAt                                                                                                                    <  1:1267955/148‑1 (MQ=255)
                           cAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACAttt                                                                                                                  >  2:1666131/1‑148 (MQ=255)
                            aGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTg                                                                                                                 >  2:1422042/1‑148 (MQ=255)
                              aTAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGca                                                                                                               >  2:377508/1‑148 (MQ=255)
                                aaGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGcaca                                                                                                             <  2:1590574/148‑1 (MQ=255)
                                 aGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACaa                                                                                                            >  2:1954998/1‑148 (MQ=255)
                                  ggCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATcc                                                                                                                         <  2:719938/134‑1 (MQ=255)
                                  ggCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAg                                                                                                           >  1:192515/1‑148 (MQ=255)
                                   gCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGa                                                                                                          <  2:2552299/148‑1 (MQ=255)
                                      ttttCGCAGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTTTCAGGCCTAAAGGGATTTCTAACTCATTGTTAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGc                                                                                                       <  1:3270818/148‑1 (MQ=255)
                                       tttCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGcc                                                                                                      >  1:794676/1‑148 (MQ=255)
                                         tCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTg                                                                                                    <  1:1996174/148‑1 (MQ=255)
                                          cgccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGa                                                                                                   <  1:1599032/148‑1 (MQ=255)
                                          cgccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGa                                                                                                   <  1:1040409/148‑1 (MQ=255)
                                              gcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCg                                                                                               <  2:1934625/148‑1 (MQ=255)
                                                  ccGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTTGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGATCAAGATGCCTGATGAGACGc                                                                                           >  2:3063954/1‑148 (MQ=39)
                                                    gACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTg                                                                                         <  1:487356/148‑1 (MQ=255)
                                                     aCATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGt                                                                                        >  1:269400/1‑148 (MQ=255)
                                                     aCATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGt                                                                                        <  1:982439/148‑1 (MQ=255)
                                                     aCATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGt                                                                                        >  1:270368/1‑148 (MQ=255)
                                                     aCATTCCCACACGATGCCTGATGCGACGCTGGAGCGTCTTAACATGCCTAAACGGATTTGTAACTCAATGATAAATTTGTTTTTGTATGTAGGATAAGGCGTTCACGCCGCATCCGACATTTTCACAATAAGCCTGATGCGACGCTGt                                                                                        >  1:1454732/1‑148 (MQ=255)
                                                           gCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGt                                                                                  <  2:2157598/148‑1 (MQ=255)
                                                              cacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCtt                                                                               <  2:500093/148‑1 (MQ=255)
                                                               acGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTa                                                                              <  2:3041631/148‑1 (MQ=255)
                                                                cGATGCGTGAGGCGACGCTGGCGCGTCTTATCAGGCCTAACCGGATTTCGAACTCATTGATAAATTTGATTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACACGATGCCTGATGCGACGCTGTCCGCGTCTTAt                                                                             <  1:2024235/148‑1 (MQ=37)
                                                                cGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTAt                                                                             <  1:1365116/148‑1 (MQ=255)
                                                                  aTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCa                                                                           <  2:1859528/148‑1 (MQ=255)
                                                                     ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGc                                                                        <  1:2893608/148‑1 (MQ=255)
                                                                     ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGc                                                                        >  2:2728197/1‑148 (MQ=255)
                                                                     ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGc                                                                        >  1:1990102/1‑148 (MQ=255)
                                                                         aTGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTAc                                                                    >  1:498158/1‑148 (MQ=255)
                                                                         aTGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTAc                                                                    >  1:2197871/1‑148 (MQ=255)
                                                                           gCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGt                                                                  <  1:2090072/148‑1 (MQ=255)
                                                                                   ggCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATc                                                          <  1:3252746/148‑1 (MQ=255)
                                                                                     cgcgTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAg                                                        <  1:2837896/148‑1 (MQ=255)
                                                                                      gcgTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGa                                                       <  2:2344313/148‑1 (MQ=255)
                                                                                          cTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                   >  2:903188/1‑148 (MQ=255)
                                                                                          cTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                   >  2:3312501/1‑148 (MQ=255)
                                                                                             aTCATGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGt                                                <  1:1706379/148‑1 (MQ=255)
                                                                                             aTCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGt                                                <  2:358904/148‑1 (MQ=255)
                                                                                               cAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCa                                              <  2:242195/148‑1 (MQ=255)
                                                                                                aGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCAc                                             <  1:262554/148‑1 (MQ=255)
                                                                                                  gCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTg                                           >  1:2502733/1‑148 (MQ=255)
                                                                                                    cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCt                                         <  2:50891/148‑1 (MQ=255)
                                                                                                    cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCt                                         <  2:875461/148‑1 (MQ=255)
                                                                                                      aaaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCAATGCttt                                       >  2:2779361/1‑148 (MQ=255)
                                                                                                        aGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTgg                                     >  1:1779649/1‑148 (MQ=255)
                                                                                                          ggATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGtt                                   <  1:1583981/148‑1 (MQ=255)
                                                                                                            aTTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttct                                 >  1:1530002/1‑148 (MQ=255)
                                                                                                            aTTTCTAACTCATTGATAAATTTGTTTTTGTAGGGAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttct                                 <  2:1225433/148‑1 (MQ=255)
                                                                                                            aTTTCTAACTCATTGATAAATTTGTGTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttct                                 <  2:1718161/148‑1 (MQ=255)
                                                                                                             tttCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttctt                                <  1:2466324/148‑1 (MQ=255)
                                                                                                              ttCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                               >  2:106188/1‑148 (MQ=255)
                                                                                                              ttCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCACATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGACACTGCTTTGGttcttc                               >  1:526392/1‑148 (MQ=255)
                                                                                                               tCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCg                              <  1:1456991/148‑1 (MQ=255)
                                                                                                               tCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCg                              <  1:1272664/148‑1 (MQ=255)
                                                                                                                 tAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCt                            <  1:404668/148‑1 (MQ=255)
                                                                                                                  aaCTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTc                           <  2:2073274/148‑1 (MQ=255)
                                                                                                                  aaCTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTc                           >  2:1277982/1‑148 (MQ=255)
                                                                                                                  aaCTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTc                           >  2:1404010/1‑148 (MQ=255)
                                                                                                                    cTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCa                         >  2:1496803/1‑148 (MQ=255)
                                                                                                                     tCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttc                                  <  2:1897675/138‑1 (MQ=255)
                                                                                                                      cATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGt                       <  1:2880912/148‑1 (MQ=255)
                                                                                                                        ttGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCa                     >  2:1241215/1‑148 (MQ=255)
                                                                                                                           aTAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCg                  >  2:1799841/1‑148 (MQ=255)
                                                                                                                           aTAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCg                  >  2:1639124/1‑148 (MQ=255)
                                                                                                                            tAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGc                 >  2:2070277/1‑148 (MQ=255)
                                                                                                                            tAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGc                 >  1:101072/1‑148 (MQ=255)
                                                                                                                             aaaTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCt                <  1:289889/148‑1 (MQ=255)
                                                                                                                             aaaTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCt                >  2:2448480/1‑148 (MQ=255)
                                                                                                                               aTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCttt              <  1:1954996/148‑1 (MQ=255)
                                                                                                                               aTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCttt              <  2:794676/148‑1 (MQ=255)
                                                                                                                                  tGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCgg           <  2:1449235/148‑1 (MQ=255)
                                                                                                                                     ttttGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaaa        >  2:2754662/1‑148 (MQ=255)
                                                                                                                                     ttttGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaaa        >  2:1634941/1‑148 (MQ=255)
                                                                                                                                       ttGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa      <  2:1609236/148‑1 (MQ=255)
                                                                                                                                        tGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGAt     >  2:2382847/1‑148 (MQ=255)
                                                                                                                                          taggtaggATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCg   >  1:2756688/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
AATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGC  >  CP009273/107896‑108182

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: