Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,730,161 G→C 5.1% L309L (CTC→CTG dapE ← succinyl‑diaminopimelate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,730,1610GC 5.1% 252.4 / 3.5 79L309L (CTC→CTGdapEsuccinyl‑diaminopimelate desuccinylase
Reads supporting (aligned to +/- strand):  ref base G (36/39);  new base C (2/2);  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 = 7.38e-01

ACCTGGTGGATGGTGGCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGTAGAGAAGCG  >  NC_002947/1730015‑1730295
                                                                                                                                                  |                                                                                                                                       
aCCTGGTGGATGGTGGCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGag                                                                                                                                       <  2:181356/149‑1 (MQ=255)
 ccTGGTGGATGGTGGCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGagc                                                                                                                                      <  3:366846/149‑1 (MQ=255)
 ccTGGTGGATGGTGGCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGag                                                                                                                                       >  7:327280/1‑148 (MQ=255)
        gATGGTGGCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGCAGTGCCGCCGCTGGTGGAGGGCTGGGTCTCGCGACCGGTGACGCCCTTGATGCTCGACGCCGTGGCGGCCagcatcaagg                                                                                                                               <  2:391061/149‑7 (MQ=255)
           ggtggCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCt                                                                                                                             <  2:229471/148‑1 (MQ=255)
                  ttGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGt                                                                                                                     <  6:112831/149‑1 (MQ=255)
                           ccGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACgg                                                                                                            <  2:297604/149‑1 (MQ=255)
                            cGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGc                                                                                                           >  1:397098/1‑149 (MQ=255)
                              aGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAg                                                                                                         >  3:7196/1‑149 (MQ=255)
                                ttCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGCGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGc                                                                                                       >  6:396665/1‑149 (MQ=255)
                                ttCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGc                                                                                                       >  7:109426/1‑149 (MQ=255)
                                ttCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGc                                                                                                       >  7:2746/1‑149 (MQ=255)
                                ttCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCAGCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGc                                                                                                       >  6:390919/1‑149 (MQ=255)
                                      caccTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAg                                                                                                                  >  7:150941/1‑132 (MQ=255)
                                      caccTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAg                                                                                                                  <  8:150941/132‑1 (MQ=255)
                                        ccTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGt                                                                                               <  8:327280/149‑1 (MQ=255)
                                        ccTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGt                                                                                               <  3:5255/149‑1 (MQ=255)
                                        ccTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGt                                                                                               <  5:294019/149‑1 (MQ=255)
                                         cTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTg                                                                                              >  8:244085/1‑149 (MQ=255)
                                         cTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTg                                                                                              <  8:51861/149‑1 (MQ=255)
                                            aGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGcc                                                                                            <  8:109426/148‑1 (MQ=255)
                                               gCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTACGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGt                                                                                        >  1:310682/1‑149 (MQ=255)
                                                   aTGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTcgac                                                                                    >  5:79929/1‑149 (MQ=255)
                                                   aTGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTcgac                                                                                    >  3:341792/1‑149 (MQ=255)
                                                       ttGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTcgaccgac                                                                                >  3:175682/1‑149 (MQ=255)
                                                          cGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAcca                                                                              <  6:356748/148‑1 (MQ=255)
                                                           gATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAccagt                                                                            >  3:203550/1‑149 (MQ=255)
                                                           gATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAccagt                                                                            >  6:221412/1‑149 (MQ=255)
                                                           gATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAccagt                                                                            >  1:311839/1‑149 (MQ=255)
                                                            aTGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAccagtc                                                                           <  7:21224/149‑1 (MQ=255)
                                                            aTGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAccagtc                                                                           <  4:7196/149‑1 (MQ=255)
                                                            aTGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGAccagtc                                                                           <  7:252596/149‑1 (MQ=255)
                                                                 gCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGtt                                                                      >  3:360734/1‑149 (MQ=255)
                                                                       aTCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCt                                                                >  1:362825/1‑149 (MQ=255)
                                                                        tCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCtt                                                               <  4:175682/149‑1 (MQ=255)
                                                                        tCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCtt                                                               <  4:341792/149‑1 (MQ=255)
                                                                          ggAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGt                                                             <  5:221412/149‑1 (MQ=255)
                                                                           gAAGTGCGGCCGCTGGTGGAGTGCTGTGTATCGGGACCGGTGACGCGCTTGATGCTGGCGGCCACGGCGTCGAGCAGCTCGCGTGGTTCTGTCAGGAACGGGAGGCCCGAGAGTGGCCAGTCGACCGAGGAGTCCAGTTc                                                                     <  5:94523/140‑1 (MQ=255)
                                                                           gAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTc                                                                     >  8:175875/1‑140 (MQ=255)
                                                                           gAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTc                                                                     <  7:175875/140‑1 (MQ=255)
                                                                           gAAGTGCCGCCGCTAGTGGAGGGCTGGGTATCGCGACCGCTGACGCCCTAGATGCTGGCCCCCACGGCGTCGAGCAGCACGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCCAccagtcaacatc                                                                     >  6:94523/1‑134 (MQ=255)
                                                                             aGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCga                                                          <  4:6942/149‑1 (MQ=255)
                                                                             aGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCga                                                          <  2:397098/149‑1 (MQ=255)
                                                                                gccgccGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCAt                                                                   <  1:376602/137‑1 (MQ=255)
                                                                                gccgccGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCAt                                                                   >  2:376602/1‑137 (MQ=255)
                                                                                 ccgccgCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAAt                                                      >  6:268713/1‑149 (MQ=255)
                                                                                   gccgcTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCg                                                                                                                       <  2:65624/82‑1 (MQ=255)
                                                                                   gccgcTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCg                                                                                                                       >  1:65624/1‑82 (MQ=255)
                                                                                   gccgcTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCg                                                    >  8:16108/1‑149 (MQ=255)
                                                                                      gcTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCg                                                 >  4:19299/1‑149 (MQ=255)
                                                                                           tgtaGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACt                                            <  2:310682/146‑1 (MQ=255)
                                                                                             gAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTc                                          <  2:362825/148‑1 (MQ=255)
                                                                                             gAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTACCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCg                                         >  7:89543/1‑149 (MQ=255)
                                                                                              aGGGCTGGGTCGCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGGCGAGCAGCTCGCCTGGTTCGGGCAGGAACGGCAGGCCCGACAGTGCCCAGGCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCgg                                        <  5:396665/149‑1 (MQ=255)
                                                                                              aGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCgg                                        <  5:390919/149‑1 (MQ=255)
                                                                                                     ggTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTg                                   >  1:246982/1‑147 (MQ=255)
                                                                                                      gTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGcc                                                  <  2:53764/132‑1 (MQ=255)
                                                                                                      gTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGcc                                                  >  1:53764/1‑132 (MQ=255)
                                                                                                      gTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAg                                <  3:19299/149‑1 (MQ=255)
                                                                                                              aCCGGTGACGCCGTTGTTGCTGGCCGCCACGTCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGTCCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGa                        <  8:89543/149‑1 (MQ=255)
                                                                                                               ccGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGAc                       <  2:246982/149‑1 (MQ=255)
                                                                                                                 ggTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACgg                     <  6:79929/149‑1 (MQ=255)
                                                                                                                  gTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGt                    <  7:244085/149‑1 (MQ=255)
                                                                                                                  gTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGt                    >  8:369470/1‑149 (MQ=255)
                                                                                                                     aCGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGa                 <  4:360734/149‑1 (MQ=255)
                                                                                                                       gCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACt               >  8:321301/1‑149 (MQ=255)
                                                                                                                        cccTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTc              <  5:268713/149‑1 (MQ=255)
                                                                                                                        cccTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTc              >  4:293350/1‑149 (MQ=255)
                                                                                                                         ccTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCg             >  5:377335/1‑149 (MQ=255)
                                                                                                                          cTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCCGGGCCTGCAGGCCTTCGACGGTGGACTCg             <  1:31193/149‑1 (MQ=255)
                                                                                                                           ttGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCgg            >  2:6342/1‑148 (MQ=255)
                                                                                                                           ttGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGt           >  5:248454/1‑149 (MQ=255)
                                                                                                                           gagctGCTCGACGCCGTGGCGGCCAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCgg            >  6:147319/5‑148 (MQ=255)
                                                                                                                            agctGCTCGACGCCGTGGCGGCCAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGTa          <  5:147319/146‑1 (MQ=255)
                                                                                                                             gATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGTa          <  1:6342/148‑1 (MQ=255)
                                                                                                                              ctGCTCGACGCCGTGGCGGCCAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGTa          >  6:165051/2‑147 (MQ=255)
                                                                                                                                    gccgccACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGt                    <  2:285087/131‑1 (MQ=255)
                                                                                                                                    gccgccACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGt                    >  1:285087/1‑131 (MQ=255)
                                                                                                                                    gccgccACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGTAGAGAAGCg  >  2:307333/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
ACCTGGTGGATGGTGGCGTTGACCGGGCCGAGTTCGACCACCTGAGTGCCCATGGTTGCGATGAAGCGGCCATCGGAAGTGCCGCCGCTGGTGGAGGGCTGGGTATCGCGACCGGTGACGCCCTTGATGCTGGCCGCCACGGCGTCGAGCAGCTCGCCTGGTTCGGTCAGGAACGGCAGGCCCGACAGTGCCCAGTCGACCGACCAGTCCAGTTCATGCTTGTCGAGAATCGCCGACACTCGGGCCTGCAGGCCTTCGACGGTGGACTCGGTAGAGAAGCG  >  NC_002947/1730015‑1730295

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: