Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 644,144 G→C 6.6% A292A (GCC→GCG acoA ← acetoin:2,6‑dichlorophenolindophenol oxidoreductase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947644,1440GC 6.6% 598.9 / 22.4 181A292A (GCC→GCGacoAacetoin:2,6‑dichlorophenolindophenol oxidoreductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (91/78);  new base C (5/7);  total (96/85)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

ATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGAAATGGCCGTA  >  NC_002947/643995‑644293
                                                                                                                                                     |                                                                                                                                                     
aTCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCgg                                                                                                                                                       <  2:2721178/150‑1 (MQ=255)
aTCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCgg                                                                                                                                                       <  2:135222/150‑1 (MQ=255)
aTCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCgg                                                                                                                                                       >  1:4895231/1‑150 (MQ=255)
  ctttcttGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGc                                                                                                                                                     >  2:4894859/1‑150 (MQ=255)
   tttcttGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCg                                                                                                                                                    >  2:2347453/1‑150 (MQ=255)
    ttcttGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGt                                                                                                                                                   <  2:2638712/150‑1 (MQ=255)
       ttGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGa                                                                                                                                                >  2:3287737/1‑150 (MQ=255)
        tGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCCCGCGTCGag                                                                                                                                               <  2:1459243/150‑2 (MQ=255)
         gCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATg                                                                                                                                              <  2:1995500/150‑1 (MQ=255)
         gCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATg                                                                                                                                              <  1:2899962/150‑1 (MQ=255)
         gCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCCCGCGTCGagg                                                                                                                                              <  1:127857/150‑3 (MQ=255)
         gCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCCCGCGTCGagg                                                                                                                                              <  2:2892776/150‑3 (MQ=255)
           cATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCATCGACGCCCGCGTCGaggac                                                                                                                                            <  1:4784818/150‑5 (MQ=255)
                tggtTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCa                                                                                                                                       <  2:1989604/150‑1 (MQ=255)
                 ggtTCTCCTGTTGTTCGTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAg                                                                                                                                      <  1:2820273/150‑1 (MQ=255)
                  gtTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGt                                                                                                                                     <  1:843099/150‑1 (MQ=255)
                   ttCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGtt                                                                                                                                    >  1:1694008/1‑150 (MQ=255)
                    tctcCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTg                                                                                                                                   >  1:4398168/1‑150 (MQ=255)
                       ccTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTc                                                                                                                                <  2:2735339/150‑1 (MQ=255)
                         tgttgtTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCgg                                                                                                                              >  1:62224/1‑150 (MQ=255)
                          gttgttCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGAATCGGATTTGGCGCGGCGAACGACGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGc                                                                                                                             >  2:4509369/1‑150 (MQ=255)
                           ttgttCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTc                                                                                                                                >  1:4358209/1‑146 (MQ=255)
                           ttgttCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGcc                                                                                                                            <  2:3067420/150‑1 (MQ=255)
                              ttCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTg                                                                                                                         >  2:3085767/1‑150 (MQ=255)
                              ttCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTg                                                                                                                         <  2:3190964/150‑1 (MQ=255)
                               tCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGa                                                                                                                        <  2:4454572/150‑1 (MQ=255)
                                cTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAg                                                                                                                       <  2:38006/150‑1 (MQ=255)
                                 ttGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGc                                                                                                                      >  1:1758780/1‑150 (MQ=255)
                                 ttGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGc                                                                                                                      >  1:4854933/1‑150 (MQ=255)
                                   gTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAg                                                                                                                    >  1:3715017/1‑150 (MQ=255)
                                    tAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAgg                                                                                                                   >  2:4320707/1‑150 (MQ=255)
                                     aaTCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTCGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGc                                                                                                                  >  1:3607545/1‑150 (MQ=255)
                                             gCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCgg                                                                                                          <  1:4499850/150‑1 (MQ=255)
                                                cAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGggtg                                                                                                       <  1:4720164/150‑1 (MQ=255)
                                                   gggTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTc                                                                                                    >  1:2144628/1‑150 (MQ=255)
                                                   gggTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTc                                                                                                    >  2:2454601/1‑150 (MQ=255)
                                                   gggTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTc                                                                                                    >  1:3121677/1‑150 (MQ=255)
                                                   gggTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTc                                                                                                    <  1:2335740/150‑1 (MQ=255)
                                                     gTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCtt                                                                                                  <  2:298579/150‑1 (MQ=255)
                                                     gTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCtt                                                                                                  <  2:2312354/150‑1 (MQ=255)
                                                       aGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCttgt                                                                                                >  2:4057582/1‑150 (MQ=255)
                                                        ggCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCttgtt                                                                                               <  1:4550996/150‑1 (MQ=255)
                                                        ggCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCttgtt                                                                                               <  2:293085/150‑1 (MQ=255)
                                                        ggCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCttgtt                                                                                               <  2:4551746/150‑1 (MQ=255)
                                                        accGACGTCTACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCttgtt                                                                                               >  2:2352256/3‑150 (MQ=255)
                                                              gtctaCGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGaaa                                                                                         >  1:1221137/5‑150 (MQ=255)
                                                              gTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGaaa                                                                                         >  1:4791978/1‑150 (MQ=255)
                                                                aGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAAtt                                                                                       <  2:1695540/150‑1 (MQ=255)
                                                                aGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAAtt                                                                                       >  2:342336/1‑150 (MQ=255)
                                                                aGACGTCGGTGAGCAGGTCGGCGGGCTGGGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAAtt                                                                                       <  2:4365686/150‑1 (MQ=255)
                                                                  aCGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGACCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGc                                                                                     >  2:2653186/1‑150 (MQ=255)
                                                                   cGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCa                                                                                    <  1:4543711/150‑1 (MQ=255)
                                                                   cGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCa                                                                                    <  2:2415875/150‑1 (MQ=255)
                                                                    gTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGa                                                                                                                                                <  2:351640/89‑1 (MQ=255)
                                                                    gTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCAt                                                                                   >  2:674527/1‑150 (MQ=255)
                                                                    gTCGGTGAGCAGGTCAGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCAt                                                                                   >  1:3900887/1‑150 (MQ=255)
                                                                       ggTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAg                                                                                <  2:1867372/150‑1 (MQ=255)
                                                                        gTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAgg                                                                               <  1:3337248/150‑1 (MQ=255)
                                                                                gTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTcgcg                                                                       >  1:2825791/1‑150 (MQ=255)
                                                                                gTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGAAGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTcgcg                                                                       <  1:130540/150‑1 (MQ=255)
                                                                                 tCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCgg                                                                      <  1:93451/150‑1 (MQ=255)
                                                                                 tCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCgg                                                                      <  2:2668134/150‑1 (MQ=255)
                                                                                   ggcggGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGa                                                                    >  1:556260/1‑150 (MQ=255)
                                                                                      gggCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGTGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTc                                                                 <  1:512990/150‑1 (MQ=255)
                                                                                      gggCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTc                                                                 <  1:1373707/150‑1 (MQ=255)
                                                                                            tgGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGc                                                                                     <  1:3016408/124‑1 (MQ=255)
                                                                                            tgGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGTATTCGCGCAg                                                           >  1:3095731/1‑150 (MQ=255)
                                                                                            tgGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAg                                                           >  1:1311674/1‑150 (MQ=255)
                                                                                                  gggATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGttctt                                                     <  2:331005/150‑1 (MQ=255)
                                                                                                   ggATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGttcttc                                                    <  2:3832402/150‑1 (MQ=255)
                                                                                                     aTCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTcttcac                                                  >  2:2435677/1‑150 (MQ=255)
                                                                                                      tCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTcttcact                                                 >  1:382871/1‑150 (MQ=255)
                                                                                                        ggATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTcttcacttc                                               >  1:1410001/1‑150 (MQ=255)
                                                                                                        ggATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTcttcacttc                                               <  1:3518762/150‑1 (MQ=255)
                                                                                                         gATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTcttcacttca                                              >  1:2289217/1‑150 (MQ=255)
                                                                                                         gATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTcttcacttca                                              >  1:2289277/1‑150 (MQ=255)
                                                                                                           tctGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTtcatc                                            <  1:2675964/148‑1 (MQ=255)
                                                                                                              ggCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAgg                                         <  1:1885825/150‑1 (MQ=255)
                                                                                                              ggCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAgg                                         <  2:4254763/150‑1 (MQ=255)
                                                                                                                cgcgGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAgggg                                       >  1:2489399/1‑150 (MQ=255)
                                                                                                                     cgAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCgg                                  >  1:773993/1‑150 (MQ=255)
                                                                                                                      gAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGt                                 >  2:5030257/1‑150 (MQ=255)
                                                                                                                       aaCGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTa                                <  1:378603/150‑1 (MQ=255)
                                                                                                                              tCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTgc                         <  1:4048995/150‑1 (MQ=255)
                                                                                                                                   gATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATc                    <  1:679474/150‑1 (MQ=255)
                                                                                                                                   gATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATc                    <  1:5074766/150‑1 (MQ=255)
                                                                                                                                   gATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATc                    <  1:2347453/150‑1 (MQ=255)
                                                                                                                                        ggTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCtt               >  2:3622325/1‑150 (MQ=255)
                                                                                                                                         gccatCGACGCCCGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTc              >  2:4899776/5‑150 (MQ=255)
                                                                                                                                         gTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGGTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTc              <  2:1268267/150‑1 (MQ=255)
                                                                                                                                           catCGACGCCCGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGa            >  1:1910581/3‑150 (MQ=255)
                                                                                                                                            cTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGaa           >  2:116081/1‑150 (MQ=255)
                                                                                                                                            cTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGaa           >  1:498173/1‑150 (MQ=255)
                                                                                                                                            atCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGaa           >  2:842039/2‑150 (MQ=255)
                                                                                                                                            atCGACGCCCGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGaa           >  2:957623/2‑150 (MQ=255)
                                                                                                                                             tCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGaaa          >  1:3044900/1‑150 (MQ=255)
                                                                                                                                             tCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGaaa          >  1:3891151/1‑150 (MQ=255)
                                                                                                                                              cGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGAAAt         <  2:40059/150‑1 (MQ=255)
                                                                                                                                                aCGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGAAATgg       <  2:1099510/150‑1 (MQ=255)
                                                                                                                                                     ggCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGAAATGGCCgta  <  2:4398168/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTTTCTTGCCATGGTGGTTCTCCTGTTGTTCTTGTAATCCGCGGCTCAGGGGTAGGCGACGTAGACGTCGGTGAGCAGGTCGGCGGGCTGTGGCTTGGGATCGGATTTGGCGCGGCGAACGGCGTCTTCGATCAGGTCCTCGACGCGGGCGTCGATGGCGTCCAGTTGCTCGGCCGTGAGCAGGCCGGCGCGGGTGGTCTTGTTGCGAAATTGCATGAGGCAGTCGCGGGATTCGCGCAGGTTCTTCACTTCATCAGGGGCGCGGTAGGTTTGCGCATCGCCTTCGAAATGGCCGTA  >  NC_002947/643995‑644293

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: