Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,566 T→A 10.9% L128Q (CTG→CAG)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,5660TA10.9% 169.6 / 10.4 64L128Q (CTG→CAG) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base A (3/4);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGG  >  NC_002947/3528423‑3528705
                                                                                                                                               |                                                                                                                                            
cGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCctgctg                                                                                                                                          <  6:8103/148‑1 (MQ=255)
  aCGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGcc                                                                                                                                       <  7:127141/149‑1 (MQ=255)
        ccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCCACAAGGCCAGCAGGCCTt‑ggtg                                                                                                                                <  6:49566/149‑4 (MQ=255)
              tGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTt‑ggtg                                                                                                                                <  1:317234/143‑4 (MQ=255)
              tGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTt‑ggtg                                                                                                                                >  2:317234/1‑140 (MQ=255)
                gTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGc                                                                                                                        >  5:362620/1‑149 (MQ=255)
                   ctgctgGTCGAAAACCTGCCGGTTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGa                                                                                                                     <  5:83109/149‑1 (MQ=255)
                       tgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGt                                                                                                                 <  3:11499/149‑1 (MQ=255)
                        ggTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTg                                                                                                                <  3:197896/149‑1 (MQ=255)
                          tCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCTCACGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGc                                                                                                              <  5:312984/149‑1 (MQ=255)
                          tCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGc                                                                                                              >  1:138548/1‑149 (MQ=255)
                          tCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGc                                                                                                              <  2:138548/149‑1 (MQ=255)
                          tCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCACCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGc                                                                                                              >  2:18694/1‑149 (MQ=255)
                            gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGa                                                                                                            >  8:51101/1‑149 (MQ=255)
                            gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGa                                                                                                            >  4:214864/1‑149 (MQ=255)
                            gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGa                                                                                                            >  3:90737/1‑149 (MQ=255)
                             aaaaCCCGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGa                                                                                                            >  6:155265/1‑148 (MQ=255)
                                  atGACGGGTTACAAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAATGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGACTGCTGGACTT‑GTGGGTGGACGATGACAGTTGCGATACTGc                                                                                                      >  2:312006/2‑149 (MQ=255)
                                     aCGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTAc                                                                                                   <  3:214864/149‑1 (MQ=255)
                                             ccAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGt                                                                                           >  2:284493/1‑149 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGc                                                                                                      <  7:15557/137‑1 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGc                                                                                                      >  8:15557/1‑137 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTg                                                                                          >  7:209144/1‑149 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTg                                                                                          >  7:220784/1‑149 (MQ=255)
                                                  aCCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCt                                                                                       <  8:220784/148‑1 (MQ=255)
                                                       tgCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAg                                                                                 >  6:95425/1‑149 (MQ=255)
                                                            gcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTTGCGATACTGCTACCGGCGCGTGGCTGCGCAgtgtgt                                                                            <  7:7657/149‑1 (MQ=255)
                                                            gcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAgtgtgt                                                                            >  8:231965/1‑149 (MQ=255)
                                                                aTCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGGCTGCTGAACTT‑GTGGAAGAGCGCTGACAGTGGCGATACTGCAACCGGCGCGTGGCTGCGCAGTGTGTTCGc                                                                        >  1:115647/1‑149 (MQ=255)
                                                                   aCCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa                                                                     <  8:209144/149‑1 (MQ=255)
                                                                      cGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgcg                                                                  >  8:9806/1‑149 (MQ=255)
                                                                         gCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                               >  3:351779/1‑149 (MQ=255)
                                                                          cccGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCtg                                                              <  4:351779/149‑1 (MQ=255)
                                                                            cggcggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCtgtg                                                            >  3:306907/1‑149 (MQ=255)
                                                                                ggcggGCAGAGAAGGTGGCGTAACAGCTGTTTCGTTATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGACGATGACAGTGGCGATACTGCTACAGGCGAGTGGCTGCTCACTGTGTTCGCTGAGCGCCttagatc                                                         >  4:322220/1‑142 (MQ=255)
                                                                                ggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                               <  3:322220/142‑1 (MQ=255)
                                                                                   gggCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGcc                                                     <  7:231965/149‑1 (MQ=255)
                                                                                     gCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGt                                                   >  5:337349/1‑149 (MQ=255)
                                                                                     gCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGt                                                   <  5:155265/149‑1 (MQ=255)
                                                                                      cAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTg                                                  <  1:284493/149‑1 (MQ=255)
                                                                                             gggggCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGc                                           <  7:51101/149‑1 (MQ=255)
                                                                                                    tACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCgggg                                    <  4:90737/149‑1 (MQ=255)
                                                                                                      ccAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGGGCGCCTGTGGCTGGCCGTGCACCTGCATCgggg                                    >  6:252379/1‑147 (MQ=255)
                                                                                                       cAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGGGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAg                                 <  5:252379/149‑1 (MQ=255)
                                                                                                                   gatgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGTGCAGTGCCGATGCCGt                     >  4:211576/1‑149 (MQ=255)
                                                                                                                        cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                               >  5:275718/1‑102 (MQ=255)
                                                                                                                        cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                               <  6:275718/102‑1 (MQ=255)
                                                                                                                        cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGcc                <  6:362620/149‑1 (MQ=255)
                                                                                                                         ttGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCt               >  6:238425/1‑149 (MQ=255)
                                                                                                                           gctgctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGa                           <  6:48831/135‑1 (MQ=255)
                                                                                                                           gctgctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGa                           >  5:48831/1‑135 (MQ=255)
                                                                                                                           gctgctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCt               <  1:312006/147‑1 (MQ=255)
                                                                                                                              gctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGt                     <  1:263207/138‑1 (MQ=255)
                                                                                                                              gctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGt                     >  2:263207/1‑138 (MQ=255)
                                                                                                                                  caccacCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGtt      >  4:145246/1‑149 (MQ=255)
                                                                                                                                   acccacaAGGCCAGCAGGCCTT‑GGTGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGAc                             <  3:283074/120‑1 (MQ=255)
                                                                                                                                   acccacaAGGCCAGCAGGCCTT‑GGTGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGAc                             >  4:283074/6‑125 (MQ=255)
                                                                                                                                   accaccAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTg     >  2:19547/1‑149 (MQ=255)
                                                                                                                                      accaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg  >  8:137066/1‑149 (MQ=255)
                                                                                                                                      accaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg  >  8:196826/1‑149 (MQ=255)
                                                                                                                                      accaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg  <  4:306907/149‑1 (MQ=255)
                                                                                                                                       ccaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGAt                          >  6:167692/1‑124 (MQ=255)
                                                                                                                                       ccaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGAt                          <  5:167692/124‑1 (MQ=255)
                                                                                                                                       acaAGGCCAGCAGGCCTTGGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtggctg                                                        <  7:61241/94‑1 (MQ=255)
                                                                                                                                       acaAGGCCAGCAGGCCTTGGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtggctg                                                        >  8:61241/2‑95 (MQ=255)
                                                                                                                                               |                                                                                                                                            
CGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGG  >  NC_002947/3528423‑3528705

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: ‑