Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,567,255 A→T 5.7% intergenic (+59/‑262) 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,2550AT 5.7% 274.9 / 4.8 87intergenic (+59/‑262)pcaR/pcaKpca regulon transcriptional regulator/4‑hydroxybenzoate transporter
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base T (3/2);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GGGCGGGTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCAC  >  NC_002947/1567107‑1567400
                                                                                                                                                    |                                                                                                                                                  
gggcgggTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGa                                                                                                                                                    >  6:266540/1‑149 (MQ=255)
 ggcgggTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGa                                                                                                                                                    >  2:75149/1‑148 (MQ=255)
 ggcgggTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGAt                                                                                                                                                   <  3:346872/149‑1 (MQ=255)
    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                                                                                                                                >  4:228338/1‑141 (MQ=255)
                            agcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGcc                                                                                                                                <  3:228338/141‑1 (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                                                                                                  >  5:57026/1‑149 (MQ=255)
                                                 tGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGa                                                                                                  >  6:7856/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                                                                 >  1:56852/1‑149 (MQ=255)
                                                                                  ttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgacc                                                                 >  4:371970/1‑149 (MQ=255)
                                                                                     gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                              <  7:189802/149‑3 (MQ=255)
                                                                                     gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                              <  6:57026/149‑3 (MQ=255)
                                                                                     gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                              <  8:46047/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                            >  5:265571/1‑149 (MQ=255)
                                                                                                                       tCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATa                            <  1:13118/149‑1 (MQ=255)
                                                                                                                         aGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATAc                           <  3:87628/148‑1 (MQ=255)
                                                                                                                         aGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACt                          <  6:372427/149‑1 (MQ=255)
                                                                                                                          gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGcactca                                      >  7:94890/1‑136 (MQ=255)
                                                                                                                          gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGcactca                                      <  8:94890/136‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
GGGCGGGTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCAC  >  NC_002947/1567107‑1567400

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: