Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,924,777 G→A 7.0% H13Y (CAC→TAC)  peaB ← quinohemoprotein amine dehydrogenase modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,924,7770GA 7.0% 255.6 / 6.3 85H13Y (CAC→TAC) peaBquinohemoprotein amine dehydrogenase modification protein
Reads supporting (aligned to +/- strand):  ref base G (41/38);  new base A (3/3);  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 = 7.10e-01

CTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTG  >  NC_002947/3924628‑3924924
                                                                                                                                                     |                                                                                                                                                   
tgctcggCCGAGGAGCTGATGCAGCGCCTTGGCCAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTg                                                                                                                                                     >  3:246567/6‑149 (MQ=37)
   ctgGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCTCCTGTCACCGCGTCAGCGGTTAACAGCGAACTGCTGGGGATATGCAACAGCATGCGGTCGGCATCAACCTGTACTTCGTg                                                                                                                                                     <  4:268544/147‑1 (MQ=255)
   ctgGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGc                                                                                                                                                    <  1:257885/148‑1 (MQ=255)
   ctgGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCa                                                                                                                                                   <  8:4846/149‑1 (MQ=255)
       ccAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGtt                                                                                                                                               >  7:119409/1‑149 (MQ=255)
        cAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGGATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTg                                                                                                                                              >  5:325632/1‑149 (MQ=255)
        cAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTg                                                                                                                                              >  8:23222/1‑149 (MQ=255)
         aaGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGc                                                                                                                                             >  3:19378/1‑149 (MQ=255)
                      aGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTc                                                                                                                                         <  7:147086/140‑1 (MQ=255)
                      aGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTc                                                                                                                                         >  8:147086/1‑140 (MQ=255)
                       gCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGtt                                                                                                                                >  6:3030/1‑148 (MQ=255)
                        ctcctcGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCa                                                                                                                              <  1:132405/149‑1 (MQ=255)
                           ctcGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCag                                                                                                                             >  8:230969/1‑147 (MQ=255)
                           ctcGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCag                                                                                                                             <  8:119409/147‑1 (MQ=255)
                           ctcGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCag                                                                                                                             <  3:278312/147‑1 (MQ=255)
                              ggCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAgtagta                                                                                                                        >  8:18533/1‑149 (MQ=255)
                                 cGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGc                                                                                                                     <  3:146729/149‑1 (MQ=255)
                                  gAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGcc                                                                                                                    <  7:230969/149‑1 (MQ=255)
                                            gCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGAcc                                                                                                          >  4:245484/1‑149 (MQ=255)
                                                gcgcAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTccc                                                                                                      <  3:294511/149‑1 (MQ=255)
                                                gcgcAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTcc                                                                                                       <  6:58796/148‑1 (MQ=255)
                                                  gcAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGAcc                                                                                                          <  1:34513/143‑1 (MQ=255)
                                                  gcAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGAcc                                                                                                          >  2:34513/1‑143 (MQ=255)
                                                  gcAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGa                                                                                                    >  7:105789/1‑149 (MQ=255)
                                                  gcAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGa                                                                                                    <  7:18533/149‑1 (MQ=255)
                                                     aGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATaa                                                                                                 <  6:35476/149‑1 (MQ=255)
                                                      ggTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAAt                                                                                                <  2:314115/149‑1 (MQ=255)
                                                      ggTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAAt                                                                                                <  7:60669/149‑1 (MQ=255)
                                                             atgatgCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAgg                                                                                                                                  <  7:65777/108‑1 (MQ=255)
                                                             atgatgCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAgg                                                                                                                                  >  8:65777/1‑108 (MQ=255)
                                                                 tgCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTg                                                                                      >  8:206541/1‑148 (MQ=255)
                                                                 tgCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTg                                                                                      <  7:206541/148‑1 (MQ=255)
                                                                   ccgccgGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGa                                                                                    >  5:134274/1‑148 (MQ=255)
                                                                            aCCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTg                                                                           >  4:271999/1‑148 (MQ=255)
                                                                            aCCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTg                                                                           <  5:250654/148‑1 (MQ=255)
                                                                            aCCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTg                                                                           >  6:199011/1‑148 (MQ=255)
                                                                              cgcgTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGc                                                                                                              >  4:315501/1‑111 (MQ=255)
                                                                              cgcgTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGc                                                                                                              <  3:315501/111‑1 (MQ=255)
                                                                              cgcgTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAg                                                                        >  6:316448/1‑149 (MQ=255)
                                                                                             cAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGtt                                                                            <  1:79238/130‑1 (MQ=255)
                                                                                             cAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGtt                                                                            >  2:79238/1‑130 (MQ=255)
                                                                                               gCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAg                                                        >  7:114865/1‑148 (MQ=255)
                                                                                               gCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAg                                                        <  8:114865/148‑1 (MQ=255)
                                                                                                  aaCTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAgcggc                                                    <  3:245484/149‑1 (MQ=255)
                                                                                                  aaCTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAgcggc                                                    >  3:153535/1‑149 (MQ=255)
                                                                                                     tgctgGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCggg                                                 >  7:244381/1‑149 (MQ=255)
                                                                                                     tgctgGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCggg                                                 >  4:204492/1‑149 (MQ=255)
                                                                                                       ctgGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCggg                                                 <  7:23222/147‑1 (MQ=255)
                                                                                                          gggATATGGAACAGCATGCGGTCGACATAAACCTTTAATTCGTGCAGGTTGCATTAAACCAGGTTCAGTAATACGCCCATGGCGACCTCCCGATAATCGAACCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGgttgtt                                            >  7:50171/1‑149 (MQ=255)
                                                                                                                tGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAgc                                       <  4:19378/148‑1 (MQ=255)
                                                                                                                  gAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAgcgc                                     >  2:319660/1‑148 (MQ=255)
                                                                                                                  gAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAgcgc                                     <  1:319660/148‑1 (MQ=255)
                                                                                                                  gAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCt                                    <  6:134274/149‑1 (MQ=255)
                                                                                                                   aaCAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTg                                   >  4:107248/1‑149 (MQ=255)
                                                                                                                     cccgaTGCGGTCGTCATCAACCTGTACTTCGTGCAGGTTGCGTGCGACCAGGTTCAGTAGTAGGCCCCTGGCGACCTCGCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGt                                  <  7:125451/144‑1 (MQ=255)
                                                                                                                     cAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTa                                 >  6:244996/1‑149 (MQ=255)
                                                                                                                     cAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTa                                 <  6:325632/149‑1 (MQ=255)
                                                                                                                        cATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACgg                              >  7:134133/1‑149 (MQ=255)
                                                                                                                         aTGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGt                             <  5:3030/149‑1 (MQ=255)
                                                                                                                            cGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGa                           <  4:246567/148‑1 (MQ=255)
                                                                                                                             ggTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAgc                                       >  8:242846/1‑135 (MQ=255)
                                                                                                                             ggTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAgc                                       <  7:242846/135‑1 (MQ=255)
                                                                                                                                 ggCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCt                                    <  6:73808/134‑1 (MQ=255)
                                                                                                                                 ggCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCt                                    >  5:73808/1‑134 (MQ=255)
                                                                                                                                 ggCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCa                     >  4:180452/1‑149 (MQ=255)
                                                                                                                                 ggCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCa                     <  5:244996/149‑1 (MQ=255)
                                                                                                                                 gaacgcAACCTGCACGAAGTACAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAgcggc                                                    <  6:282783/113‑1 (MQ=255)
                                                                                                                                 gaacgcAACCTGCACGAAGTACAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAgcggc                                                    >  5:282783/6‑118 (MQ=255)
                                                                                                                                  aacgcAACCTGCACGAAGTACAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGgttgt                                             <  1:175299/120‑1 (MQ=255)
                                                                                                                                  aacgcAACCTGCACGAAGTACAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGgttgt                                             >  2:175299/5‑124 (MQ=255)
                                                                                                                                  aacgcAACCTGCACGAAGTACAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGgttgt                                             >  8:187352/5‑124 (MQ=255)
                                                                                                                                  aacgcAACCTGCACGAAGTACAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGgttgt                                             <  7:187352/120‑1 (MQ=255)
                                                                                                                                   cATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTACAGCGCTGTACGGTGACGATCAtg                   >  5:279650/1‑149 (MQ=255)
                                                                                                                                    aTCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCAtgt                  >  6:97084/1‑149 (MQ=255)
                                                                                                                                    aTCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCAtgt                  >  7:143037/1‑149 (MQ=255)
                                                                                                                                     tCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCAtgtg                 >  8:257276/1‑149 (MQ=255)
                                                                                                                                        aCCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATTATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGccc              >  3:275596/1‑149 (MQ=255)
                                                                                                                                        aCCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTcc                                          <  7:236962/121‑1 (MQ=255)
                                                                                                                                        aCCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTcc                                          >  8:236962/1‑121 (MQ=255)
                                                                                                                                          cTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGAc                          <  7:231590/135‑1 (MQ=255)
                                                                                                                                          cTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGAc                          >  8:231590/1‑135 (MQ=255)
                                                                                                                                              aCTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACACCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGAc                          <  2:223640/131‑1 (MQ=255)
                                                                                                                                              aCTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACACCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGAc                          >  1:223640/1‑131 (MQ=255)
                                                                                                                                                 tCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGt     >  7:39448/1‑149 (MQ=255)
                                                                                                                                                  cGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGt     <  6:227675/148‑1 (MQ=255)
                                                                                                                                                    tGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTg  <  7:257276/149‑1 (MQ=255)
                                                                                                                                                    tGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACAGTGACGATCATGTGCCCTTCGCCGGTCTg  >  2:30573/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCAGCGGTGAACAGCGAACTGCTGGGGATATGGAACAGCATGCGGTCGGCATCAACCTGTACTTCGTGCAGGTTGCGTTCGACCAGGTTCAGTAGTACGCCCATGGCGACCTCCCGATAATCGATCCGTTGGAAAACCGTTGCAGCCTGCCGTTACGGCAGCGGCGGGTTGTTCCAGCGCTGTACGGTGACGATCATGTGCCCTTCGCCGGTCTG  >  NC_002947/3924628‑3924924

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: