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

CATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAT  >  NC_000913/4456955‑4457250
                                                                                                                                                      |                                                                                                                                                 
cATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGa                                                                                                                                                   <  2:421396/151‑1 (MQ=255)
      ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATc                                                                                                                                             <  2:592497/151‑1 (MQ=255)
          ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGGGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATc                                                                                                                                             >  1:943939/1‑147 (MQ=255)
           gcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTCCATCAGaa                                                                                                                                         >  2:531028/1‑150 (MQ=255)
             aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTctggaaa                                                                                                                                       >  1:242049/1‑144 (MQ=255)
             aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct                                                                                                                                      <  2:128450/151‑1 (MQ=255)
             aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct                                                                                                                                      <  2:731412/151‑1 (MQ=255)
                tggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCCGAAActg                                                                                                                                     >  2:194898/1‑149 (MQ=255)
                    agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCCTCAGAAAATGCTGGa                                                                                                                                >  2:799601/1‑150 (MQ=255)
                    agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTCCATCAGAAAATGCTgg                                                                                                                                 >  1:949487/1‑149 (MQ=255)
                        gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGACCATCAGGAACTTCCGGATggg                                                                                                                            >  1:1012539/1‑150 (MQ=255)
                          ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCt                                                                                                                          <  1:821850/150‑1 (MQ=255)
                            cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACCTCTGggtggggtgg                                                                                                                        >  2:537483/1‑150 (MQ=255)
                            cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAAATGCTGGATGGGGTgg                                                                                                                        >  2:254631/1‑150 (MQ=255)
                            cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCCTCAGAAAATGCtgggtgggggggg                                                                                                                       >  1:417475/1‑146 (MQ=255)
                                 gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaaga                                                                                                                  <  1:77767/151‑1 (MQ=255)
                                 gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCTCAGAAAATGCTGGGTGGGGTGGCg                                                                                                                      >  1:379162/1‑147 (MQ=255)
                                  ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa                                                                                                                 <  2:1026158/151‑1 (MQ=255)
                                   cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCgaaga                                                                                                                  >  2:238389/1‑149 (MQ=255)
                                      cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGCTATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGTGGGGACCAACACAACATCCTCGAAGGGCTGGGgaaaaaaagg                                                                                                             >  1:1038914/1‑147 (MQ=255)
                                          cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAg                                                                                                         <  1:887354/151‑1 (MQ=255)
                                              cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTTCATCCGAAAATGCTGGGTGGGGTGGCGGAGAAGGCGGGAGCcgc                                                                                                     >  2:919277/1‑151 (MQ=255)
                                                gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg                                                                                                    <  1:958699/150‑1 (MQ=255)
                                                   cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGt                                                                                                <  2:666532/151‑1 (MQ=255)
                                                    acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAACGCTGGGTGGGGGGGCGGAGAAAGAGGGAGCCgcgcgg                                                                                                 >  1:572403/1‑147 (MQ=255)
                                                     cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCCTCAGAAACTGCTGGGTGGGGGGGGGAAGAAGGGGGGAGGCGCGCAgga                                                                                               >  1:1093160/1‑148 (MQ=255)
                                                     cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGTGGGGTCCATTAGAAAACGCTGGATGGGGTGGGGGAGAAAGGGGGAGGCGCGCAg                                                                                                 >  2:1027041/1‑148 (MQ=255)
                                                       ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTaa                                                                                              <  1:736258/149‑1 (MQ=255)
                                                        gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                           <  1:860037/151‑1 (MQ=255)
                                                        gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                           <  1:732039/151‑1 (MQ=255)
                                                              cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGATCCATCAGAAAATGCTGGATGGGGGGGCGGAGAAGGCGGGAGGCGCGCCGTAATTTGGGCtt                                                                                     >  2:245573/1‑151 (MQ=255)
                                                              cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCtt                                                                                     <  1:102588/151‑1 (MQ=255)
                                                               gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGccct                                                                                     >  2:594673/1‑148 (MQ=255)
                                                                                aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAAAAAATGCTGGGTGGGCTGGCGAAGAAAGGGGAAGACGCGCAGGAATTTGGCCTTAGCCAGATATAGTATGCg                                                                   >  1:511028/1‑151 (MQ=255)
                                                                                  atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCg                                                                   >  1:927393/1‑149 (MQ=255)
                                                                                     tggtggTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGCAGCCGCGCAGTAAATTGGGCTTTGCCAGACATAGAATGCGTAg                                                                >  1:305176/1‑149 (MQ=255)
                                                                                        tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAAGGGGAAACCGCGCCGTAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCgg                                                           >  2:190270/1‑151 (MQ=255)
                                                                                         ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAAGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                           >  2:876616/1‑150 (MQ=255)
                                                                                          gtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGAt                                                         <  1:942576/151‑1 (MQ=255)
                                                                                            gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCCTCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGCGGGAGGCGCGCAGTAATTCGGCCCTAGCCCGATATAGTATGCGGAGGCCgggca                                                        >  2:617277/1‑147 (MQ=255)
                                                                                             aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                      <  2:1090089/151‑1 (MQ=255)
                                                                                               aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                    <  2:1039970/151‑1 (MQ=255)
                                                                                               aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                    <  1:737106/151‑1 (MQ=255)
                                                                                               aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                    <  1:322373/151‑1 (MQ=255)
                                                                                                 ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                  <  2:185905/151‑1 (MQ=255)
                                                                                                     tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGGGAAGAAGGCGGAAGCCGCGCAGTAATTTGGCCCTAGCCAGATATAGTATGGGTAGGCCGGGTAAgggggtccc                                              >  1:1199402/1‑144 (MQ=255)
                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCa                                               <  1:912972/149‑1 (MQ=255)
                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcg                                             <  1:1067951/151‑1 (MQ=255)
                                                                                                      cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                              <  2:601460/150‑1 (MQ=255)
                                                                                                        gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgcc                                           <  1:637304/151‑1 (MQ=255)
                                                                                                        gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAAGCCGGAGCCCCGCAGGAATTCGGCCTTAGCCCAATATAGTATGCGTAGGCCGGGTAAGGGGGTTAcgc                                            >  2:289337/1‑150 (MQ=255)
                                                                                                         ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccg                                          <  1:846771/151‑1 (MQ=255)
                                                                                                         ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGAGGGATACATCAAAAAATGGTGGAGGGGATGGCGAAGGAGGCGGAAGCCGCGCAGCAATTCGGCCCTAGGCAGCTATTGTATGCGGAGGGCGGGGAAGGCGGTCAagcc                                           >  2:215169/1‑150 (MQ=255)
                                                                                                          cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGGAATTTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCAcgccgc                                         >  2:716845/1‑151 (MQ=255)
                                                                                                           tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGATGGGGTGGCGAAGAAAGGGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAATATGCGTAGGCCGGGTAAGGCGTTCAcgccgc                                         >  1:993817/1‑150 (MQ=255)
                                                                                                                aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCCTAGCCCGGTATAGTATGCGTAGGCCGGGTAAAGCGGTCACGCCCGATCCgg                                   >  2:379434/1‑151 (MQ=255)
                                                                                                                 ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                   <  1:858476/150‑1 (MQ=255)
                                                                                                                   atatTCTGGTGATGAGCAACGGCGGGTTTGGTGGGAACCCTCAGAAACTGAGGGATGGGGTGGGGAAGAAGGCGGAAGGCGCCCAGGAATTAGGCCTTAGCCAGATATAGAATGAGTCGGCCCGATAAGGGGTTCACCCCCCAccc                                     >  2:1194848/1‑146 (MQ=255)
                                                                                                                    tatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGGAGCCGCGCCGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTTCCGCCGCATCCGGCAtt                               >  1:598315/1‑151 (MQ=255)
                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                              <  2:885731/151‑1 (MQ=255)
                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                               <  2:938325/150‑1 (MQ=255)
                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                               <  1:263489/150‑1 (MQ=255)
                                                                                                                     atTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                              >  2:13362/1‑151 (MQ=255)
                                                                                                                       tCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGTGGGGTGGCGGAGAAGGGGGGAGCCGCGCAGTAATTCTGCCTTTGCCCGATATAGTATGCGGAGGCCGGGTAAAGGGGTCACGCCGCATCCCGCAttt                              >  1:567341/1‑149 (MQ=255)
                                                                                                                       tCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTg                             <  2:1006562/150‑1 (MQ=255)
                                                                                                                              atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGt                     <  1:1040766/151‑1 (MQ=255)
                                                                                                                              atgaGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAAGGGGAAGCCGCGCAGGAAATTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTTAGCATAg                      >  1:1109472/1‑150 (MQ=37)
                                                                                                                                  gCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                 <  1:55543/151‑1 (MQ=255)
                                                                                                                                      cggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATg             <  1:820792/151‑1 (MQ=255)
                                                                                                                                           gTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg        <  2:227281/151‑1 (MQ=255)
                                                                                                                                                 gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt  <  2:897017/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAT  >  NC_000913/4456955‑4457250

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