Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,276 G→T 6.3% G391V (GGC→GTC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,2760GT 6.3% 253.2 / 5.1 79G391V (GGC→GTC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base T (2/3);  total (38/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 = 6.74e-01

GCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGCC  >  NC_002947/5962129‑5962423
                                                                                                                                                    |                                                                                                                                                   
gCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGc                                                                                                                                                    >  1:144108/1‑149 (MQ=255)
   agaCTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGAc                                                                                                                                                 <  1:73840/149‑1 (MQ=255)
           gACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGtt                                                                                                                                         >  8:107927/1‑149 (MQ=255)
             cTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCTAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAggatttca                                                                                                                                       <  5:15661/149‑7 (MQ=255)
                ccTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                     <  6:122077/148‑1 (MQ=255)
                ccTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                     <  1:158018/148‑1 (MQ=255)
                 cTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACAGGGCGt                                                                                                                                          <  4:19548/142‑1 (MQ=255)
                 cTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACAGGGCGt                                                                                                                                          >  3:19548/1‑142 (MQ=255)
                       aaGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccg                                                                                                                             >  1:236703/1‑149 (MQ=255)
                         ggACCGCGTCCGCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgcc                                                                                                                           <  4:115358/149‑1 (MQ=255)
                              gcgTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgccgc                                                                                                                         >  8:114936/1‑146 (MQ=255)
                                   cTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGt                                                                                                                 >  3:199485/1‑149 (MQ=255)
                                            gCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAAc                                                                                                         <  5:43544/148‑1 (MQ=255)
                                              cGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTg                                                                                                       >  4:286272/1‑148 (MQ=255)
                                              cGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGt                                                                                                      <  4:235888/149‑1 (MQ=255)
                                               gAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGc                                                                                                              <  8:262798/140‑1 (MQ=255)
                                               gAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGc                                                                                                              >  7:262798/1‑140 (MQ=255)
                                                aaGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTAt                                                                                                    >  6:114829/1‑149 (MQ=255)
                                                    tGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCACCTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                <  7:114936/149‑1 (MQ=255)
                                                    tGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGAc                                                                                                                                                 <  8:205403/100‑1 (MQ=255)
                                                    tGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGAc                                                                                                                                                 >  7:205403/1‑100 (MQ=255)
                                                    tGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                <  2:144108/149‑1 (MQ=255)
                                                      aCCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGc                                                                                              <  8:108161/149‑1 (MQ=255)
                                                       ccTGCTGGCGGTAGAGTGCGCGGGCGCCTATGGTTTTGTGATGAGCTCCAACTACTACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGc                                                                                              <  7:287150/148‑1 (MQ=255)
                                                       ccTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGc                                                                                              <  7:107927/148‑1 (MQ=255)
                                                        ctgctgGCCGTAGAGTCCGCGGCGCCTCATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCg                                                                                             >  1:190237/1‑148 (MQ=255)
                                                              gCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCct                                                                                      <  6:170494/149‑1 (MQ=255)
                                                                cGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCctgc                                                                                    >  6:285877/1‑149 (MQ=255)
                                                                   agagTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCtgctgc                                                                                 >  5:165594/1‑149 (MQ=255)
                                                                     agTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCg                                                                               <  1:242914/149‑1 (MQ=255)
                                                                     agTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCg                                                                               <  5:285877/149‑1 (MQ=255)
                                                                       tCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                             <  8:203401/149‑1 (MQ=255)
                                                                                  cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCtgctgc                                                                                 <  3:311252/134‑1 (MQ=255)
                                                                                  cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCtgctgc                                                                                 >  4:311252/1‑134 (MQ=255)
                                                                                  cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTaa                                                                         >  2:339953/1‑142 (MQ=255)
                                                                                  cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTaa                                                                         <  1:339953/142‑1 (MQ=255)
                                                                                  cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtg                                                                  >  8:40452/1‑149 (MQ=255)
                                                                                                aTGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTa                                                     <  4:112384/149‑1 (MQ=255)
                                                                                                     cTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctgct                                                             >  5:86815/1‑136 (MQ=255)
                                                                                                     cTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctgct                                                             <  6:86815/136‑1 (MQ=255)
                                                                                                       ccAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAg                                                <  6:29632/147‑1 (MQ=255)
                                                                                                       ccAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGaa                                              >  5:29632/1‑148 (MQ=255)
                                                                                                       ccAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGAt                                              >  6:47045/1‑149 (MQ=255)
                                                                                                        cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGttt                                                       >  5:140592/1‑139 (MQ=255)
                                                                                                        cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGttt                                                       <  6:140592/139‑1 (MQ=255)
                                                                                                        cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATg                                             >  3:243027/1‑149 (MQ=255)
                                                                                                        cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATg                                             >  4:4238/1‑149 (MQ=255)
                                                                                                            tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtg                                         >  5:148043/1‑149 (MQ=255)
                                                                                                             acaacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgg                                        >  7:103591/1‑149 (MQ=255)
                                                                                                             acaacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgg                                        >  3:158511/1‑149 (MQ=255)
                                                                                                               aacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggc                                      <  4:243027/149‑1 (MQ=255)
                                                                                                                 cacCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctg                                    <  3:4238/149‑1 (MQ=255)
                                                                                                                  acCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                   >  1:114452/1‑149 (MQ=255)
                                                                                                                  acCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                   <  6:165594/149‑1 (MQ=255)
                                                                                                                  acCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACAGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                   >  4:222553/1‑149 (MQ=255)
                                                                                                                    ccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGCCCGTCGCCGCACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTACGCCCATGCGGCTGCGGTTTACaaa                                                 >  8:119116/1‑133 (MQ=255)
                                                                                                                      gcgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCaa                               <  3:222553/149‑1 (MQ=255)
                                                                                                                       cgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAAc                              <  6:148043/149‑1 (MQ=255)
                                                                                                                           cGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctg                                    <  4:30730/139‑1 (MQ=255)
                                                                                                                           cGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctg                                    >  3:30730/1‑139 (MQ=255)
                                                                                                                            gTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACtt                         >  4:55559/1‑149 (MQ=255)
                                                                                                                                 gcTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATgg                    >  3:215413/1‑149 (MQ=255)
                                                                                                                                   tGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTc                  <  5:47045/149‑1 (MQ=255)
                                                                                                                                   tGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTAGCCGAACTGTATGCTCGCGAAAGCCTGCTGCCGGAGTAAGCCGATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGCCTTCATGGTc                  >  3:314468/1‑149 (MQ=255)
                                                                                                                                      aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAt                                          >  6:186018/1‑122 (MQ=255)
                                                                                                                                      aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAt                                          <  5:186018/122‑1 (MQ=255)
                                                                                                                                      aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTg               >  3:302351/1‑149 (MQ=255)
                                                                                                                                       cgccTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGtttt                                                      >  3:93407/3‑109 (MQ=255)
                                                                                                                                       cgccTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGtttt                                                      <  4:93407/107‑1 (MQ=255)
                                                                                                                                       aTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTgg              <  5:114829/149‑1 (MQ=255)
                                                                                                                                          cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCAGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctgctgc                                                           <  6:202458/101‑1 (MQ=255)
                                                                                                                                          cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCAGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctgctgc                                                           >  5:202458/1‑101 (MQ=255)
                                                                                                                                          cTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGAcc           <  3:286272/149‑1 (MQ=255)
                                                                                                                                            ggTCGACGGCGACCAGGCATTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTg         >  8:234709/1‑149 (MQ=255)
                                                                                                                                             gTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTg         <  4:158511/148‑1 (MQ=255)
                                                                                                                                             gTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTg         <  4:199485/148‑1 (MQ=255)
                                                                                                                                                gACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGt       >  6:106988/1‑147 (MQ=255)
                                                                                                                                                gACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTca     <  8:103591/149‑1 (MQ=255)
                                                                                                                                                gACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTca     <  2:114452/149‑1 (MQ=255)
                                                                                                                                                   ggCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTcagcc  <  5:106988/149‑1 (MQ=255)
                                                                                                                                                   ggCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTcagc   >  3:180049/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGCC  >  NC_002947/5962129‑5962423

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: