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

CCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456961‑4457251
                                                                                                                                                |                                                                                                                                                  
ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAt                                                                                                                                               <  2:967800/150‑1 (MQ=255)
  gtggcagggggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGGGATGAGCAACGGCGGGTTTTGTGGGGTcc                                                                                                                                                 >  1:282565/1‑146 (MQ=255)
   tggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGAtgc                                                                                                                                                 >  2:884894/1‑143 (MQ=255)
      caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAc                                                                                                                                        <  2:989465/151‑1 (MQ=255)
      caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGGGGGGTCCATCAGaaa                                                                                                                                         >  1:255855/1‑150 (MQ=255)
                aagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCCGAAAATGCTGGGTgg                                                                                                                              >  2:82162/1‑151 (MQ=255)
                  gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATggg                                                                                                                             <  1:1075287/150‑1 (MQ=255)
                  gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGc                                                                                                                            <  1:987142/151‑1 (MQ=255)
                  gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGGGGGGTCCATCAGGAAATGCtgggtggg                                                                                                                             >  2:759254/1‑150 (MQ=255)
                    ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCt                                                                                                                           <  1:651189/150‑1 (MQ=255)
                       agaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCg                                                                                                                       <  2:1112645/151‑1 (MQ=255)
                        gaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAAATGCTGggtggggtgggg                                                                                                                       >  1:726707/1‑148 (MQ=255)
                        gaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGCTgg                                                                                                                         >  2:172300/1‑148 (MQ=255)
                          agCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaag                                                                                                                    <  2:472086/151‑1 (MQ=255)
                                  ttCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGGGGAGAAGGCgg                                                                                                             >  2:137230/1‑150 (MQ=255)
                                  ttCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCATCAGAAAATGCTGGATGGGGTGGCggagaaagggg                                                                                                             >  2:242515/1‑145 (MQ=255)
                                    cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa                                                                                                           <  1:1069508/150‑1 (MQ=255)
                                     aGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAg                                                                                                          <  2:66201/150‑1 (MQ=255)
                                      gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                        <  2:1159686/151‑1 (MQ=255)
                                      gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                        <  1:87919/151‑1 (MQ=255)
                                       ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcg                                                                                                       <  1:1215226/151‑1 (MQ=255)
                                          gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                      <  2:283345/149‑1 (MQ=255)
                                            acacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTTGGGTCCCTCAGAAACTGCCGGAGGGGGTGGCGTAGAAAGCGGGAGGcg                                                                                                       >  2:1062290/1‑146 (MQ=255)
                                                 ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTaa                                                                                               <  2:73292/149‑1 (MQ=255)
                                                 ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTaa                                                                                               <  2:23297/149‑1 (MQ=255)
                                                 ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAtt                                                                                             <  1:1017711/151‑1 (MQ=255)
                                                 ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAt                                                                                              <  2:1057720/150‑1 (MQ=255)
                                                     tGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTTCCTCCAAAACCTCCGGGTGGGCTGGCGGAGAAGGGGGGAGCCgcgccggaaat                                                                                             >  1:1106351/1‑140 (MQ=255)
                                                      ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGGGGAGAAAGCGGGAGCCGCGCAGTAAATTGGcc                                                                                        >  1:546291/1‑151 (MQ=255)
                                                      ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTTGGcc                                                                                        >  1:25269/1‑151 (MQ=255)
                                                            gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATCCAAAAATGCTGGATGGGGTGGGGGAAAAAGCGGGAGGCGCGCAGGAATTTGGCCTTAggc                                                                                  >  1:950985/1‑149 (MQ=255)
                                                             tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGcc                                                                                  >  2:1279010/1‑150 (MQ=255)
                                                               gATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAAAGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGa                                                                               <  1:1109825/151‑1 (MQ=255)
                                                                aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGat                                                                              <  2:982062/151‑1 (MQ=255)
                                                                  aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCCTCAGAAACTGCTGGGTGGGGTGGCGAAGAAGGGGGGAGCCGCGCAGGAAATTGGCCCTAggc                                                                                  >  1:1111161/1‑143 (MQ=255)
                                                                  aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCCTCAGAAACTGCTGGGTGGGGTGGGGGAGAAAGGGGGAGCCGCGCCGTAAATTGGCCCTAGCCCGatat                                                                            >  1:756942/1‑151 (MQ=255)
                                                                    gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGatat                                                                            >  2:373604/1‑149 (MQ=255)
                                                                       ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                       <  2:1132596/151‑1 (MQ=255)
                                                                       ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                       <  1:150444/151‑1 (MQ=255)
                                                                       ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGGCCCGCAGAAAATGCCGGGGGGGGTGGCGGAGAAAGGAGAGGCGCCGCCGTAAATCGGCCTTTGccccatttaat                                                                         >  1:467825/1‑140 (MQ=255)
                                                                        gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg                                                                      <  1:356337/151‑1 (MQ=255)
                                                                          aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAAATCGGCCTTAGCCAGATATAGTATGCg                                                                    >  2:160336/1‑151 (MQ=255)
                                                                               tggtggTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGTAGGGGGCCCACACAAACTGTCGGTTGGGCTTGCGCAGAAGGGGGGATACGCGCCGGAATTCGGCCCTACCACCTTATAGTATGCGTAGGc                                                               >  1:627560/1‑151 (MQ=255)
                                                                                ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGAGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                              >  2:360660/1‑151 (MQ=255)
                                                                                ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGc                                                               <  1:361240/150‑1 (MQ=255)
                                                                                ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGc                                                               <  2:1124018/150‑1 (MQ=255)
                                                                                  tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTTGGCCTTTGCCCGATATAGTATGCGGAGGGCgg                                                            >  2:449682/1‑151 (MQ=255)
                                                                                   ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCgg                                                            >  1:387443/1‑150 (MQ=255)
                                                                                    gtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGAt                                                          <  1:1062281/151‑1 (MQ=255)
                                                                                     tgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAACTTCTGGGTGGGCTGGGGGAGAAAGCGGGAGGCGCGCAGGAATTCGGCCTTAGGCCGATTAAGTATGCGGAGGCCgg                                                            >  1:250955/1‑148 (MQ=255)
                                                                                      gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                         >  2:1357754/1‑150 (MQ=255)
                                                                                      gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                         <  2:871997/150‑1 (MQ=255)
                                                                                      gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                         <  2:177466/150‑1 (MQ=255)
                                                                                       aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                         >  2:104046/1‑149 (MQ=255)
                                                                                        aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAgg                                                      <  2:1127412/151‑1 (MQ=255)
                                                                                         aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                     <  2:1247766/151‑1 (MQ=255)
                                                                                         aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGc                                                     <  1:1307225/151‑1 (MQ=255)
                                                                                             gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCg                                                    >  1:311098/1‑148 (MQ=255)
                                                                                             gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                 <  2:1111170/151‑1 (MQ=255)
                                                                                              cTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTTCCTCCGAAACTTCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCCCGCCGGAATTTGGGCCTTGCCCGAGATAGTATGGGGCGGCCgggaaaagcg                                                    >  1:19504/1‑139 (MQ=255)
                                                                                                cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                               <  2:145643/150‑1 (MQ=255)
                                                                                                cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                               <  2:43790/150‑1 (MQ=255)
                                                                                                 aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgc                                             >  2:637750/1‑151 (MQ=255)
                                                                                                     tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGGAGAAGGCGGAAGCCGCGCAGTAATTTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGGGTTCAcgccgc                                          >  2:923826/1‑150 (MQ=255)
                                                                                                      ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCCTCAAAAACTGCTGGGTGGGGTTGCGAAGAAGGCGGGAGCCGCGCAGGAAATTGGCCCTAGCCAGAGATAGTATGCGTAGGCCGGGGAAGGCGTTCaaacc                                            >  1:1125019/1‑143 (MQ=255)
                                                                                                        cGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGCGGAAGGCGCGCAGTAATTTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATcc                                      >  2:972746/1‑151 (MQ=255)
                                                                                                         gACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATTCGTAGGCCGGGTAAGGCGTTCACGCCGCATcc                                      >  1:1095694/1‑150 (MQ=255)
                                                                                                             atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGCGGAGAAGGCGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGc                                   >  1:1160534/1‑149 (MQ=255)
                                                                                                              tatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGGAGGCCGGGTAAAGCGGTCACGCCGCATCCGGCAtt                                >  2:935610/1‑151 (MQ=255)
                                                                                                               atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                               <  2:534803/151‑1 (MQ=255)
                                                                                                                 tCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTg                              <  2:448254/150‑1 (MQ=255)
                                                                                                                    ggTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACCGCTGGATGGGGTGGCGGAGAAAGGGGAAGCCGCGCCGTAATTTGGCCCTTGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAttt                               >  1:960666/1‑146 (MQ=37)
                                                                                                                      tgatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCAt                         <  1:1142219/150‑1 (MQ=255)
                                                                                                                        atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGt                      <  1:475132/151‑1 (MQ=255)
                                                                                                                         tgaGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGGGGAGAAGGGGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTTGAGCATa                        >  1:958593/1‑148 (MQ=38)
                                                                                                                         tgaGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACAGCTGGATGGGGTGGCGGAGAAAGCGGGAGGCGCCCCGGAATTTGGCCTTAGCCAGATATAGTATGCGGAGGCCGGGTAAGGGGTTCACGCCCCATCCGGCCttt                               >  1:1104308/1‑141 (MQ=255)
                                                                                                                              aaCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTg                 <  1:360894/150‑1 (MQ=255)
                                                                                                                               aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt               <  2:1243994/151‑1 (MQ=255)
                                                                                                                               aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCGGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt               <  2:1272954/151‑1 (MQ=255)
                                                                                                                                cggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt               <  1:759247/150‑1 (MQ=255)
                                                                                                                                  gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            <  1:971731/151‑1 (MQ=255)
                                                                                                                                  gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            <  2:1287126/151‑1 (MQ=255)
                                                                                                                                  gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg            <  1:362739/151‑1 (MQ=255)
                                                                                                                                    ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc          <  1:154878/151‑1 (MQ=255)
                                                                                                                                     gTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg         <  1:57477/151‑1 (MQ=255)
                                                                                                                                            tggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATg  >  1:1303999/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATG  >  NC_000913/4456961‑4457251

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: