Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,567,259 A→T 5.8% 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.8% 200.7 / 2.9 69intergenic (+63/‑258)pcaR/pcaKpca regulon transcriptional regulator/4‑hydroxybenzoate transporter
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

GGTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCA  >  NC_002947/1567112‑1567405
                                                                                                                                                   |                                                                                                                                                  
ggTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAAc                                                                                                                                                   <  7:274650/149‑1 (MQ=255)
   gACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGca                                                                                                                                                <  3:337801/149‑1 (MQ=255)
    aCCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGca                                                                                                                                                >  6:287780/1‑148 (MQ=255)
      ccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGTTTATCGCACAg                                                                                                                                             <  6:226332/149‑1 (MQ=255)
      ccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAACCGCACAg                                                                                                                                             >  8:135955/1‑149 (MQ=255)
      ccGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAg                                                                                                                                             >  5:330100/1‑149 (MQ=255)
       cGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAgt                                                                                                                                            >  2:99145/1‑149 (MQ=255)
          aGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACAAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTc                                                                                                                                         <  8:223373/149‑1 (MQ=255)
              aaCTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGc                                                                                                                                     >  3:90502/1‑149 (MQ=255)
                   gAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGt                                                                                                                                >  2:71673/1‑149 (MQ=255)
                        gcGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTAGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAAtt                                                                                                                           <  1:182709/149‑1 (MQ=255)
                           gTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCg                                                                                                                        >  5:380986/1‑149 (MQ=255)
                                tGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCAcc                                                                                                                   >  4:167222/1‑149 (MQ=255)
                                 gCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCt                                                                                                                  >  3:142261/1‑149 (MQ=255)
                                     aTTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTcggc                                                                                                              <  4:97976/149‑1 (MQ=255)
                                              gCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGa                                                                                                       <  4:142261/147‑1 (MQ=255)
                                              gCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGAt                                                                                                      >  8:14768/1‑148 (MQ=255)
                                                agccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTaa                                                                                                   >  4:98867/1‑149 (MQ=255)
                                                agccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTaa                                                                                                   >  8:356543/1‑149 (MQ=255)
                                                 gccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAAt                                                                                                  <  5:54184/149‑1 (MQ=255)
                                                 gccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAAt                                                                                                  >  2:343874/1‑149 (MQ=255)
                                                 gccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAAt                                                                                                  >  8:31645/1‑149 (MQ=255)
                                                  ccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTcgg                                                                                                               >  7:204289/1‑135 (MQ=255)
                                                  ccagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTcgg                                                                                                               <  8:204289/135‑1 (MQ=255)
                                                   cagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGt                                                                                                <  5:207493/149‑1 (MQ=255)
                                                   cagccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGt                                                                                                <  4:355524/149‑1 (MQ=255)
                                                      ccGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGtctc                                                                                             <  1:198688/149‑1 (MQ=255)
                                                        gggACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGATCAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGtctctc                                                                                           <  4:308590/149‑1 (MQ=255)
                                                              tttGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACg                                                                                     >  5:311118/1‑149 (MQ=255)
                                                                  ccaccaGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCg                                                                                 >  8:308407/1‑149 (MQ=255)
                                                                    accaGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGc                                                                               >  6:62334/1‑149 (MQ=255)
                                                                        gttgttTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGc                                                                           <  3:364941/149‑1 (MQ=255)
                                                                         ttgttTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGcc                                                                          >  8:322037/1‑149 (MQ=255)
                                                                            tttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgac                                                                       >  3:166490/1‑149 (MQ=255)
                                                                            tttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgac                                                                       >  1:42832/1‑149 (MQ=255)
                                                                            tttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccga                                                                        <  3:98867/148‑1 (MQ=255)
                                                                             ttGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgacc                                                                      >  4:320944/1‑149 (MQ=255)
                                                                                gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTGGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccga                                                                    >  2:232336/1‑148 (MQ=255)
                                                                                gTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGccgaccgtt                                                                   <  5:62334/149‑3 (MQ=255)
                                                                                     gaCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAg                                                               >  5:134139/1‑148 (MQ=255)
                                                                                      aCCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaa                                                             <  1:343874/149‑1 (MQ=255)
                                                                                       ccGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACTGTCGGCTAGCCGACCGAACAAGaa                                                             <  7:31645/148‑1 (MQ=255)
                                                                                         gCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTGGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaaaaa                                                          <  1:232336/149‑1 (MQ=255)
                                                                                         gCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGaaaaa                                                          <  7:308407/149‑1 (MQ=255)
                                                                                           aGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAaca                                                        <  7:135955/149‑1 (MQ=255)
                                                                                           aGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAaca                                                        <  4:119800/149‑1 (MQ=255)
                                                                                            gggAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAacac                                                       <  7:239625/149‑1 (MQ=255)
                                                                                             ggAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACt                                                      <  3:320944/149‑1 (MQ=255)
                                                                                                 cACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTagag                                                  <  6:330100/149‑1 (MQ=255)
                                                                                                 cACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTagag                                                  <  6:380986/149‑1 (MQ=255)
                                                                                                   cTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAgg                                                 >  6:69139/1‑148 (MQ=255)
                                                                                                          gCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGcactcac                                          >  5:189632/1‑148 (MQ=255)
                                                                                                                   cAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGTGTGAATTGCGTCACCTCGGCCTGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATa                                 <  1:71673/148‑1 (MQ=255)
                                                                                                                    aGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACt                               <  2:42832/149‑1 (MQ=255)
                                                                                                                     gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACt                               <  4:272929/148‑1 (MQ=255)
                                                                                                                     gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACt                               >  3:272929/1‑148 (MQ=255)
                                                                                                                     gAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTg                              >  7:102657/1‑149 (MQ=255)
                                                                                                                         aGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCAc                          >  1:352632/1‑149 (MQ=255)
                                                                                                                          ggCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCAcc                         <  7:322037/149‑1 (MQ=255)
                                                                                                                             aaaaCCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCg                      >  1:392407/1‑149 (MQ=255)
                                                                                                                                aCCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCgctg                   <  5:69139/149‑1 (MQ=255)
                                                                                                                                 cacgacacTGTGCGTTTATCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCgctgc                  >  2:148152/8‑149 (MQ=255)
                                                                                                                                       acTGTGCGTTTATCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTa                                                     >  6:327912/2‑108 (MQ=255)
                                                                                                                                       acTGTGCGTTTATCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTa                                                     <  5:327912/107‑1 (MQ=255)
                                                                                                                                        cTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGcc           <  4:90502/149‑1 (MQ=255)
                                                                                                                                        cTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGc            <  2:392407/148‑1 (MQ=255)
                                                                                                                                              gATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACGTCCGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCCCACAAAAAAAACACAAGAGGCACTCACTCATGAATACTCTCACCTCGCTGCTGCGGCCGCcccc     >  6:77878/1‑149 (MQ=255)
                                                                                                                                              gATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCAccc     >  3:386630/1‑149 (MQ=255)
                                                                                                                                                tAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGc   <  1:99145/149‑1 (MQ=255)
                                                                                                                                                  aaCGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCa  <  3:167222/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTGACCCGCAGTGAACTGGAGCAGCGGTTCTTGCCGATTTTGCTGGCAGCCAGCCGGGACCTTTGCCACCAGTTGTTTGGTTGAGACCGCAGGGAGCACTTTGCGGCCCCGGATCAGAACAGGCAAAACCGTTTCCTGTGCGATAAACGCACAGTGTCGTGCCGGGTGAATTGCGTCACCTCGGCCGGCTGATTAATGTCTCTCGCAACGGTCGGCTAGCCGACCGAACAAGAAAAACACTAGAGGCACTCACTCATGAATACTGTCACCTCGCTGCTGCGGCCGCACCCGCA  >  NC_002947/1567112‑1567405

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