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

GGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  NC_000913/4456965‑4457254
                                                                                                                                             |                                                                                                                                                      
ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGa                                                                                                                                              <  1:977107/150‑1 (MQ=255)
        gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGGTCC‑ATCAGAAACTGCTgg                                                                                                                                     >  2:988104/1‑151 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GTGGGGGCC‑ATCATGAACTGCTgg                                                                                                                                     >  2:989160/1‑149 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGAt                                                                                                                                   <  2:405966/151‑1 (MQ=255)
          agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGAt                                                                                                                                   <  2:987665/151‑1 (MQ=255)
                ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GTGGGGTCC‑ATCCGAAAATTCtgggtggg                                                                                                                                >  1:446543/1‑148 (MQ=255)
                ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTg                                                                                                                             <  2:459455/151‑1 (MQ=255)
                  cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGc                                                                                                                           <  2:925532/151‑1 (MQ=255)
                     aagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAACCGCCGGATGGGCTGGCg                                                                                                                          >  2:1080243/1‑149 (MQ=255)
                         cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGGTCC‑ATCAGAAACTGCTGGGTGGGCTGGCgga                                                                                                                        >  2:968851/1‑145 (MQ=255)
                            cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg                                                                                                                 <  1:125671/151‑1 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTT‑GTGGGAACC‑ATCACGAACTGGTGGGTGGGGTGGGGGAGGAGgggg                                                                                                                >  1:1169792/1‑151 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTT‑GGGGGGTCC‑ATCAGAAAACGGTGGGTGGGCGGGCgaagaaagggg                                                                                                                >  1:529087/1‑146 (MQ=255)
                             gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTG‑GGGGGGGCC‑CTCAGAAACTTCTGGATGGGGTGGCGAAGAAGGCgg                                                                                                                >  1:263079/1‑151 (MQ=255)
                              ttCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGGTCC‑ATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCgg                                                                                                                >  2:857203/1‑150 (MQ=255)
                               tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGa                                                                                                               <  2:917269/150‑1 (MQ=255)
                                 aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGTGGGGGAT‑C‑ATCAGTAACCGCTTGTTGGGGCGGTGGAGAAAGCgg                                                                                                                >  1:90457/1‑147 (MQ=255)
                                   ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGGTAC‑ATCAAAAAATTCTGGGTGGGGTGGCGAAGAAgggg                                                                                                                 >  1:391078/1‑142 (MQ=255)
                                   ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGGGCC‑ATCAGAAACCGCTGGGTGGGGTGGGGGAGAAGGGGGAAGcc                                                                                                           >  1:16967/1‑150 (MQ=255)
                                      gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                         <  1:96473/149‑1 (MQ=255)
                                          acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAAATGCTGGGTGGGGTGGGGAAGAAGGCGGGAGCCgcgc                                                                                                       >  2:918059/1‑147 (MQ=255)
                                             ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGGGCC‑ATCAGAAAATGCTGGGTGGGGGGGGGAAGAAGGCGGGAGCCGGGCAGGaaat                                                                                                >  1:771869/1‑149 (MQ=255)
                                               aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTT‑GGGGGGGCC‑ATCAGAAAATGGTGGATGGGGGGGCGGAGAAGGGGGAAGACGCGCAGGaa                                                                                                  >  1:117189/1‑147 (MQ=255)
                                                    cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGGTGC‑ATCAAAAACTGCTGGGGGGGCGGGGGAAGAAAGGGGGAGGCGCGCAGAAAATCTGcc                                                                                           >  1:1206301/1‑149 (MQ=255)
                                                      aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                       <  2:701603/151‑1 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTTGGCCTTAGcc                                                                                     >  2:441392/1‑151 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGcc                                                                                     <  2:990637/151‑1 (MQ=255)
                                                        gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGc                                                                                      <  1:156029/150‑1 (MQ=255)
                                                          ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGcc                                                                                     >  2:561073/1‑149 (MQ=255)
                                                          ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTT‑GTGGGGACC‑ATCAGAAAATGCTGGGAGGGCTGGGGAAGAAGGCGGGAGGCGCGCAGGAAATCGGCCTTAGcccg                                                                                   >  1:1246272/1‑149 (MQ=255)
                                                            aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAAAAACTGCTGGGAGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCCTAGcc                                                                                     >  2:504905/1‑147 (MQ=255)
                                                             tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGGCCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGACGCGGCGGAATTCCGCCTTAGcc                                                                                     >  2:433952/1‑146 (MQ=255)
                                                               cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTT‑GTGGGGGCC‑GTTAGAAAATTCTGGATAGGCGGGGGGAGAAAGGGGAAGCCCCGCAGTAAGTTGGCCGTAGCCGGgtat                                                                               >  1:27588/1‑150 (MQ=255)
                                                                 cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTTGGCCTTAGCCAGatat                                                                               >  2:100467/1‑148 (MQ=255)
                                                                  tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTG‑GTGGGGGCCAATCA‑AAACCACGGGAGGGGGTGGCCAAGAACACGGGAGCCGCGCAGGAACTCGGACtcaccgaca                                                                                  >  1:634973/1‑136 (MQ=255)
                                                                   ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGGGCC‑CTCAGAAAATGCTGGATGGGGTGGGGGAGAAGGCGGGAGGCGCGCAGTCATTTGGCCTTTGCCAggaata                                                                              >  1:351586/1‑142 (MQ=255)
                                                                    gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                         <  2:931271/151‑1 (MQ=255)
                                                                    gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                         <  1:927843/151‑1 (MQ=255)
                                                                     cAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGc                                                                        <  1:2934/151‑1 (MQ=255)
                                                                      aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGGTTC‑ATCAGAAAATGCTGGGTGGGCTGGGGGAGAAAGGGGGAGCCCCGCAGGAAATCTGCCCTTGCCCGATATAGTATGCg                                                                       >  1:469616/1‑151 (MQ=255)
                                                                      aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGGTCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCCTAGCCAGATATAGTATGCg                                                                       >  2:1172351/1‑151 (MQ=255)
                                                                       gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTTGGCCCTAGCCAGATATAGTATGCGt                                                                      >  2:916863/1‑151 (MQ=255)
                                                                        atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGt                                                                      >  2:735103/1‑150 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGGGGAAGCCGCGCAGGAATTTGGCCTTAGCCAGATATAGTATGCGGAgg                                                                   >  1:932229/1‑151 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAgg                                                                   <  1:270771/151‑1 (MQ=255)
                                                                          atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAg                                                                    <  2:584401/150‑1 (MQ=255)
                                                                               ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                               <  1:524698/150‑1 (MQ=255)
                                                                                  gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                            <  1:959828/150‑1 (MQ=255)
                                                                                   aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GGGGGATCC‑ATCAGAAACCGTTGGATGGGCTGGCGCAGAAAGCGGAAGCCGCCCAGAAATTCTGCCCTAGACAGATATAGGATg                                                                         >  2:105214/1‑136 (MQ=255)
                                                                                     aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                        <  2:561720/151‑1 (MQ=255)
                                                                                       ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                      <  2:1035727/151‑1 (MQ=255)
                                                                                        cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGtt                                                     <  2:1074166/151‑1 (MQ=255)
                                                                                         gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                    <  2:77237/151‑1 (MQ=255)
                                                                                          cTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GTGGGGGCC‑ATCAGAAACTCCTGGATGGGGTGGCGGAGAAAGCGGGAGCCGCGCCGGAAATCGGCCTTAGCCAGGTATAGTATGCGGAGGCCGGGTAAGGGGtt                                                     >  1:207270/1‑149 (MQ=255)
                                                                                                  ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTT‑GTGGGGTCC‑CTCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCCTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgcc                                               >  1:544715/1‑147 (MQ=255)
                                                                                                  ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                           <  1:973893/151‑1 (MQ=255)
                                                                                                    cGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATcc                                         >  2:1191303/1‑151 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                         <  2:399629/150‑1 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                         <  2:161735/150‑1 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                        <  1:317867/151‑1 (MQ=255)
                                                                                                     gACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GGGGGGTCC‑ATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCCATAAATCGGCCTTAGACCGATATAATATGCGGAGGGCCGGTAAGGCGGTCACGCCCCATcc                                         >  1:750809/1‑150 (MQ=255)
                                                                                                       ccATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGGAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                       >  2:637169/1‑150 (MQ=255)
                                                                                                        cATATTCTGGTGATGAGCAACGGCGGTTTTT‑GTGGGGTCC‑ATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCgg                                       >  1:1107660/1‑149 (MQ=255)
                                                                                                         atatTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGGTCC‑ATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGGAGGCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGGTAAGGCGGTTACGCCGCATCCgg                                       >  1:1189956/1‑148 (MQ=255)
                                                                                                         atatTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                    <  1:777292/151‑1 (MQ=255)
                                                                                                           atTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAttt                                  >  2:494334/1‑151 (MQ=255)
                                                                                                               tGGTGATGAGCAACGGCGGTTTTG‑GGGGGATCC‑ACTAGAAAATGCTGGGTGGGCTGGCGAAGAAAGCGGGAGCCGCGCCGTAATTCGGCCCTCGCCCGATATAGTATGCGTAGGCCGGAAAAGGCGTTCACCCCGCATCCGGCCttt                                  >  2:1067944/1‑147 (MQ=255)
                                                                                                                ggTGATGAGCAACGGCGGTTTTT‑GTGGGGTCC‑ATCAGAAACTGCTGGATGGGCTGGCGGAGAAAGCGGGAGCCGCGCCGTAAATCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACCCCGCCTCCGGCATTTTAg                               >  1:754875/1‑149 (MQ=255)
                                                                                                                  tgatgaGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTTGAGc                              >  2:179227/1‑148 (MQ=255)
                                                                                                                  tgatgaGCAACGGCGGTTTTG‑GTGGGAGCC‑ATCAGAAAATGCTGGGTGGGCTGGGGGAGAAGGCGGAAGCCGCGCAGTTAATCGGGCTTAGCCCGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCttt                                  >  1:1038492/1‑144 (MQ=38)
                                                                                                                   gatgaGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAg                          <  2:1124568/151‑1 (MQ=255)
                                                                                                                          aaCGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGa                   <  2:891115/151‑1 (MQ=255)
                                                                                                                           aCGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt                  <  1:916324/151‑1 (MQ=255)
                                                                                                                            cggcggTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATg                 <  2:578735/151‑1 (MQ=255)
                                                                                                                             ggcggTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGc                <  2:475822/151‑1 (MQ=255)
                                                                                                                                 gTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg            <  1:405891/151‑1 (MQ=255)
                                                                                                                                  ttttG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc           >  2:527480/1‑151 (MQ=255)
                                                                                                                                          gggATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc    <  2:446549/151‑1 (MQ=255)
                                                                                                                                            gATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGt  >  2:899206/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
GGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTG‑GTGGGATCC‑ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGT  >  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: ‑