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

Only 100 of 101 total aligned reads displayed.

ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456958‑4457255
                                                                                                                                                   |                                                                                                                                                        
atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTTGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATcc                                                                                                                                                       <  1:953159/151‑1 (MQ=255)
atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATcc                                                                                                                                                       <  2:670659/151‑1 (MQ=255)
   ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGTAACGGCGGTTTTGGTGGGATCCATc                                                                                                                                                    <  1:355747/151‑1 (MQ=255)
   ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATc                                                                                                                                                    >  2:142379/1‑151 (MQ=255)
    cGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTc                                                                                                                                                        >  1:310544/1‑146 (MQ=255)
        gcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGCGGGATCCATCAGaaa                                                                                                                                               <  1:745787/151‑1 (MQ=255)
         caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCATCCGaaa                                                                                                                                               >  1:1157780/1‑150 (MQ=255)
          aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGaaa                                                                                                                                               >  2:78080/1‑149 (MQ=255)
               gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAAATGCT‑gg                                                                                                                                       >  1:347674/1‑151 (MQ=255)
                cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGa                                                                                                                                      <  1:434460/151‑1 (MQ=255)
                cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCATCAGAAACTGCT‑ggg                                                                                                                                      >  2:945753/1‑150 (MQ=255)
                 agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCT‑gg                                                                                                                                       >  2:906913/1‑149 (MQ=255)
                 agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGATCCATCAGAAACTGCT‑gg                                                                                                                                       >  1:902458/1‑149 (MQ=255)
                 agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCT‑gg                                                                                                                                       >  1:674762/1‑149 (MQ=255)
                 agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGAt                                                                                                                                     <  2:715666/151‑1 (MQ=255)
                 agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGAt                                                                                                                                     <  1:625266/151‑1 (MQ=255)
                    agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAAATGCT‑ggggggg                                                                                                                                  >  2:771052/1‑146 (MQ=255)
                     gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCCGAAACTGGt‑gggtggg                                                                                                                                  >  1:1194629/1‑150 (MQ=255)
                       ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTg                                                                                                                               <  2:1201825/151‑1 (MQ=255)
                         cagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCACCCGAAAATGCt‑gggtggggttg                                                                                                                              >  1:1070984/1‑146 (MQ=255)
                           gaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATCCATCAGAAAATGGT‑GGATGGGCTgg                                                                                                                              >  1:1067956/1‑148 (MQ=255)
                              gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGGCATCAAAAACTGCt‑gggtggggcgg                                                                                                                              >  2:958179/1‑141 (MQ=255)
                               ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCT‑GGGTGGGCTGGCgaaga                                                                                                                        >  1:506492/1‑150 (MQ=255)
                                       cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGATCCATCAGAAAAAGCT‑GGGTGGGGTGGCGAAGAAGGCgg                                                                                                                  >  1:892219/1‑148 (MQ=255)
                                       cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGaa                                                                                                                <  1:716510/150‑1 (MQ=255)
                                       cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAg                                                                                                               <  1:60504/151‑1 (MQ=255)
                                         gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                             <  2:1055577/151‑1 (MQ=255)
                                         gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACAGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGcc                                                                                                             <  1:687446/151‑1 (MQ=255)
                                           cTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCT‑GGGTGGGGTGGGGAAGAAAGCgggaacc                                                                                                             >  1:691138/1‑144 (MQ=255)
                                              cacaCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCa                                                                                                       <  2:680250/151‑1 (MQ=255)
                                               acacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATTCT‑GGGTGGGCTGGCGAAGAAGGCGGGAGgcgc                                                                                                           >  2:783128/1‑147 (MQ=255)
                                                cacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGt                                                                                                     <  2:340600/151‑1 (MQ=255)
                                                 acTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTa                                                                                                    <  2:1078265/151‑1 (MQ=255)
                                                  cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCT‑GGATGGGGTGGCGGAAAAAGCGGGAGCcgc‑g                                                                                                         >  1:1078509/1‑145 (MQ=255)
                                                   tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGTGGGAGCCAACCGAAAATGCT‑GGGTGGGGTGGCGGAGAAAGCGGAAGCCGC‑GCAGGaa                                                                                                   >  2:463611/1‑150 (MQ=255)
                                                    ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAAtt                                                                                                 <  2:756143/151‑1 (MQ=255)
                                                         ggCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCT‑GGATGGGGTGGCGAAGAAGGCGGGAGCCGC‑GCAGGAATTCGGcc                                                                                            >  1:852254/1‑151 (MQ=255)
                                                          gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCt                                                                                           <  1:492316/151‑1 (MQ=255)
                                                          gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGACCACCAGAAACTGCT‑CGATGGGCTGGGCCGGAAGGGGGAAGACGC‑GCAGTAAATCCGccc                                                                                           >  2:192515/1‑150 (MQ=255)
                                                            gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTa                                                                                         <  1:1002093/151‑1 (MQ=255)
                                                            gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTa                                                                                         <  1:297658/151‑1 (MQ=255)
                                                             aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCT‑GGGTGGGCTGGCGGAGAAGGCGGGAGCCGC‑GCAGTAATTCGGCCtt                                                                                          >  2:7947/1‑149 (MQ=255)
                                                               gtgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGc                                                                                       <  1:1107517/150‑1 (MQ=255)
                                                                tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCT‑GGGTGGGCTGGCGGAGAAGGCGGGAGCCGC‑GCAGTAATTCGGCCTTAGcc                                                                                      >  1:961296/1‑150 (MQ=255)
                                                                tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGGAGAAGGCGGGAGCCGC‑GCCGTAATTCGGCCTTAGcc                                                                                      >  1:1168698/1‑150 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCTGAAAATGGT‑GGTTGGGCTGGGGCACAAGGGGGAAGCCCC‑GCCGTAATTCCGCCTTTGCCCGat                                                                                  >  2:233458/1‑151 (MQ=255)
                                                                    tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTTCCTCAAAAAATCCTGGGGTGGGCGGCGGAAGAAGG‑GGAAAGCGC‑GCCGTaaatttggcttatt                                                                                       >  2:931628/1‑133 (MQ=255)
                                                                    tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCT‑GGGTGGGGTGGGGGAGAAAGCGGGAGGCGC‑GCAGGAAATCGGCCTTAGCCAGat                                                                                  >  2:868917/1‑150 (MQ=255)
                                                                     aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAAATGCT‑GGGTGGGCTGGGGGAGAAAGCGGGAGCCGC‑GCAGTAATTCGGCCTTAGcc                                                                                      >  1:133696/1‑145 (MQ=255)
                                                                         tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCACCAGAAAATGCT‑GGATGGGCTGGCGAAGAAGGGGGAAGCCGC‑GCCGTAATTCGGCCTTAGCCAGatat                                                                                >  1:590464/1‑147 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGt                                                                             <  2:514210/149‑1 (MQ=255)
                                                                             aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGGAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTAt                                                                           >  1:152596/1‑148 (MQ=255)
                                                                               atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGGAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCg                                                                        >  1:339925/1‑149 (MQ=255)
                                                                               atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCG‑GGATGGGCTGGGGTAGAAGGCGGGAGCCGC‑CCAGTAATTTCGCCTTATCCAGATATAGTATGCGt                                                                       >  1:1201091/1‑150 (MQ=255)
                                                                                 atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGACCATCAGAAAATGCT‑GGATGGGCTGGCGAAGAAGGGGGAAGCCGC‑GCAGTAATTTGGCCCTAGCCAGATATAGTATGCGCAgg                                                                    >  1:1049975/1‑151 (MQ=255)
                                                                                 atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGTGGGAGCCATCAGAAACTGCT‑GGGTGGGGTGGCGGAGAAAGCGGAAGCCGC‑CCCGGAATTCGGCCTTTGCCAGGTATAGTATGCGTAgg                                                                    >  1:626287/1‑151 (MQ=255)
                                                                                     tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGGAGCCGC‑GCAGTAAATTGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                                  >  1:1172283/1‑149 (MQ=255)
                                                                                      ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGa                                                               <  1:544310/151‑1 (MQ=255)
                                                                                       gtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGAt                                                              <  1:412838/151‑1 (MQ=255)
                                                                                         gAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                             <  1:1105671/150‑1 (MQ=255)
                                                                                          aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCCTCAGAAAATGCT‑GGGTGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATa                                                             >  1:438872/1‑149 (MQ=255)
                                                                                          aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGTGGGGTGCAGCCGAAACTGCT‑TGATGGGGTTGCGGAAAAGGCGGGAAGCGC‑GCAGGAATTTGGCCTTAATCCCAGATAGCATGCGGTGGCCg                                                                 >  1:965532/1‑145 (MQ=255)
                                                                                           aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCT‑GGGTGGGGTGGCGAAGAAGGCGGGAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAgg                                                          >  1:985674/1‑151 (MQ=255)
                                                                                           aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCT‑GGGTGGGCTGGCGAAGAAGGGGGGAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAgg                                                          >  1:849363/1‑151 (MQ=255)
                                                                                              ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGt                                                       <  2:1194611/151‑1 (MQ=255)
                                                                                                gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACCGCT‑GGATGGGCTGGCGAAGAAAGCGGAAGCCGC‑GCAGGAATTCGGCCCTAGCCAGATATAGTATTCGGAGGCCGGGTAAGGCggt                                                      >  1:1181780/1‑148 (MQ=255)
                                                                                                gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAAAGGCGGTTTTTGGGGGATACCGCATAAAATGCT‑GGGTGGGCCGGCCGAGAAGGCGGAAGGCGCGGCGGCATTGCGGCC‑TAGGCAGATAGAGGATTGGTAGGCCGGACaacgggg                                                       >  2:895369/1‑144 (MQ=255)
                                                                                                  tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCT‑GGGTGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGTTCa                                                    >  1:249441/1‑150 (MQ=255)
                                                                                                  tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCT‑GGGTGGGCTGGCGAAAAAGGCGGAAGCCGC‑GCAGAAATTCGGCCTTAGCCAGAAATAGTATGGGTAGGCCGGAAAAGGGGTTccc                                                   >  1:838072/1‑149 (MQ=255)
                                                                                                   cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCa                                                    <  2:935715/149‑1 (MQ=255)
                                                                                                   cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                                   <  2:691135/150‑1 (MQ=255)
                                                                                                    aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCT‑GGATGGGGTGGCGGAGGAGGCGGAAGCCGC‑GCAGGAATTCGGCCTTAGCCAGATATAGGATGCGGAGGCCGGATAAAGCGGTCAcgc                                                 >  2:765292/1‑151 (MQ=255)
                                                                                                    aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgc                                                 <  2:915911/151‑1 (MQ=255)
                                                                                                         ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcgccgc                                              <  1:490725/149‑1 (MQ=255)
                                                                                                          gCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                           <  1:526508/151‑1 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGGT‑GGTTGGGCTCGCGGAAAAGGCCGAAGCCGC‑GGAGGAATTTTGCCTTAAGCAGGGATAGTATGCGGTGGCCCGCAAAGGCGTTCACGCCGCac                                            >  1:447099/1‑148 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAAGCGGAAGCCGC‑GCAGGAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCGGATAAGGGGGTTACGCCGCATcc                                          >  1:1136349/1‑151 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGCTTTTGTGGGAACCATCAGAAACTGCT‑GGGTGGGCCGGGGCAGAAGGCGGAAGACGC‑GGAATAAATCGGCCTTAGCCAGATATAATATGCGGAGGCCGGATAAGGGGTTCACGCCGCATcc                                          >  1:454512/1‑151 (MQ=255)
                                                                                                              ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                       <  2:720065/151‑1 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAAAAAATGCT‑GGATGGGCTGGCGAAGAAGGCGGGAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCAt                                     >  1:4164/1‑151 (MQ=255)
                                                                                                                     cTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGGGGAGAAGGGGGGAGGCGC‑GCCGGAAATTGGCCTTAGCCAGAAATAGGATGCGTAGGCCGGAAAAGGCGGTCACCCCGCATCCGGCAttt                                   >  2:1040847/1‑148 (MQ=255)
                                                                                                                     cTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGGAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCAttt                                   >  1:987151/1‑148 (MQ=255)
                                                                                                                         tgatgaGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGGAGCCGC‑GCAGGAATTCGGCCTTTGCCAGATATAGTAAGCGTAAGCCGGATAAGGCGTTCACGCCCCCTCCCGCCTTTGAg                                >  2:892318/1‑147 (MQ=255)
                                                                                                                               gCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCt                      <  2:532955/151‑1 (MQ=255)
                                                                                                                                cAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTg                     <  2:363265/151‑1 (MQ=255)
                                                                                                                                   cggcggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATg                  <  1:954825/151‑1 (MQ=255)
                                                                                                                                    ggcggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGc                 <  1:831013/151‑1 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg                <  1:381455/151‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa               <  2:736678/150‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa               <  2:1093558/150‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc              >  2:923253/1‑151 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc              >  1:928985/1‑151 (MQ=255)
                                                                                                                                            tggtggGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCtt          <  2:674760/150‑1 (MQ=255)
                                                                                                                                              gtggGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt       >  1:588733/1‑151 (MQ=255)
                                                                                                                                               tggGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAt       <  2:567441/150‑1 (MQ=255)
                                                                                                                                               tggGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATg      >  2:1051320/1‑151 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc     <  2:347672/151‑1 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc     >  1:815796/1‑151 (MQ=255)
                                                                                                                                                  gATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGt   >  2:614291/1‑151 (MQ=255)
                                                                                                                                                   aTCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTc  >  1:99614/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCT‑GGATGGGCTGGCGAAGAAGGCGGAAGCCGC‑GCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTC  >  NC_000913/4456958‑4457255

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: ‑