Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,827,659 T→C 9.4% V191A (GTG→GCG)  mtgA → monofunctional biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,827,6590TC 9.4% 141.8 / 5.7 53V191A (GTG→GCG) mtgAmonofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base C (3/2);  total (30/23)
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.06e-01

CTGGTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGC  >  NC_002947/5827513‑5827806
                                                                                                                                                  |                                                                                                                                                   
cTGGTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGtg                                                                                                                                                    <  6:96667/148‑1 (MQ=255)
   gTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGtgctgc                                                                                                                                                <  3:12606/149‑1 (MQ=255)
   gTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGtgctg                                                                                                                                                 <  8:137926/148‑1 (MQ=255)
   gTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCACCGCTGCCAgctgcgctgc                                                                                                                                                <  8:234590/149‑1 (MQ=255)
   gTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCACCGCTGCCAgctgcgctgc                                                                                                                                                <  2:41942/149‑1 (MQ=255)
       aaGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAg                                                                                                                                            >  2:148263/1‑149 (MQ=255)
       aaGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAg                                                                                                                                            >  3:210365/1‑149 (MQ=255)
        aGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGc                                                                                                                                           <  6:226359/149‑1 (MQ=255)
        aGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGc                                                                                                                                           >  6:84861/1‑149 (MQ=255)
        aGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGc                                                                                                                                           >  1:255622/1‑149 (MQ=255)
         ggAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGcc                                                                                                                                          <  3:262648/149‑1 (MQ=255)
           aaCGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGccc                                                                                                                                         >  2:105693/1‑148 (MQ=255)
            aCGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCaa                                                                                                                                       >  4:113489/1‑149 (MQ=255)
                 ttCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAg                                                                                                                                  <  5:84861/149‑1 (MQ=255)
                    tcGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGTCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTg                                                                                                                                >  7:178693/1‑148 (MQ=255)
                         gTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGc                                                                                                                          <  3:113489/149‑1 (MQ=255)
                          tCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGc                                                                                                                          >  3:289864/1‑148 (MQ=255)
                            tATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCACCGCTGCCag                                                                                                                                             <  2:51944/127‑2 (MQ=255)
                            tATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCACCGCTGCCag                                                                                                                                             >  1:51944/1‑126 (MQ=255)
                                 gAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTg                                                                                                                                <  2:153888/135‑1 (MQ=255)
                                 gAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTg                                                                                                                                >  1:153888/1‑135 (MQ=255)
                                 gAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGccgac                                                                                                                  >  7:200191/1‑149 (MQ=255)
                                  aaCAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGccgacc                                                                                                                 >  6:109396/1‑149 (MQ=255)
                                        gcCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGc                                                                                                           <  3:274972/149‑1 (MQ=255)
                                          cGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGTCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCt                                                                                                         <  8:178693/149‑1 (MQ=255)
                                                  gCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGc                                                                                                  >  8:224935/1‑148 (MQ=255)
                                                     aGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgc                                                                                              >  7:54258/1‑149 (MQ=255)
                                                        gCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgccgg                                                                                           <  4:210365/149‑1 (MQ=255)
                                                                ggTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTg                                                                                   >  8:146141/1‑149 (MQ=255)
                                                                     ccAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTc                                                                              >  1:250746/1‑149 (MQ=255)
                                                                       aGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGc                                                                            <  8:54258/149‑1 (MQ=255)
                                                                       aGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGCGCCGGCTGGATTCGc                                                                            >  4:276186/1‑149 (MQ=255)
                                                                            gcggcgCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCa                                                                       <  7:224935/149‑1 (MQ=255)
                                                                               gcgcGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGccgaccga                                                                                                               <  4:29296/106‑1 (MQ=255)
                                                                               gcgcGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGccgaccga                                                                                                               >  3:29296/1‑106 (MQ=255)
                                                                                   gTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGc                                                                >  6:75855/1‑149 (MQ=255)
                                                                                         aCTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCt                                                           <  7:146141/148‑1 (MQ=255)
                                                                                         aCTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTg                                                          >  2:303222/1‑149 (MQ=255)
                                                                                          ggctGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgg                                                         <  5:273085/146‑1 (MQ=255)
                                                                                               gCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcgg                                                     <  2:250746/148‑1 (MQ=255)
                                                                                                     aCGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGc                                              >  7:264220/1‑149 (MQ=255)
                                                                                                     aCGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGc                                              <  3:140364/149‑1 (MQ=255)
                                                                                                       gCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTa                                            >  8:140215/1‑149 (MQ=255)
                                                                                                        ccAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGGGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCCCGCCCAGCtgg                                           <  1:148263/149‑3 (MQ=255)
                                                                                                            ccGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGa                                       <  4:289864/149‑1 (MQ=255)
                                                                                                                 tCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAg                                  >  8:218636/1‑149 (MQ=255)
                                                                                                                     ccgccAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTg                               >  3:148468/1‑148 (MQ=255)
                                                                                                                          gcttggcagCACCGCTGCCAGCTGCGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActcc                         >  5:27098/5‑149 (MQ=255)
                                                                                                                             agggagCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                      <  2:255622/149‑1 (MQ=255)
                                                                                                                              ggcagCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActcctc                       >  3:3923/1‑147 (MQ=255)
                                                                                                                              ggcagCACCGCTGCCAGCTGCGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                      <  6:27098/148‑1 (MQ=255)
                                                                                                                                   cgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCa                                                >  2:159908/1‑117 (MQ=255)
                                                                                                                                       cttGGCAGCACCGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt            >  5:248961/3‑149 (MQ=255)
                                                                                                                                             aGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGc      <  5:109396/149‑1 (MQ=255)
                                                                                                                                               cGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAg    <  5:75855/149‑1 (MQ=255)
                                                                                                                                                ggTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAgg   >  3:220071/1‑149 (MQ=255)
                                                                                                                                                 gTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggc  >  7:56704/1‑149 (MQ=255)
                                                                                                                                                 gTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggc  <  8:200191/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGGTCGAAGGAACGGATTCTCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGC  >  NC_002947/5827513‑5827806

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