Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,961,631 C→G 8.1% A176G (GCC→GGC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,961,6310CG 8.1% 192.6 / 6.4 62A176G (GCC→GGC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (3/2);  total (33/29)
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.95e-01

CTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCTGGCC  >  NC_002947/5961492‑5961779
                                                                                                                                           |                                                                                                                                                     
cTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgccga                                                                                                                                              <  3:301778/149‑1 (MQ=255)
   gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGGCAAtggcgatgccgaac                                                                                                                                            <  8:86251/148‑3 (MQ=255)
   gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCACTGCGGGTCAACCCGGACGTCGATGCCGGTACACACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGc                                                                                                                                           >  1:90173/1‑149 (MQ=255)
         cTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCta                                                                                                                                     >  6:286956/1‑149 (MQ=255)
                  gtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACgcgc                                                                                                                            >  8:196525/1‑149 (MQ=255)
                   tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                           <  2:136020/149‑1 (MQ=255)
                    cgttttcgAAATTGGCAATGTTTCCCCGGTGTCGCTTCGGGTCAACCCGGACGTCGATTCCGGTACCCACCCGTACATCACCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                           <  3:266338/142‑1 (MQ=255)
                     gccgccGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgccg                                                                                                                          >  8:150923/1‑148 (MQ=255)
                       cgccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTAACCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGATGCCATCTACTTCAGCgcagca                                                                                                                        >  3:266487/1‑148 (MQ=255)
                             aaTGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                                 >  7:288537/1‑149 (MQ=255)
                                gggCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgcc                                                                                                                                                <  6:321171/115‑1 (MQ=255)
                                gggCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgcc                                                                                                                                                >  5:321171/1‑115 (MQ=255)
                                  gCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAAc                                                                                                            >  7:91145/1‑149 (MQ=255)
                                      ggTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTg                                                                                                         >  5:7391/1‑148 (MQ=255)
                                        ttGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                                      <  4:77469/149‑1 (MQ=255)
                                        atGCCCCGGTGGCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTAACAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                                      <  8:322879/148‑1 (MQ=255)
                                            cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                                  >  6:177251/1‑149 (MQ=255)
                                            cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                                  <  6:7391/149‑1 (MQ=255)
                                            cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                                  >  4:227746/1‑149 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGGTGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                                 <  7:239224/149‑1 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCt                                                                                                                                      >  6:187951/1‑112 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCt                                                                                                                                      <  5:187951/112‑1 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                                 <  8:91145/149‑1 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                                 <  3:227746/149‑1 (MQ=255)
                                               ggTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                               <  7:150923/149‑1 (MQ=255)
                                                   tCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgtg                                                                                            <  8:46568/148‑1 (MQ=255)
                                                   tCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACATGGAAGTGGTCGgtgt                                                                                           >  4:56280/1‑149 (MQ=255)
                                                                  ccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCgg                                                                            <  4:16556/149‑1 (MQ=255)
                                                                   cGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGc                                                                           <  3:56280/149‑1 (MQ=255)
                                                                             tGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑Ga                                                                <  7:196525/149‑1 (MQ=255)
                                                                               ccGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑Gac                                                               <  5:177251/148‑1 (MQ=255)
                                                                                     cccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                        >  6:164939/1‑149 (MQ=255)
                                                                                      ccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                        >  4:180065/1‑148 (MQ=255)
                                                                                       cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                        >  5:128874/1‑147 (MQ=255)
                                                                                       cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTc                                                        >  3:241627/1‑147 (MQ=255)
                                                                                        acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                                   >  8:306981/1‑120 (MQ=255)
                                                                                        acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                                   <  7:306981/120‑1 (MQ=255)
                                                                                        acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAg                                                     <  3:180065/149‑1 (MQ=255)
                                                                                        acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAg                                                     >  2:199558/1‑149 (MQ=255)
                                                                                         cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGc                                                    <  5:286956/149‑1 (MQ=255)
                                                                                         cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGc                                                    >  6:270281/1‑149 (MQ=255)
                                                                                          ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                              >  7:132340/1‑123 (MQ=255)
                                                                                          ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                              <  8:132340/123‑1 (MQ=255)
                                                                                                   tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGa                                          >  5:342464/1‑149 (MQ=255)
                                                                                                   tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGa                                          >  5:18880/1‑149 (MQ=255)
                                                                                                         ggCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCt                                    <  4:241627/149‑1 (MQ=255)
                                                                                                          gCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTcc                                              >  5:3448/1‑138 (MQ=255)
                                                                                                          gCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTcc                                              <  6:3448/138‑1 (MQ=255)
                                                                                                          gCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTc                                   >  3:110131/1‑149 (MQ=255)
                                                                                                            cTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGa                                 <  2:90173/149‑1 (MQ=255)
                                                                                                             ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGAc                                         <  5:345511/140‑1 (MQ=255)
                                                                                                             ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGACCCATCGGCTAGCAGCT‑GACCACCGTCGAGCCCTTCCTCGAc                                         >  6:345511/1‑140 (MQ=255)
                                                                                                               aaaGAAAACAAGTTCGGCATCGCCATTGCCGCTGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCg                              >  1:19299/1‑149 (MQ=255)
                                                                                                                 aGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTCGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGc                             <  1:243970/149‑1 (MQ=255)
                                                                                                                 aGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGcc                            <  6:128874/149‑1 (MQ=255)
                                                                                                                     aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgct                        >  3:256870/1‑149 (MQ=255)
                                                                                                                        aaGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggt                     <  5:270281/149‑1 (MQ=255)
                                                                                                                          gTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGAAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc                   >  1:15650/1‑149 (MQ=255)
                                                                                                                          cctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc                   >  2:265513/3‑149 (MQ=255)
                                                                                                                           ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTc                                            >  2:39221/2‑123 (MQ=255)
                                                                                                                           ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTc                                            <  1:39221/122‑1 (MQ=255)
                                                                                                                            tcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg                 >  8:334842/1‑149 (MQ=255)
                                                                                                                                 aTCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGa            <  1:265513/149‑1 (MQ=255)
                                                                                                                                    gCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCg         >  6:174590/1‑149 (MQ=255)
                                                                                                                                     ccATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGc        >  5:190892/1‑149 (MQ=255)
                                                                                                                                       aTTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCt      <  6:342464/149‑1 (MQ=255)
                                                                                                                                         tgccgatgccgaGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc                   >  5:235061/1‑134 (MQ=255)
                                                                                                                                         tgccgatgccgaGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc                   <  6:235061/134‑1 (MQ=255)
                                                                                                                                           ccgatgccgaGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCTGGcc  >  6:325213/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
CTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCT‑GACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGCCTGGCC  >  NC_002947/5961492‑5961779

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