Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,273 A→C 6.3% D390A (GAC→GCC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,2730AC 6.3% 233.7 / 4.9 76D390A (GAC→GCC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base C (3/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGA  >  NC_002947/5962132‑5962412
                                                                                                                                             |                                                                                                                                           
agaCTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGAc                                                                                                                                      <  5:345005/149‑1 (MQ=255)
      gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAgga                                                                                                                                 <  5:187050/148‑2 (MQ=255)
      gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAgga                                                                                                                                 >  6:187050/1‑147 (MQ=255)
             ccTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGAcc                                                                                                                                     <  5:312777/137‑1 (MQ=255)
             ccTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGAcc                                                                                                                                     >  6:312777/1‑137 (MQ=255)
                 ggCAAGGACCGCGTCCCCAACCTGGCCGAAGGTGACCTGCTGGTCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAc                                                                                                                      <  3:215601/148‑1 (MQ=255)
                  gCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgc                                                                                                                    <  1:104034/149‑1 (MQ=255)
                      ggACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgcc                                                                                                                >  2:48483/1‑149 (MQ=255)
                       gACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgc                                                                                                                 >  8:153720/1‑147 (MQ=255)
                       gACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgc                                                                                                                 <  7:153720/147‑1 (MQ=255)
                          cgcgTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgccgc                                                                                                              >  6:64147/1‑147 (MQ=255)
                              tCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGAcc                                                                                                        >  3:358107/1‑149 (MQ=255)
                                   aaCCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTc                                                                                                     <  4:173522/147‑1 (MQ=255)
                                   aaCCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGc                                                                                                   >  7:310940/1‑149 (MQ=255)
                                         tccGAAGGTGACCTGCTGGCGGTATTGTCCGCGGGCGCCTATTGTTTTGTCATGAGCTCCATCTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGAGCGTCTCCGAAc                                                                                              <  4:384088/147‑1 (MQ=255)
                                           cGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACAAGGCGTTAGAAGTACGCCGCCGCGAGACCGTCGCCGAAc                                                                                              >  3:384088/1‑146 (MQ=255)
                                                gTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCt                                                                                      <  5:64147/149‑1 (MQ=255)
                                                 tGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                     <  8:310940/149‑1 (MQ=255)
                                                  gACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTgg                                                                                    >  5:54025/1‑149 (MQ=255)
                                                  gACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTgg                                                                                    <  1:309447/149‑1 (MQ=255)
                                                  gACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTgg                                                                                    <  3:3272/149‑1 (MQ=255)
                                                     ctgctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCg                                                                                  >  6:387514/1‑148 (MQ=255)
                                                        ctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGaaa                                                                               <  4:358107/148‑1 (MQ=255)
                                                         tgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCga                                                                                                            >  7:351867/1‑118 (MQ=255)
                                                         tgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCga                                                                                                            <  8:351867/118‑1 (MQ=255)
                                                           gCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCct                                                                           <  6:54025/149‑1 (MQ=255)
                                                               tAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCctgctg                                                                       <  3:319951/149‑1 (MQ=255)
                                                               tAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCctgctg                                                                       >  7:86819/1‑149 (MQ=255)
                                                                    tCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                  <  5:90880/149‑1 (MQ=255)
                                                                    tCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                  >  2:198548/1‑149 (MQ=255)
                                                                        cgGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTaa                                                              <  1:48483/149‑1 (MQ=255)
                                                                          ggCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCa                                                                                                                                    <  5:83462/77‑1 (MQ=255)
                                                                          ggCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCa                                                                                                                                    >  6:83462/1‑77 (MQ=255)
                                                                             gcCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTaa                                                              <  6:191640/144‑1 (MQ=255)
                                                                             gcCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTaa                                                              >  5:191640/1‑144 (MQ=255)
                                                                               cTATGGTTTTGTCATGAGATCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtg                                                       >  6:296475/1‑149 (MQ=255)
                                                                                tATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGaaa                                                                               >  4:291123/1‑124 (MQ=255)
                                                                                tATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGaaa                                                                               <  3:291123/124‑1 (MQ=255)
                                                                                 aTGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgct                                                     >  1:20593/1‑149 (MQ=255)
                                                                                  tGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctg                                                    <  8:86819/149‑1 (MQ=255)
                                                                                          tCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                  <  4:162969/127‑1 (MQ=255)
                                                                                          tCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                  >  3:162969/1‑127 (MQ=255)
                                                                                          tCATGAGCTCCAACTACAACACCCGCGGCCGATGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACAGTCGCCGAACTGTATGCTGGCGAAAGCCTTCTGCCGGAGTAAGACCAAGCAGCTGCGttt                                            >  4:96828/1‑149 (MQ=255)
                                                                                          gcATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGttt                                            <  6:265591/148‑1 (MQ=255)
                                                                                            aTGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTa                                          >  3:285847/1‑149 (MQ=255)
                                                                                            aTGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTa                                          >  3:226900/1‑149 (MQ=255)
                                                                                               atcTCCAACTACAACACCCGCGGCCGTTGCGCGGAAATCCTGGTCGACGGCGACCAGGCGTTCGTAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTAcc                                        <  1:181214/146‑1 (MQ=255)
                                                                                                gCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTTCGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCGGGAGTAAGCCCATGCTGCTGCGTTTTACCaa                                      <  8:224718/149‑1 (MQ=255)
                                                                                                 cTCAAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTAcc                                        >  4:230439/1‑146 (MQ=255)
                                                                                                    cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGAt                                   <  2:20593/148‑1 (MQ=255)
                                                                                                     aaCTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATg                                  >  3:161075/1‑148 (MQ=255)
                                                                                                        tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtg                              >  8:105147/1‑149 (MQ=255)
                                                                                                          caacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAg                                                                              <  4:278553/99‑1 (MQ=255)
                                                                                                          caacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAg                                                                              >  3:278553/1‑99 (MQ=255)
                                                                                                             cacCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctg                         >  1:204192/1‑149 (MQ=255)
                                                                                                              acCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                        <  2:204192/149‑1 (MQ=255)
                                                                                                               cccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctggg                       <  7:320369/149‑1 (MQ=255)
                                                                                                               cccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctggg                       <  4:285847/149‑1 (MQ=255)
                                                                                                                   cgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGt                                                                                           >  3:250492/1‑77 (MQ=255)
                                                                                                                   cgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGt                                                                                           <  4:250492/77‑1 (MQ=255)
                                                                                                                      ccGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACgcc                >  6:290451/1‑147 (MQ=255)
                                                                                                                      ccGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGAc                <  5:341102/149‑1 (MQ=255)
                                                                                                                        gTTGCGCTTAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCTTTTTACCAAGATGCATGGGCTGGGCAACGACtt              <  3:96828/149‑1 (MQ=255)
                                                                                                                         ttGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGttt                                            >  2:32289/1‑118 (MQ=255)
                                                                                                                         ttGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTCCTGCCGGAGTAAGCCCATGCTGCTGCGttt                                            <  1:32289/118‑1 (MQ=255)
                                                                                                                            cgcTGAAATCCTGGTCGACGGGGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATg          <  5:24575/149‑1 (MQ=255)
                                                                                                                            cgcTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCAt           <  5:290451/148‑1 (MQ=255)
                                                                                                                              cTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGt        >  1:30355/1‑149 (MQ=255)
                                                                                                                               tGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTc       >  6:184313/1‑149 (MQ=255)
                                                                                                                                  aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCt     <  1:198548/148‑1 (MQ=255)
                                                                                                                                  aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCt     <  3:230439/148‑1 (MQ=255)
                                                                                                                                  aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTg    <  7:105147/149‑1 (MQ=255)
                                                                                                                                    tCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGa  >  6:307611/1‑149 (MQ=255)
                                                                                                                                    gccTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGa  >  3:156928/2‑149 (MQ=255)
                                                                                                                                      cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGTCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTgg   >  6:130524/1‑146 (MQ=255)
                                                                                                                                      cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGTCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTgg   <  5:130524/146‑1 (MQ=255)
                                                                                                                                             aCGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctgct                                                  >  7:330185/1‑92 (MQ=255)
                                                                                                                                             aCGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgctgct                                                  <  8:330185/92‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
AGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGA  >  NC_002947/5962132‑5962412

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: