Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,925 G→A 6.6% A140V (GCC→GTC)  dapC ← N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,780,9250GA 6.6% 178.3 / 2.5 60A140V (GCC→GTC) dapCN‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTGGCGGTCAGC  >  NC_002947/1780779‑1781073
                                                                                                                                                  |                                                                                                                                                    
gggTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCa                                                                                                                                                    >  7:267095/1‑149 (MQ=255)
    gtTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCgg                                                                                                                                                >  1:59384/1‑149 (MQ=255)
           ggCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCcagca                                                                                                                                          <  6:244508/148‑1 (MQ=255)
            gCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGACACCGGCCagcagc                                                                                                                                        >  5:297332/1‑149 (MQ=255)
                cAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTATGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTg                                                                                                                                    >  6:112218/1‑149 (MQ=255)
                  gcgcACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGc                                                                                                                                               <  7:304191/136‑1 (MQ=255)
                  gcgcACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGc                                                                                                                                               >  8:304191/1‑136 (MQ=255)
                       ccGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTc                                                                                                                             <  3:233303/149‑1 (MQ=255)
                        cGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCg                                                                                                                            >  6:71260/1‑149 (MQ=255)
                        cGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTGGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCg                                                                                                                            >  7:297420/1‑149 (MQ=255)
                         ggTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGt                                                                                                                           <  7:4472/149‑1 (MQ=255)
                              ggTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATc                                                                                                                      >  5:169898/1‑149 (MQ=255)
                               gTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCt                                                                                                                     <  4:85981/149‑1 (MQ=255)
                                  gCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTgg                                                                                                                   <  7:106978/148‑1 (MQ=255)
                                  gCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTgg                                                                                                                   <  6:169898/148‑1 (MQ=255)
                                    cTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTgg                                                                                                                   >  7:90618/1‑146 (MQ=255)
                                      ggCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGaa                                                                                                              >  7:14090/1‑149 (MQ=255)
                                       gCGAGCACAGGAACACTATCTGGCAACGCATCCACAACTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGCGGTGGCACCGGCAAGCAGCGCTTCGCCTTCGTAGATCTGGTAGa                                                                                                               >  5:222752/1‑147 (MQ=255)
                                           gCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGTCTTCGTAGATCTGGTAGAACGGGt                                                                                                         <  7:44168/149‑1 (MQ=255)
                                              caggaacaggaTCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTgg                                                                                                      >  3:84874/1‑149 (MQ=255)
                                                   acaggaTCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGACACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTg                                                                                                 <  6:297332/149‑1 (MQ=255)
                                                    caggaTCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGa                                                                                                >  5:313774/1‑149 (MQ=255)
                                                            ggCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTgggg                                                                                                    >  5:194199/1‑137 (MQ=255)
                                                            ggCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTgggg                                                                                                    <  6:194199/137‑1 (MQ=255)
                                                            ggCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAgg                                                                                        >  2:124009/1‑149 (MQ=255)
                                                             gCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTCGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGc                                                                                       >  8:168925/1‑149 (MQ=255)
                                                              gaaggCTTCCACACCTCGGCCGGTTCGGCGTCGAAGTCCGGGTTGTAACCGTTGTTTTGCAGGCAGGGCATGTAGTGCGGGGTGGCAGCGGCCAGGAGCGCTGCGCCTTCGTAGATGTGGTAGAACGGGTTGGGGCTGATCACCAGGcc                                                                                      <  7:188577/145‑1 (MQ=255)
                                                                      ccACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCgg                                                                              >  3:21247/1‑149 (MQ=255)
                                                                              cGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGtt                                                                       >  3:235823/1‑148 (MQ=255)
                                                                                gCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGAAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGa                                                                     >  3:92539/1‑148 (MQ=255)
                                                                                 ccGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTc                                                                                                                             <  3:271815/91‑1 (MQ=255)
                                                                                 ccGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTc                                                                                                                             >  4:271815/1‑91 (MQ=255)
                                                                                  cGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCa                                                                  >  2:42425/1‑149 (MQ=255)
                                                                                       cggcggCGAAGTGCGGGTTGAATCCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGTGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTTGTAGAACGGGTTGGGGCTGATCTCCAGGCCGTCATCGGCGCGGTTGACCACGGcc                                                             <  6:222752/149‑1 (MQ=255)
                                                                                         gCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTg                                                           <  6:313774/149‑1 (MQ=255)
                                                                                           gTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTggg                                                         >  4:167856/1‑149 (MQ=255)
                                                                                           gTCGAAGTCCGGGTCGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTggg                                                         <  7:168925/149‑1 (MQ=255)
                                                                                            tCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGt                                                        >  8:88100/1‑149 (MQ=255)
                                                                                              gAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGa                                                      <  1:124009/149‑1 (MQ=255)
                                                                                               aaGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGaa                                                     <  8:90618/149‑1 (MQ=255)
                                                                                                 gTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAgg                                                   <  1:42425/149‑1 (MQ=255)
                                                                                                 gTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAgg                                                   <  4:21247/149‑1 (MQ=255)
                                                                                                            aaCCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCt                                        >  6:795/1‑149 (MQ=255)
                                                                                                            aaCCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCt                                        >  1:50804/1‑149 (MQ=255)
                                                                                                              ccGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTc                                      <  7:88100/149‑1 (MQ=255)
                                                                                                              ccGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGACACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCtt                                       >  7:81469/1‑148 (MQ=255)
                                                                                                                      ttCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGc                              >  1:203786/1‑149 (MQ=255)
                                                                                                                      ttCCAGGCAGGGCAGGTAGTGCGGGGTGACACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTg                               <  8:81469/148‑1 (MQ=255)
                                                                                                                        ccAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCg                            >  8:174338/1‑149 (MQ=255)
                                                                                                                         caggcaggGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGt                           <  4:84874/149‑1 (MQ=255)
                                                                                                                           ggcagggcagGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGtt                          >  7:22856/1‑148 (MQ=255)
                                                                                                                            gcagggcaggTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGa                        >  1:128593/1‑149 (MQ=255)
                                                                                                                                ggcaggTAGTGCGGGGGGGCACCGGCCAGCAGCGCTGCGCCTTCGTCGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCAGCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCgg                    <  8:297420/149‑1 (MQ=255)
                                                                                                                                 gcaggTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCggg                   <  8:14090/149‑1 (MQ=255)
                                                                                                                                  caggTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGgga                  <  2:50804/149‑1 (MQ=255)
                                                                                                                                          gCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTgg          >  2:60753/1‑149 (MQ=255)
                                                                                                                                          gCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTgg          >  6:160029/1‑149 (MQ=255)
                                                                                                                                           cGGGGTGGCACCGTCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAAGGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTgg          <  7:174338/148‑1 (MQ=255)
                                                                                                                                               gTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTGGCGGt      >  4:21538/1‑148 (MQ=255)
                                                                                                                                                tGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACTGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTGGCGGTCa    <  8:22856/149‑1 (MQ=255)
                                                                                                                                                  gCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTGGCGGTCAGc  >  5:54221/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTTGTTTTCCAGGCAGGGCAGGTAGTGCGGGGTGGCACCGGCCAGCAGCGCTGCGCCTTCGTAGATCTGGTAGAACGGGTTGGGGCTGATCACCAGGCCGTCATCGGCGCGGTTGACCACGGCCTGGGTGAAGGCGAACAGCGCTTCACGGGTGCCGTTGACCGGGAGGATGTGGCGGTCAGC  >  NC_002947/1780779‑1781073

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