Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG46.0% 148.7 / 44.5 113I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (10/51);  new base G (52/0);  total (62/51)
Fisher's exact test for biased strand distribution p-value = 2.19e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456956‑4457251
                                                                                                                                                       |                                                                                                                                                  
    tCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCa                                                                                                                                                <  1:1600771/150‑1 (MQ=255)
         ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGa                                                                                                                                           <  1:1302017/150‑1 (MQ=255)
          gcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGaaa                                                                                                                                         <  2:617238/151‑1 (MQ=255)
             ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActg                                                                                                                                      <  1:884766/151‑1 (MQ=255)
                  cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAaaaa                                                                                                                                         >  1:1690888/1‑143 (MQ=255)
                    gaagtggcagaagCCTGCGTTCAGCTTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                               <  2:364256/151‑1 (MQ=255)
                    gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAATGCTGGGTg                                                                                                                               >  2:565111/1‑151 (MQ=255)
                    gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGGACATCAGGAAATGCTgg                                                                                                                                  >  1:1248641/1‑148 (MQ=255)
                    gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGTCCATCAGAAAATGCTGGATg                                                                                                                               >  2:635459/1‑151 (MQ=255)
                     aagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAACTGCTGGATgg                                                                                                                              >  1:1576439/1‑151 (MQ=255)
                      agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCtgggtggg                                                                                                                             >  2:355341/1‑151 (MQ=255)
                         ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTg                                                                                                                          <  2:990172/151‑1 (MQ=255)
                            agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTGCTGGATGGGCTgggg                                                                                                                       >  1:690150/1‑149 (MQ=255)
                                gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAAATGCTGGGTGGGCTGGCg                                                                                                                       >  2:1129631/1‑147 (MQ=255)
                                 ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCgaaga                                                                                                                   >  2:1494350/1‑150 (MQ=255)
                                 ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCCGAAACTGCTGGAGGGGCTgg                                                                                                                         >  1:1343298/1‑144 (MQ=255)
                                  cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATCAAAAAATGCTGGATggg                                                                                                                             >  1:1032800/1‑139 (MQ=255)
                                  cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGATCCAACAGAAAATGCTGGATGGGGTGGCggagaa                                                                                                                  >  2:1062516/1‑150 (MQ=255)
                                     cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg                                                                                                              <  1:1331734/151‑1 (MQ=255)
                                      gTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCgg                                                                                                             >  2:932971/1‑151 (MQ=255)
                                       ttCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCCGGGTGGGCGGGGGGAGAAGgggg                                                                                                             >  1:1710750/1‑150 (MQ=255)
                                        tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGa                                                                                                            <  2:692379/150‑1 (MQ=255)
                                          aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTTCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAg                                                                                                          >  2:343695/1‑150 (MQ=255)
                                          aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGTCATCAGAAACTTCGGGGTGGGGCGGGGGAGAAGGCgg                                                                                                             >  1:642346/1‑147 (MQ=255)
                                          aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGCCCTCAGAAAATGCTGGATGGGGTGGGgaagaaagggggag                                                                                                          >  1:1016802/1‑142 (MQ=255)
                                               gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg                                                                                                     <  2:702729/150‑1 (MQ=255)
                                                  cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTTCATCAGAAACTGCTTGGTGGGCTGGCGGAAAAAGCGGAAGGcgcg                                                                                                     >  2:40842/1‑147 (MQ=255)
                                                  cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGCTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAGGCGGAAggc                                                                                                        >  2:923611/1‑142 (MQ=255)
                                                   acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGATCCTCAAAAAATGCCGGGAGGGGCGGCGGAGAAGGGGGGAGACgcgccg                                                                                                  >  1:467227/1‑147 (MQ=255)
                                                      ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAt                                                                                              <  2:925776/150‑1 (MQ=255)
                                                           ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCCGGAATTTGGcc                                                                                        >  2:1699862/1‑151 (MQ=255)
                                                             cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTGCATCAGAAAATGCTGGGTGGGGTGGGGGAGAAGGGGGGAGCCGCGCAGGAATTAGGGCtt                                                                                      >  1:696260/1‑151 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                     <  1:1059764/151‑1 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                     <  2:1115728/151‑1 (MQ=255)
                                                              gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa                                                                                     <  2:196878/151‑1 (MQ=255)
                                                               aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg                                                                                    <  2:640690/151‑1 (MQ=255)
                                                                  tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGGGGAGAAGGCGGGAGCCGCGCAGTAAATCGGCCTTAGcc                                                                                  >  2:646754/1‑150 (MQ=255)
                                                                   ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCCATCGAAAATGGTGGATGGGGGGGCGGAGAAAGGGGAAGCGGGGCCGTAAttgcggcttt                                                                                     >  1:159533/1‑138 (MQ=255)
                                                                    gATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGa                                                                               <  2:1679113/151‑1 (MQ=255)
                                                                    gATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGa                                                                               <  2:1316589/151‑1 (MQ=255)
                                                                       aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGatat                                                                            <  1:1406203/151‑1 (MQ=255)
                                                                       aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata                                                                             <  1:591228/150‑1 (MQ=255)
                                                                       aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata                                                                             <  2:921075/150‑1 (MQ=255)
                                                                          cTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTTCATCAGAAAATGCTGGGTGGGGTGGCGAAGAAGGCGGAAGCCGCGGCGTAATTTGGCCTTAGCCAGa                                                                               >  2:805943/1‑145 (MQ=255)
                                                                            ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGt                                                                         <  2:967761/149‑1 (MQ=255)
                                                                            ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTa                                                                        <  1:872083/150‑1 (MQ=255)
                                                                            ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                       <  2:1353350/151‑1 (MQ=255)
                                                                               aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGGGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCg                                                                    >  2:386657/1‑151 (MQ=255)
                                                                                gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCg                                                                    <  2:1261707/150‑1 (MQ=255)
                                                                                 atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCAACAGAAAATGCTGGGTGGGGTGGCGAAGAAAGCGGGAGGCGCGCAGGAATTTGGCCTTAGCCAGATATAATATGCg                                                                    >  1:1595737/1‑149 (MQ=255)
                                                                                  tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCCTAGCCCGATATAGTATGCGt                                                                   >  2:479255/1‑149 (MQ=255)
                                                                                  tatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGt                                                                   >  2:575689/1‑149 (MQ=255)
                                                                                   atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAg                                                                 <  2:1357004/150‑1 (MQ=255)
                                                                                     ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGc                                                               <  2:353878/150‑1 (MQ=255)
                                                                                       tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAAAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                            >  1:265347/1‑151 (MQ=255)
                                                                                        ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCTGGATGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCTGCCTTAGCCAGATATAGTATGCGGAGGCCgg                                                            >  1:1339701/1‑150 (MQ=255)
                                                                                        ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                            <  1:1336153/150‑1 (MQ=255)
                                                                                        ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGa                                                           <  1:598505/151‑1 (MQ=255)
                                                                                           gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                         >  2:1626046/1‑150 (MQ=255)
                                                                                            aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATTCTGGATGGGGTGGGGGAGGAGGCGGGAGGCGCGCGGTAATTCGGCCTTTGGCAGATATAGGATGCGGAGGCCGGGTa                                                         >  1:1256310/1‑149 (MQ=255)
                                                                                            aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                       <  2:1084614/151‑1 (MQ=255)
                                                                                                ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                   <  1:1507245/151‑1 (MQ=255)
                                                                                                 cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCggt                                                  >  2:1128551/1‑149 (MQ=255)
                                                                                                 cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGtt                                                  <  2:1130119/151‑1 (MQ=255)
                                                                                                  gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGtt                                                  >  1:20252/1‑150 (MQ=255)
                                                                                                  gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCggt                                                  >  2:462776/1‑148 (MQ=255)
                                                                                                   cTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAACGCTGGATTGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTTGGCCTTAGCCCGATATAATATGCGTAGGCCGGATAAAGCggt                                                  >  2:314503/1‑147 (MQ=255)
                                                                                                    tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCggt                                                  >  2:528877/1‑146 (MQ=255)
                                                                                                    tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTGGCGGAGAAAGCGGAAGCCCCGCAGTAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                               >  1:474431/1‑151 (MQ=255)
                                                                                                     cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcg                                              <  1:212478/151‑1 (MQ=255)
                                                                                                      aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGGTGGGGGAGAAGGCGGGAGCCGCGCAGGAATTTCGGCTTAGCCCGATATAGTATGCGTAGGCCGGGGAAGGCGGTCAcgc                                             >  1:1473611/1‑151 (MQ=255)
                                                                                                       gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGGTGGATGGGGTGGGGAAGAAAGCGGGAGCCGCGCAGGAAATCGGCCTTAGCCAGATATAGTATGGGTAGGCCGGGTAAGGGGTTCAcgcc                                            >  2:1403016/1‑151 (MQ=255)
                                                                                                        ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccg                                           <  1:1133385/151‑1 (MQ=255)
                                                                                                         cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGGGAAGAAAGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGAAATAGTATGCGTAGGCCGGGTAAGGGGGTCAcgcc                                            >  2:713533/1‑149 (MQ=255)
                                                                                                          tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgccgcc                                         >  2:1278877/1‑150 (MQ=255)
                                                                                                           ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:932903/150‑1 (MQ=255)
                                                                                                           ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCa                                         <  1:409842/150‑1 (MQ=255)
                                                                                                           ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                        <  1:392974/151‑1 (MQ=255)
                                                                                                             cGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAAAAAATGCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTTGGCCTTAGCCCGATATAATATGCGTAGGCCGGGTAAGGCGGTCCCGCCGCATcc                                      >  2:924843/1‑151 (MQ=255)
                                                                                                              gACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGGTAAGGCGGTCACGCCGCATcc                                      >  1:463923/1‑150 (MQ=255)
                                                                                                               aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAAGCGGGAGCCGCGCAGTAAATCGGCCTTAGCCAGATATAGTATGCGTAAGCCGGGTAAAGCGGTCACGCCCCATCCg                                     >  1:1422244/1‑150 (MQ=255)
                                                                                                               aCCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                    <  2:1664984/151‑1 (MQ=255)
                                                                                                                 cATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACCGCTGGATGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCCTCCGGc                                   >  1:299033/1‑150 (MQ=255)
                                                                                                                  atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                 <  1:229579/151‑1 (MQ=255)
                                                                                                                    atTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATCCATCAGAAAATGCTGGGTGGGGTGGCGAAAAAAGCGGGAGCCGCGCCGCAATTCGGCCTTTGCCAGATATAATATTCGTAGGCCGGGTAAGGCGTTCACGCCGCCTCCCGCCttt                               >  1:1181258/1‑151 (MQ=255)
                                                                                                                    atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGAGTCAACACGAACAGCTGGGTGGGCGGGGGAAAAAAGAGGAAAACGCACAGTAAATTGGAGTTAGGCAAAGATAATATGCGTAAGACGGCTAAGGGGTTTTACCCTAATCCGGCAtt                                >  1:1511529/1‑150 (MQ=255)
                                                                                                                       cTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAAGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                               >  1:92922/1‑148 (MQ=255)
                                                                                                                            gatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATa                        <  1:3648/150‑1 (MQ=255)
                                                                                                                              tgaGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAAGCGGGAGCCGCGCCGTAATTTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCCGCATTTTAGCATAGTg                     >  1:848356/1‑151 (MQ=255)
                                                                                                                                  cAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTg                 <  2:1248738/151‑1 (MQ=255)
                                                                                                                                        cggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa           <  2:193513/151‑1 (MQ=255)
                                                                                                                                         ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc          >  2:741941/1‑151 (MQ=255)
                                                                                                                                          gTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg         >  2:779458/1‑151 (MQ=255)
                                                                                                                                            ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg         <  1:1494249/150‑1 (MQ=255)
                                                                                                                                              ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCt       <  2:1087408/150‑1 (MQ=255)
                                                                                                                                               tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCtt      <  1:1713962/150‑1 (MQ=255)
                                                                                                                                               tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg     >  2:588607/1‑151 (MQ=255)
                                                                                                                                                 gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt   >  1:1579822/1‑151 (MQ=255)
                                                                                                                                                 gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt   <  2:1343391/151‑1 (MQ=255)
                                                                                                                                                    gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATg  <  2:1691003/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456956‑4457251

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