Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG41.8% 99.5 / 30.8 79I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (12/34);  new base G (33/0);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.57e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456965‑4457254
                                                                                                                                             |                                                                                                                                                      
ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGa                                                                                                                                              <  1:977107/150‑1 (MQ=255)
        gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTgg                                                                                                                                     >  2:988104/1‑151 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCATGAACTGCTgg                                                                                                                                     >  2:989160/1‑149 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGAt                                                                                                                                   <  2:405966/151‑1 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGAt                                                                                                                                   <  2:987665/151‑1 (MQ=255)
                ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCCGAAAATTCtgggtggg                                                                                                                                >  1:446543/1‑148 (MQ=255)
                ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTg                                                                                                                             <  2:459455/151‑1 (MQ=255)
                  cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGc                                                                                                                           <  2:925532/151‑1 (MQ=255)
                     aagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACCGCCGGATGGGCTGGCg                                                                                                                          >  2:1080243/1‑149 (MQ=255)
                         cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCgga                                                                                                                        >  2:968851/1‑145 (MQ=255)
                            cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg                                                                                                                 <  1:125671/151‑1 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGAACCATCACGAACTGGTGGGTGGGGTGGGGGAGGAGgggg                                                                                                                >  1:1169792/1‑151 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTCCATCAGAAAACGGTGGGTGGGCGGGCgaagaaagggg                                                                                                                >  1:529087/1‑146 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGGCCCTCAGAAACTTCTGGATGGGGTGGCGAAGAAGGCgg                                                                                                                >  1:263079/1‑151 (MQ=255)
                              ttCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCgg                                                                                                                >  2:857203/1‑150 (MQ=255)
                               tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGa                                                                                                               <  2:917269/150‑1 (MQ=255)
                                 aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGTGGGGGATCATCAGTAACCGCTTGTTGGGGCGGTGGAGAAAGCgg                                                                                                                >  1:90457/1‑147 (MQ=255)
                                   ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTACATCAAAAAATTCTGGGTGGGGTGGCGAAGAAgggg                                                                                                                 >  1:391078/1‑142 (MQ=255)
                                   ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATCAGAAACCGCTGGGTGGGGTGGGGGAGAAGGGGGAAGcc                                                                                                           >  1:16967/1‑150 (MQ=255)
                                      gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                         <  1:96473/149‑1 (MQ=255)
                                          acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGGGAAGAAGGCGGGAGCCgcgc                                                                                                       >  2:918059/1‑147 (MQ=255)
                                             ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCATCAGAAAATGCTGGGTGGGGGGGGGAAGAAGGCGGGAGCCGGGCAGGaaat                                                                                                >  1:771869/1‑149 (MQ=255)
                                               aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGGCCATCAGAAAATGGTGGATGGGGGGGCGGAGAAGGGGGAAGACGCGCAGGaa                                                                                                  >  1:117189/1‑147 (MQ=255)
                                                    cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTGCATCAAAAACTGCTGGGGGGGCGGGGGAAGAAAGGGGGAGGCGCGCAGAAAATCTGcc                                                                                           >  1:1206301/1‑149 (MQ=255)
                                                      aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                       <  2:701603/151‑1 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTTGGCCTTAGcc                                                                                     >  2:441392/1‑151 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGcc                                                                                     <  2:990637/151‑1 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGc                                                                                      <  1:156029/150‑1 (MQ=255)
                                                          ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGcc                                                                                     >  2:561073/1‑149 (MQ=255)
                                                          ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGTGGGGACCATCAGAAAATGCTGGGAGGGCTGGGGAAGAAGGCGGGAGGCGCGCAGGAAATCGGCCTTAGcccg                                                                                   >  1:1246272/1‑149 (MQ=255)
                                                            aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAAAAACTGCTGGGAGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCCTAGcc                                                                                     >  2:504905/1‑147 (MQ=255)
                                                             tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGCCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGACGCGGCGGAATTCCGCCTTAGcc                                                                                     >  2:433952/1‑146 (MQ=255)
                                                               cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGCCGTTAGAAAATTCTGGATAGGCGGGGGGAGAAAGGGGAAGCCCCGCAGTAAGTTGGCCGTAGCCGGgtat                                                                               >  1:27588/1‑150 (MQ=255)
                                                                 cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTTGGCCTTAGCCAGatat                                                                               >  2:100467/1‑148 (MQ=255)
                                                                  tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGGCCAATCAAAACCACGGGAGGGGGTGGCCAAGAACACGGGAGCCGCGCAGGAACTCGGACtcaccgaca                                                                                  >  1:634973/1‑136 (MQ=255)
                                                                   ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCCTCAGAAAATGCTGGATGGGGTGGGGGAGAAGGCGGGAGGCGCGCAGTCATTTGGCCTTTGCCAggaata                                                                              >  1:351586/1‑142 (MQ=255)
                                                                    gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                         <  2:931271/151‑1 (MQ=255)
                                                                    gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                         <  1:927843/151‑1 (MQ=255)
                                                                     cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGc                                                                        <  1:2934/151‑1 (MQ=255)
                                                                      aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTTCATCAGAAAATGCTGGGTGGGCTGGGGGAGAAAGGGGGAGCCCCGCAGGAAATCTGCCCTTGCCCGATATAGTATGCg                                                                       >  1:469616/1‑151 (MQ=255)
                                                                      aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCCTAGCCAGATATAGTATGCg                                                                       >  2:1172351/1‑151 (MQ=255)
                                                                       gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTTGGCCCTAGCCAGATATAGTATGCGt                                                                      >  2:916863/1‑151 (MQ=255)
                                                                        atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGt                                                                      >  2:735103/1‑150 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGGGGAAGCCGCGCAGGAATTTGGCCTTAGCCAGATATAGTATGCGGAgg                                                                   >  1:932229/1‑151 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                                   <  1:270771/151‑1 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAg                                                                    <  2:584401/150‑1 (MQ=255)
                                                                               ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                               <  1:524698/150‑1 (MQ=255)
                                                                                  gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                            <  1:959828/150‑1 (MQ=255)
                                                                                   aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATCCATCAGAAACCGTTGGATGGGCTGGCGCAGAAAGCGGAAGCCGCCCAGAAATTCTGCCCTAGACAGATATAGGATg                                                                         >  2:105214/1‑136 (MQ=255)
                                                                                     aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                        <  2:561720/151‑1 (MQ=255)
                                                                                       ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                      <  2:1035727/151‑1 (MQ=255)
                                                                                        cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGtt                                                     <  2:1074166/151‑1 (MQ=255)
                                                                                         gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                    <  2:77237/151‑1 (MQ=255)
                                                                                          cTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTCCTGGATGGGGTGGCGGAGAAAGCGGGAGCCGCGCCGGAAATCGGCCTTAGCCAGGTATAGTATGCGGAGGCCGGGTAAGGGGtt                                                     >  1:207270/1‑149 (MQ=255)
                                                                                                  ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCCTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgcc                                               >  1:544715/1‑147 (MQ=255)
                                                                                                  ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                           <  1:973893/151‑1 (MQ=255)
                                                                                                    cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATcc                                         >  2:1191303/1‑151 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                         <  2:399629/150‑1 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                         <  2:161735/150‑1 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                        <  1:317867/151‑1 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCCATAAATCGGCCTTAGACCGATATAATATGCGGAGGGCCGGTAAGGCGGTCACGCCCCATcc                                         >  1:750809/1‑150 (MQ=255)
                                                                                                       ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGGAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                       >  2:637169/1‑150 (MQ=255)
                                                                                                        cATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCgg                                       >  1:1107660/1‑149 (MQ=255)
                                                                                                         atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGGAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGGTAAGGCGGTTACGCCGCATCCgg                                       >  1:1189956/1‑148 (MQ=255)
                                                                                                         atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                    <  1:777292/151‑1 (MQ=255)
                                                                                                           atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAttt                                  >  2:494334/1‑151 (MQ=255)
                                                                                                               tGGTGATGAGCAACGGCGGTTTTGGGGGGATCCACTAGAAAATGCTGGGTGGGCTGGCGAAGAAAGCGGGAGCCGCGCCGTAATTCGGCCCTCGCCCGATATAGTATGCGTAGGCCGGAAAAGGCGTTCACCCCGCATCCGGCCttt                                  >  2:1067944/1‑147 (MQ=255)
                                                                                                                ggTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAAGCGGGAGCCGCGCCGTAAATCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACCCCGCCTCCGGCATTTTAg                               >  1:754875/1‑149 (MQ=255)
                                                                                                                  tgatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTTGAGc                              >  2:179227/1‑148 (MQ=255)
                                                                                                                  tgatgaGCAACGGCGGTTTTGGTGGGAGCCATCAGAAAATGCTGGGTGGGCTGGGGGAGAAGGCGGAAGCCGCGCAGTTAATCGGGCTTAGCCCGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCttt                                  >  1:1038492/1‑144 (MQ=38)
                                                                                                                   gatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAg                          <  2:1124568/151‑1 (MQ=255)
                                                                                                                          aaCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGa                   <  2:891115/151‑1 (MQ=255)
                                                                                                                           aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt                  <  1:916324/151‑1 (MQ=255)
                                                                                                                            cggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATg                 <  2:578735/151‑1 (MQ=255)
                                                                                                                             ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGc                <  2:475822/151‑1 (MQ=255)
                                                                                                                                 gTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg            <  1:405891/151‑1 (MQ=255)
                                                                                                                                  ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  2:527480/1‑151 (MQ=255)
                                                                                                                                          gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc    <  2:446549/151‑1 (MQ=255)
                                                                                                                                            gATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGt  >  2:899206/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
GGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456965‑4457254

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: