Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,567,259 A→T 5.9% intergenic (+63/‑258) pcaR → / → pcaK pca regulon transcriptional regulator/4‑hydroxybenzoate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,567,2590AT 5.9% 256.1 / 4.8 85intergenic (+63/‑258)pcaR/pcaKpca regulon transcriptional regulator/4‑hydroxybenzoate transporter
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base T (3/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

GGGTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCATG  >  NC_002947/1567111‑1567407
                                                                                                                                                    |                                                                                                                                                     
gggTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATaaa                                                                                                                                                       <  4:184750/149‑1 (MQ=255)
 ggTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAAc                                                                                                                                                      <  1:115508/149‑1 (MQ=255)
     aCCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGca                                                                                                                                                   <  8:15646/148‑1 (MQ=255)
      cccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGTTTATCGcaca                                                                                                                                                 <  6:376726/149‑1 (MQ=255)
      cccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGcaca                                                                                                                                                 <  1:176983/149‑1 (MQ=255)
       ccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGTTTATCGCACAg                                                                                                                                                <  8:243864/149‑1 (MQ=255)
       ccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTAGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAg                                                                                                                                                >  6:154114/1‑149 (MQ=255)
        cGCATTGAACTCGAGCAGCGGTTCTTGCCGAATTTGCTCGCAGCCACCCGGGACCATTCACCACAGATGTTTGGTTGAGACCGCAGCCAGCACTATCCGGCCCCGGAACAGAACAGGCAAACCCGTTACCTCTGACATAAACCAACAg                                                                                                                                                >  8:47652/1‑148 (MQ=255)
               aaCTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGc                                                                                                                                        >  1:52463/1‑149 (MQ=255)
                    gAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCggg                                                                                                                                   <  4:6088/148‑1 (MQ=255)
                    gAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCgg                                                                                                                                    <  1:186378/147‑1 (MQ=255)
                     agcagcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTg                                                                                                                                 <  7:290351/149‑1 (MQ=255)
                      gcagcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTg                                                                                                                                 >  1:340775/1‑148 (MQ=255)
                        agcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGcc                                                                                                                                       <  3:228338/141‑1 (MQ=255)
                        agcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGcc                                                                                                                                       >  4:228338/1‑141 (MQ=255)
                        agcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAAt                                                                                                                              >  4:239995/1‑149 (MQ=255)
                        agcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAAt                                                                                                                              <  1:263918/149‑1 (MQ=255)
                           ggTTCTTGGCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGc                                                                                                                           <  2:340775/149‑1 (MQ=255)
                             ttCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGt                                                                                                                         >  3:159942/1‑149 (MQ=255)
                                ttGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCAc                                                                                                                      >  8:189802/1‑149 (MQ=255)
                                  gCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCa                                                                                                                       <  3:236479/146‑1 (MQ=255)
                                   ccGCTTTTGCTGGCAGCCAGCCGGGACCTGTGCCACCAGGTGTTGGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATACACGCACAGTGTCGTGCCGGGTGAAGTGCGTCCCCTc                                                                                                                   >  6:298103/1‑149 (MQ=255)
                                   ccGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTc                                                                                                                   >  5:215824/1‑149 (MQ=255)
                                     gATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTcgg                                                                                                                 <  6:215824/149‑1 (MQ=255)
                                         ttGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGCACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTcggccg                                                                                                              <  5:266540/148‑1 (MQ=255)
                                            cTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCt                                                                                                           >  7:46047/1‑148 (MQ=255)
                                             tGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGa                                                                                                         >  6:7856/1‑149 (MQ=255)
                                             tGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGa                                                                                                         >  5:57026/1‑149 (MQ=255)
                                               gCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGAt                                                                                                        >  4:202340/1‑148 (MQ=255)
                                                  gccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACAGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAAt                                                                                                    <  6:283119/149‑1 (MQ=255)
                                                    cagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGt                                                                                                  >  4:87628/1‑149 (MQ=255)
                                                     agccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGtc                                                                                                 >  1:373188/1‑149 (MQ=255)
                                                         gggACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGtctctc                                                                                             <  1:75149/149‑1 (MQ=255)
                                                         gggACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGtctctc                                                                                             >  8:71881/1‑149 (MQ=255)
                                                           gACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCg                                                                                            >  8:345263/1‑148 (MQ=255)
                                                                   ccaccaGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCg                                                                                   >  3:368830/1‑149 (MQ=255)
                                                                   ccaccaGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCg                                                                                   >  8:60609/1‑149 (MQ=255)
                                                                      ccaGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGc                                                                                 >  2:323935/1‑148 (MQ=255)
                                                                         gttgttTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGc                                                                             <  3:202340/149‑1 (MQ=255)
                                                                            gttTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGCCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccga                                                                          >  2:325583/1‑149 (MQ=255)
                                                                             tttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgac                                                                         <  4:159942/149‑1 (MQ=255)
                                                                              ttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgacc                                                                        >  4:371970/1‑149 (MQ=255)
                                                                              ttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgacc                                                                        >  1:56852/1‑149 (MQ=255)
                                                                                 gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                                     <  8:46047/149‑3 (MQ=255)
                                                                                 gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                                     <  7:189802/149‑3 (MQ=255)
                                                                                 gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                                     <  6:57026/149‑3 (MQ=255)
                                                                                 gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGaa                                                                     >  2:13118/1‑149 (MQ=255)
                                                                                 gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGaa                                                                     <  7:345263/149‑1 (MQ=255)
                                                                                     agaCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACaa                                                                  <  7:60609/148‑1 (MQ=255)
                                                                                      gaCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAg                                                                 >  4:79359/1‑148 (MQ=255)
                                                                                      gaCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAg                                                                 >  6:230549/1‑148 (MQ=255)
                                                                                       aCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGa                                                                >  3:362037/1‑148 (MQ=255)
                                                                                       aCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATCTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaa                                                               <  4:333852/149‑1 (MQ=255)
                                                                                        ccGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaa                                                               >  5:372427/1‑148 (MQ=255)
                                                                                         cGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaaa                                                              >  6:285072/1‑148 (MQ=255)
                                                                                          gCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaaaaa                                                            <  5:7856/149‑1 (MQ=255)
                                                                                              ggAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAacac                                                         >  1:97867/1‑148 (MQ=255)
                                                                                              ggAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACt                                                        <  8:140541/149‑1 (MQ=255)
                                                                                                            ccccGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCactcactc                                          <  2:52463/149‑1 (MQ=255)
                                                                                                              ccGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCAt                                        <  6:166754/149‑1 (MQ=255)
                                                                                                               cGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATg                                       <  2:373188/149‑1 (MQ=255)
                                                                                                                ggATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGa                                      <  5:154114/149‑1 (MQ=255)
                                                                                                                  aTCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACACACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAAt                                    <  6:138523/149‑1 (MQ=255)
                                                                                                                   tCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATa                                   <  1:13118/149‑1 (MQ=255)
                                                                                                                   tCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATa                                   >  5:265571/1‑149 (MQ=255)
                                                                                                                     aGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATAc                                  <  3:87628/148‑1 (MQ=255)
                                                                                                                     aGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACt                                 <  6:372427/149‑1 (MQ=255)
                                                                                                                      gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGcactca                                             <  8:94890/136‑1 (MQ=255)
                                                                                                                      gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGcactca                                             >  7:94890/1‑136 (MQ=255)
                                                                                                                      gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTg                                >  2:263404/1‑149 (MQ=255)
                                                                                                                      gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCCGCTGATTAATGTCTCTCGCAACTGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTAATCAAAACTg                                >  6:296132/1‑149 (MQ=255)
                                                                                                                        aCAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTc                              <  3:371970/149‑1 (MQ=255)
                                                                                                                           ggCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCAcc                           >  1:363020/1‑149 (MQ=255)
                                                                                                                           ggCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATAATGTCAcc                           >  4:168839/1‑149 (MQ=255)
                                                                                                                              aaaaCCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTAACCTc                         >  6:252127/1‑148 (MQ=255)
                                                                                                                               aaaCCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCgc                       <  3:239995/149‑1 (MQ=255)
                                                                                                                                        acTGTGCGTTTATCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGc              >  4:46634/2‑149 (MQ=255)
                                                                                                                                         cTGTGCGTTTATCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGcc             >  7:159944/1‑149 (MQ=255)
                                                                                                                                         cTGTGCGTTTATCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGcc             >  7:315301/1‑149 (MQ=255)
                                                                                                                                         cTGTGCGATAAACGCACAGTGTCGTGCCGGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGaa                                                                     >  1:23395/1‑94 (MQ=255)
                                                                                                                                         cTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGaa                                                                     <  2:23395/93‑1 (MQ=255)
                                                                                                                                         cTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAac                                                           >  2:136709/1‑103 (MQ=255)
                                                                                                                                         cTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAac                                                           <  1:136709/103‑1 (MQ=255)
                                                                                                                                              cGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCAc         <  5:230549/148‑1 (MQ=255)
                                                                                                                                               gATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATa                                   >  2:237791/1‑121 (MQ=255)
                                                                                                                                               gATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATa                                   <  1:237791/121‑1 (MQ=255)
                                                                                                                                                  aaaCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGc     <  1:323935/148‑1 (MQ=255)
                                                                                                                                                  aaaCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCa    >  8:302762/1‑149 (MQ=255)
                                                                                                                                                   aaCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCAt   <  7:71881/149‑1 (MQ=255)
                                                                                                                                                   aaCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCAt   >  6:179541/1‑149 (MQ=255)
                                                                                                                                                   aaCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCAt   <  2:97867/149‑1 (MQ=255)
                                                                                                                                                    aCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCATg  <  6:265571/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCATG  >  NC_002947/1567111‑1567407

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