Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,573,906 G→A 5.4% intergenic (+52/‑275) pcaK → / → pcaF‑I 4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,573,9060GA 5.4% 939.2 / 26.7 294intergenic (+52/‑275)pcaK/pcaF‑I4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (130/148);  new base A (11/5);  total (141/153)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 309 total aligned reads displayed.

TTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACCC  >  NC_002947/1573759‑1574055
                                                                                                                                                   |                                                                                                                                                     
ttCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggc                                                                                                                                                     >  2:6311582/1‑150 (MQ=255)
ttCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggc                                                                                                                                                     >  2:3268337/1‑150 (MQ=255)
  cGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggcgg                                                                                                                                                   <  2:2896335/150‑1 (MQ=255)
   gAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggcggc                                                                                                                                                  <  1:4841484/150‑1 (MQ=255)
     gCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGt                                                                                                                                                >  2:2689066/1‑150 (MQ=255)
      cAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGCCCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGtt                                                                                                                                               <  1:5106698/150‑1 (MQ=255)
      cAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGtt                                                                                                                                               <  2:1997635/150‑1 (MQ=255)
        ggTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgt                                                                                                                                             <  1:334256/150‑1 (MQ=255)
         gTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                            >  1:4660076/1‑150 (MQ=255)
         gTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                            >  2:5112783/1‑150 (MQ=255)
          tGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgc                                                                                                                                           <  1:3008889/150‑1 (MQ=255)
          tGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgc                                                                                                                                           <  1:6198343/150‑1 (MQ=255)
           gCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGAGCGGCGttacgct                                                                                                                                          <  1:101106/150‑6 (MQ=38)
               aCGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCaa                                                                                                                                      >  1:118478/1‑150 (MQ=255)
                  gCGTTGCTGGTACCGGCGGCGCTGGCGACGGTGGGAGTGGTCGTCAAGGGGCTGGCCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCa                                                                                                                                   <  2:2446516/150‑1 (MQ=255)
                   cGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAg                                                                                                                                  <  1:2494576/150‑1 (MQ=255)
                    gTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAgg                                                                                                                                 <  1:4407206/150‑1 (MQ=255)
                         tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                                            >  1:2916337/1‑150 (MQ=255)
                         tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                                            <  1:893052/150‑1 (MQ=255)
                         tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                                            >  1:5003958/1‑150 (MQ=255)
                         tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                                            >  2:5766356/1‑150 (MQ=255)
                           gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGcaca                                                                                                                          <  1:983404/150‑1 (MQ=255)
                           gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCAGCGGTGTGCTGCAATCAGGGCGcaca                                                                                                                          >  2:4272910/1‑150 (MQ=255)
                            tACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAg                                                                                                                         >  1:60238/1‑150 (MQ=255)
                              ttttcggcgCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGtt                                                                                                                       >  1:6799772/5‑150 (MQ=255)
                              ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGtt                                                                                                                       >  2:4120392/1‑150 (MQ=255)
                              ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGtt                                                                                                                       <  1:4010149/150‑1 (MQ=255)
                               cggcggCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTg                                                                                                                      <  1:1376689/150‑1 (MQ=255)
                                 gcggcgCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGcc                                                                                                                    <  1:4496849/150‑1 (MQ=255)
                                  cggcgCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGccc                                                                                                                   >  1:1720318/1‑150 (MQ=255)
                                      gcTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgcg                                                                                                               >  1:4254821/1‑150 (MQ=255)
                                        tGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGcgcg                                                                                                             >  2:3022525/1‑150 (MQ=255)
                                         ggCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGt                                                                                                            >  2:1795478/1‑150 (MQ=255)
                                          gCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTg                                                                                                           >  1:2276441/1‑150 (MQ=255)
                                          gCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTg                                                                                                           >  2:3197786/1‑150 (MQ=255)
                                            gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGgggg                                                                                                         >  1:2213588/1‑150 (MQ=255)
                                             aCCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgggg                                                                                                        >  1:2157891/1‑150 (MQ=255)
                                             aCCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgggg                                                                                                        <  2:1984548/150‑1 (MQ=255)
                                               cGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTg                                                                                                      <  1:3452464/150‑1 (MQ=255)
                                                 tGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                                    >  1:3101740/1‑150 (MQ=255)
                                                  gggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGAt                                                                                                   <  2:1088042/150‑1 (MQ=255)
                                                   ggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATa                                                                                                  <  1:1456848/150‑1 (MQ=255)
                                                    gAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAg                                                                                                 >  2:2371698/1‑150 (MQ=255)
                                                     aGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGt                                                                                                >  1:2006735/1‑150 (MQ=255)
                                                               aaGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaacaa                                                                                      >  2:1215252/1‑150 (MQ=255)
                                                                 ggggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAAt                                                                                    <  2:4399141/150‑1 (MQ=255)
                                                                 ggggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAAt                                                                                    >  1:4804391/1‑150 (MQ=255)
                                                                  gggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATg                                                                                   <  1:3254357/150‑1 (MQ=255)
                                                                   ggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCTACAAATgt                                                                                  <  2:2191283/150‑2 (MQ=255)
                                                                     cTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                                                <  1:152910/150‑1 (MQ=255)
                                                                     cTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                                                <  1:5141684/150‑1 (MQ=255)
                                                                     cTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                                                >  2:6088891/1‑150 (MQ=255)
                                                                      tGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTc                                                                               <  1:279574/150‑1 (MQ=255)
                                                                        gTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                                                             <  2:3456205/150‑1 (MQ=255)
                                                                         tCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGtt                                                                            <  2:1720318/150‑1 (MQ=255)
                                                                          cAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTa                                                                           <  1:458490/150‑1 (MQ=255)
                                                                             ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                                                        <  1:353840/150‑1 (MQ=255)
                                                                                cGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                                                     <  2:2575526/150‑1 (MQ=255)
                                                                                 gCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCGGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGcac                                                                    >  1:1948969/1‑150 (MQ=255)
                                                                                    gACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGa                                                                 <  1:6645147/150‑1 (MQ=255)
                                                                                     aCGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGAt                                                                >  2:2782786/1‑150 (MQ=255)
                                                                                         aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgctg                                                                                                                                         <  1:5506028/73‑1 (MQ=255)
                                                                                         aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgctg                                                                                                                                         >  2:5506028/1‑73 (MQ=255)
                                                                                         aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTc                                                            >  2:1165029/1‑150 (MQ=255)
                                                                                         aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTc                                                            <  1:1831709/150‑1 (MQ=255)
                                                                                          aCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCg                                                           >  2:155447/1‑150 (MQ=255)
                                                                                           ccTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGa                                                          <  2:2240389/150‑1 (MQ=255)
                                                                                            cTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAt                                                         >  2:5910083/1‑150 (MQ=255)
                                                                                             tGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAtt                                                        >  1:5670977/1‑150 (MQ=255)
                                                                                              gagggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTa                                                       >  1:848068/1‑150 (MQ=255)
                                                                                                  gcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCACACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCgg                                                   <  2:1459747/150‑1 (MQ=255)
                                                                                                    gagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAt                                                 <  1:365863/150‑1 (MQ=255)
                                                                                                     agggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt                                                <  1:6311582/150‑1 (MQ=255)
                                                                                                       ggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGt                                              >  1:2961138/1‑150 (MQ=255)
                                                                                                         ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTaa                                            <  2:861728/150‑1 (MQ=255)
                                                                                                         ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTaa                                            <  1:3762558/150‑1 (MQ=255)
                                                                                                           gcgTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAgg                                          <  2:5052631/150‑1 (MQ=255)
                                                                                                               tGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacc                                      <  2:4506810/150‑1 (MQ=255)
                                                                                                                gCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacct                                     <  2:2523747/150‑1 (MQ=255)
                                                                                                                   gCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacctacc                                  >  1:5849373/1‑150 (MQ=255)
                                                                                                                     ccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAg                                >  2:4971200/1‑150 (MQ=255)
                                                                                                                      cATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGt                               >  1:6174008/1‑150 (MQ=255)
                                                                                                                          cacaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTcc                           <  1:2405907/150‑1 (MQ=255)
                                                                                                                           acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCt                          >  2:3374053/1‑150 (MQ=255)
                                                                                                                             agcacaaCGCCGCCCCTGCGGGAGCGGCGTTACGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTa                        >  1:4135243/6‑150 (MQ=38)
                                                                                                                                 tAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCt                    <  2:5327408/150‑1 (MQ=255)
                                                                                                                                  aaCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTg                   <  2:3824234/150‑1 (MQ=255)
                                                                                                                                  aaCGCCGCCCCTGCGGGAGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAAGGGTTCGTTAATCGCACAGAGAGGCGAGGAGCGGCTTGTACGGTCCCTACCCGTCTCCTTCAGCTg                   >  2:6459422/1‑150 (MQ=255)
                                                                                                                                   aCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTTGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                  <  2:295636/150‑1 (MQ=255)
                                                                                                                                   aCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                  <  2:6664144/150‑1 (MQ=255)
                                                                                                                                   aCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                  <  2:6799772/150‑1 (MQ=255)
                                                                                                                                   aCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                  <  2:1403068/150‑1 (MQ=255)
                                                                                                                                     gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCACCAAATGGTTCGTTAATCGCACAGATAGTCGAGTATCGGATTGTCCGGTACCGACCAGTCTCCTGAATCTGGGt                >  2:4065598/1‑150 (MQ=255)
                                                                                                                                        gcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTTGGTTTa             >  1:3304909/1‑150 (MQ=255)
                                                                                                                                          tCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTAcc           >  1:1277644/1‑150 (MQ=255)
                                                                                                                                             cGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGaa        >  1:6393636/1‑150 (MQ=255)
                                                                                                                                             cGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGaa        <  1:3887786/150‑1 (MQ=255)
                                                                                                                                             cGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGaa        <  1:2696062/150‑1 (MQ=255)
                                                                                                                                              gCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGGCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAg       >  2:3554886/1‑150 (MQ=255)
                                                                                                                                                  gggCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCAcc   <  2:1595817/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACCC  >  NC_002947/1573759‑1574055

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: