Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,667,310 A→C 6.0% intergenic (+37/+34) nhaP → / ← PP_4975 NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,667,3100AC 6.0% 355.2 / 12.1 126intergenic (+37/+34)nhaP/PP_4975NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base A (61/57);  new base C (4/4);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 126 total aligned reads displayed.

GACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACA  >  NC_002947/5667172‑5667455
                                                                                                                                          |                                                                                                                                                  
gACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAaga                                                                                                                                         >  1:1282041/1‑150 (MQ=255)
  ccTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAgg                                                                                                                                       <  2:1376240/150‑1 (MQ=255)
   ctgctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCTCTTCGCGGGTTCACCCGCGaatgggc                                                                                                                                      <  1:1847485/150‑6 (MQ=255)
   ctgctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaatgggc                                                                                                                                      <  2:301540/150‑6 (MQ=255)
   ctgctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGc                                                                                                                                      >  2:2062806/1‑150 (MQ=255)
    tgctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGcc                                                                                                                                     >  2:2052703/1‑150 (MQ=255)
     gctgctgTCGATCACCTACATCGTGGTGTTGTCGGCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCg                                                                                                                                    <  2:99381/150‑1 (MQ=255)
         ctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGcac                                                                                                                                >  2:1426882/1‑150 (MQ=255)
         ctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGcac                                                                                                                                >  2:344934/1‑150 (MQ=255)
            tCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAgg                                                                                                                             >  2:235093/1‑150 (MQ=255)
            tCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAgg                                                                                                                             >  2:955536/1‑150 (MQ=255)
              gATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCa                                                                                                                           >  1:1775095/1‑150 (MQ=255)
                tCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAAc                                                                                                                         <  1:2052703/150‑1 (MQ=255)
                    cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     <  1:955536/150‑1 (MQ=255)
                    cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     <  1:1800452/150‑1 (MQ=255)
                    cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     <  1:618395/150‑1 (MQ=255)
                    cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGAAACCAgg                                                                                                                     >  1:78829/1‑150 (MQ=255)
                      aCATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCTATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                   <  1:480667/150‑1 (MQ=255)
                       cATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     <  2:2160632/147‑1 (MQ=255)
                                   gtcgtcGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                                     >  1:1041709/1‑150 (MQ=255)
                                   gtcgtcGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                                     <  2:1282041/150‑1 (MQ=255)
                                    tcgtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGa                                                                                                          >  2:1846734/1‑144 (MQ=255)
                                    tcgtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGa                                                                                                          <  1:1846734/144‑1 (MQ=255)
                                      gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAg                                                                                                                              >  1:1043635/1‑123 (MQ=255)
                                      gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAg                                                                                                                              <  2:1043635/123‑1 (MQ=255)
                                      gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGaa                                                                                                  >  1:1243160/1‑150 (MQ=255)
                                         gATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGAc                                                                                                         >  1:605897/1‑140 (MQ=255)
                                         gATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGAc                                                                                                         <  2:605897/140‑1 (MQ=255)
                                           ttCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCgg                                                                                                                                   <  1:424472/113‑1 (MQ=255)
                                           ttCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCgg                                                                                                                                   >  2:424472/1‑113 (MQ=255)
                                           ttCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                             <  2:1775095/150‑1 (MQ=255)
                                           ttCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                             >  2:2290176/1‑150 (MQ=255)
                                                gTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                                        >  2:39786/1‑134 (MQ=255)
                                                gTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                                        <  1:39786/134‑1 (MQ=255)
                                                 tGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTg                                                                                       <  2:1337180/150‑1 (MQ=255)
                                                  gCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTgg                                                                                      <  2:650989/150‑1 (MQ=255)
                                                  gCAGGGGTGGAGCATTGGGGGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAAGGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAGCCCGCGAAGAGGCCGGCACAGGCAACCAGGGTGTATTCGACTGCCGAATCCGGGAACTgg                                                                                      <  2:69739/150‑1 (MQ=255)
                                                    aGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGACCTGGTc                                                                                    <  2:1243160/150‑1 (MQ=255)
                                                      gggTTGAGAATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCt                                                                                  <  1:2175749/150‑1 (MQ=255)
                                                          tGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCAc                                                                              <  2:2347941/150‑1 (MQ=255)
                                                             gCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTa                                                                           >  2:1007950/1‑150 (MQ=255)
                                                              cATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTaa                                                                          >  1:234650/1‑149 (MQ=255)
                                                              cATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                          <  2:680897/150‑1 (MQ=255)
                                                               aTTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTc                                                                                                                  >  2:123210/1‑110 (MQ=255)
                                                               aTTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTc                                                                                                                  <  1:123210/110‑1 (MQ=255)
                                                               aTTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAAGGCCGGGTCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACt                                                                         >  1:1862713/1‑150 (MQ=255)
                                                                  gggCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCAGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGa                                                                      <  1:999057/150‑1 (MQ=255)
                                                                      gggTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                          >  1:1446298/1‑142 (MQ=255)
                                                                      gggTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                          <  2:1446298/142‑1 (MQ=255)
                                                                      gggTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGAtctc                                                                  >  2:1452627/1‑150 (MQ=255)
                                                                       ggtggtGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                 <  2:2192582/150‑1 (MQ=255)
                                                                        gtggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGaa                                                                                          <  1:2191210/124‑1 (MQ=255)
                                                                        gtggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCgg                                                                <  1:1007950/150‑1 (MQ=255)
                                                                           gtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                   >  1:1752480/1‑97 (MQ=255)
                                                                           gtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                   <  2:1752480/97‑1 (MQ=255)
                                                                           gtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTgc                                                             <  1:2062806/150‑1 (MQ=255)
                                                                                   ggTCAGCGCCCGGGCTTGAGAATGCCGGGGCCGGTTTGCGGCCGATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATgg                                                     <  2:350492/150‑1 (MQ=255)
                                                                                   ggTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     >  1:658108/1‑87 (MQ=255)
                                                                                   ggTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     <  2:658108/87‑1 (MQ=255)
                                                                                     tCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATgg                                                     <  1:1198237/148‑1 (MQ=255)
                                                                                     tCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGc                                                   >  2:1525084/1‑150 (MQ=255)
                                                                                     ccAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGc                                                   >  2:1654173/2‑150 (MQ=255)
                                                                                      cAGCGCCCAGCGTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGCCTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCt                                                  <  1:238763/150‑1 (MQ=255)
                                                                                       aGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                             <  1:42163/106‑1 (MQ=255)
                                                                                       aGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                             >  2:42163/1‑106 (MQ=255)
                                                                                       aGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTg                                                 >  1:511344/1‑150 (MQ=255)
                                                                                             cAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGTGACCTTGCGCTTCCGGCa                                           >  2:1870310/1‑150 (MQ=39)
                                                                                             cAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCa                                           >  2:322608/1‑150 (MQ=255)
                                                                                               gCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAg                                         >  2:1781265/1‑150 (MQ=255)
                                                                                                ccTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGGTCTCGGTGCGCCCATGGGCTGCCGGCAAGc                                        >  1:2062771/1‑150 (MQ=255)
                                                                                                ccTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGc                                        <  1:1606830/150‑1 (MQ=255)
                                                                                                  ttGAGAATGCCGGGGCCGCTTTGCGGCCCCTTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCa                                      >  1:2183884/1‑150 (MQ=255)
                                                                                                  ttGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTACTTCTCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGAGACCATGGGCTGCCGGCAAGCCa                                      >  1:1955851/1‑150 (MQ=38)
                                                                                                  ttGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCa                                      >  2:100158/1‑150 (MQ=255)
                                                                                                        aTGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCg                                <  1:1532111/150‑1 (MQ=255)
                                                                                                        aTGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCg                                >  1:254708/1‑150 (MQ=255)
                                                                                                         tGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGGACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGc                               <  1:1870310/150‑1 (MQ=255)
                                                                                                                gCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGtt                        <  1:1525084/150‑1 (MQ=255)
                                                                                                                 ccGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTc                                 <  2:2214810/140‑1 (MQ=255)
                                                                                                                 ccGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTc                                 >  1:2214810/1‑140 (MQ=255)
                                                                                                                  cGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCAGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGGGCGCCCATGGGCTTCCGGCAAGCCATCCTCGCCCAGGTTGa                      >  2:2280410/1‑150 (MQ=255)
                                                                                                                   gCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                             <  1:2340144/78‑1 (MQ=255)
                                                                                                                   gCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc                     <  1:344934/150‑1 (MQ=255)
                                                                                                                      ttacGGCTCATTCGCGGGTGAACCCCGGAAGAGGCCGGCCCAGGCAACCAGGGTGTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGaa                  <  2:254708/147‑1 (MQ=255)
                                                                                                                       ggCGGCCCATGTGGGGGGGATCCCGCGAAGAGGGCGGCACGGGGAACCAGGGTCTACTCGACGGCCGAATCCGGGAACTGGTCCTGCACGTGCTTGGTCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAac                 <  2:234650/149‑1 (MQ=255)
                                                                                                                         cGGCCTCTTCGCGGGTTCACCCGCGAATGGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCaa                                          <  2:1203840/123‑1 (MQ=38)
                                                                                                                         cGGCCCATTGGCGGGGGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAacac               <  1:725443/150‑1 (MQ=255)
                                                                                                                          ggCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACAcc              >  1:944825/1‑150 (MQ=255)
                                                                                                                           gCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCa             >  1:1292472/1‑150 (MQ=255)
                                                                                                                            cccATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCAt            >  2:1480318/1‑150 (MQ=255)
                                                                                                                             ctcttCGCGGGTTCACCCGCGAATGGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTgc                                                             >  1:949380/4‑100 (MQ=38)
                                                                                                                             ctcttCGCGGGTTCACCCGCGAATGGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTgc                                                             <  2:949380/97‑1 (MQ=38)
                                                                                                                              tcttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCt          >  2:1403241/3‑150 (MQ=255)
                                                                                                                              cATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGAGCCCATGGGCTGCCGGCGAGCCATCCTCGCCCAGGTTGACGAAGACCATCt          >  1:35203/1‑150 (MQ=255)
                                                                                                                                ttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTg        >  2:2157002/1‑150 (MQ=255)
                                                                                                                                ttCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTg        <  1:235093/150‑1 (MQ=255)
                                                                                                                                  cgcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTc      <  1:1654173/150‑1 (MQ=255)
                                                                                                                                   gcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTg                                                                                       >  2:190263/1‑68 (MQ=255)
                                                                                                                                      ggTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCATGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACa  <  1:2290176/150‑1 (MQ=255)
                                                                                                                                      ggTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACTAacac               <  1:823538/137‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACA  >  NC_002947/5667172‑5667455

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