Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,569 T→A 10.9% L129Q (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,5690TA10.9% 172.5 / 10.8 64L129Q (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 = 7.82e-01

CGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCT  >  NC_002947/3528423‑3528717
                                                                                                                                                  |                                                                                                                                                      
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‑GTGGGTGGACGATGACAGTTGCGATACT‑Gc                                                                                                                  >  2:312006/2‑149 (MQ=255)
                                     aCGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTAc                                                                                                               <  3:214864/149‑1 (MQ=255)
                                             ccAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGt                                                                                                       >  2:284493/1‑149 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑Gc                                                                                                                  <  7:15557/137‑1 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑Gc                                                                                                                  >  8:15557/1‑137 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTg                                                                                                      >  7:209144/1‑149 (MQ=255)
                                              cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTg                                                                                                      >  7:220784/1‑149 (MQ=255)
                                                  aCCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCt                                                                                                   <  8:220784/148‑1 (MQ=255)
                                                       tgCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAg                                                                                             >  6:95425/1‑149 (MQ=255)
                                                            gcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTTGCGATACT‑GCTACCGGCGCGTGGCTGCGCAgtgtgt                                                                                        <  7:7657/149‑1 (MQ=255)
                                                            gcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAgtgtgt                                                                                        >  8:231965/1‑149 (MQ=255)
                                                                aTCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGGCTGCTGAACTT‑GTGGAAGAGCGCTGACAGTGGCGATACT‑GCAACCGGCGCGTGGCTGCGCAGTGTGTTCGc                                                                                    >  1:115647/1‑149 (MQ=255)
                                                                   aCCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa                                                                                 <  8:209144/149‑1 (MQ=255)
                                                                      cGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgcg                                                                              >  8:9806/1‑149 (MQ=255)
                                                                         gCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                                           >  3:351779/1‑149 (MQ=255)
                                                                          cccGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCtg                                                                          <  4:351779/149‑1 (MQ=255)
                                                                            cggcggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCtgtg                                                                        >  3:306907/1‑149 (MQ=255)
                                                                                ggcggGCAGAGAAGGTGGCGTAACAGCTGTTTCGTTATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGACGATGACAGTGGCGATACT‑GCTACAGGCGAGTGGCTGCTCACTGTGTTCGCTGAGCGCCttagatc                                                                     >  4:322220/1‑142 (MQ=255)
                                                                                ggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                                           <  3:322220/142‑1 (MQ=255)
                                                                                   gggCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGcc                                                                 <  7:231965/149‑1 (MQ=255)
                                                                                     gCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGt                                                               <  5:155265/149‑1 (MQ=255)
                                                                                     gCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGt                                                               >  5:337349/1‑149 (MQ=255)
                                                                                      cAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTg                                                              <  1:284493/149‑1 (MQ=255)
                                                                                             gggggCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGc                                                       <  7:51101/149‑1 (MQ=255)
                                                                                                    tACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCgggg                                                <  4:90737/149‑1 (MQ=255)
                                                                                                      ccAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGGGCGCCTGTGGCTGGCCGTGCACCTGCATCgggg                                                >  6:252379/1‑147 (MQ=255)
                                                                                                       cAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGGGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAg                                             <  5:252379/149‑1 (MQ=255)
                                                                                                                   gatgaCTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGTGCAGTGCCGATGCCGt                                 >  4:211576/1‑149 (MQ=255)
                                                                                                                        cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                                           <  6:275718/102‑1 (MQ=255)
                                                                                                                        cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCt                                                                           >  5:275718/1‑102 (MQ=255)
                                                                                                                        cTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGcc                            <  6:362620/149‑1 (MQ=255)
                                                                                                                         ttGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCt                           >  6:238425/1‑149 (MQ=255)
                                                                                                                           gctgctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGa                                       >  5:48831/1‑135 (MQ=255)
                                                                                                                           gctgctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGa                                       <  6:48831/135‑1 (MQ=255)
                                                                                                                           gctgctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCt                           <  1:312006/147‑1 (MQ=255)
                                                                                                                              gctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGt                                 >  2:263207/1‑138 (MQ=255)
                                                                                                                              gctgcACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGt                                 <  1:263207/138‑1 (MQ=255)
                                                                                                                                  caccacCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGtt                  >  4:145246/1‑149 (MQ=255)
                                                                                                                                   acccacaAGGCCAGCAGGCCTT‑GGTGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGAc                                         >  4:283074/6‑125 (MQ=255)
                                                                                                                                   acccacaAGGCCAGCAGGCCTT‑GGTGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGAc                                         <  3:283074/120‑1 (MQ=255)
                                                                                                                                   accaccAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTg                 >  2:19547/1‑149 (MQ=255)
                                                                                                                                      accaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg              >  8:137066/1‑149 (MQ=255)
                                                                                                                                      accaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg              <  4:306907/149‑1 (MQ=255)
                                                                                                                                      accaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCgg              >  8:196826/1‑149 (MQ=255)
                                                                                                                                       ccaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGAt                                      >  6:167692/1‑124 (MQ=255)
                                                                                                                                       ccaAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGAt                                      <  5:167692/124‑1 (MQ=255)
                                                                                                                                       acaAGGCCAGCAGGCCTTGGTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtggctg                                                                    <  7:61241/94‑1 (MQ=255)
                                                                                                                                       acaAGGCCAGCAGGCCTTGGTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtggctg                                                                    >  8:61241/2‑95 (MQ=255)
                                                                                                                                                 ctgGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACTGGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcgg    >  1:63858/1‑149 (MQ=255)
                                                                                                                                                 ctgGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTggcgg    <  7:9806/148‑1 (MQ=255)
                                                                                                                                                  tgGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCt  <  5:127985/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CGACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTT‑GTGGGTGGCCGCTGACAGTGGCGATACT‑GCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGCACTGGCGGCT  >  NC_002947/3528423‑3528717

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