Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,329,743 G→T D261Y (GAT→TAT)  dacB → peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,329,7430GT98.4% 199.0 / NA 64D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (30/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCA  >  NC_000913/3329500‑3329987
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
aGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                  >  1:304228/1‑251 (MQ=255)
      ccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCACCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTCGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAAc                                                                                                                                                                                                                                             >  1:527254/1‑250 (MQ=255)
           gCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTATGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                                                                                                                       <  1:8955/251‑1 (MQ=255)
                       ggCTTTTATACGCGTGGCATCTTATCACCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                                            <  1:174740/250‑1 (MQ=255)
                           ttgataCGCGTGGCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACGGGATGGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGc                                                                                                                                                                                                                         <  2:413833/246‑1 (MQ=255)
                                    gTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCCTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGtgacacctcggggcgcaccaccgc                                                                                                                                                                                                              >  2:348734/1‑228 (MQ=255)
                                     gggCATCTTAGTCCCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGAGGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgctg                                                                                                                                                                                                             <  2:769123/250‑1 (MQ=255)
                                               ttaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCAACTTGAGCGGGAAACTGCCGCGGCAGaatc                                                                                                                                                                                                   >  2:32447/1‑248 (MQ=255)
                                                          cGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACCGGCGGGGCTCCCCTGGAGCGGGACACTGCTGCGCCACACTC‑ACGGTAACgc                                                                                                                                                                                        >  2:613070/1‑249 (MQ=255)
                                                          cGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCAACTC‑AGGTTAACGa                                                                                                                                                                                        >  1:160272/1‑250 (MQ=255)
                                                           gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAAc                                                                                                                                                                                      <  1:316467/251‑1 (MQ=255)
                                                                    ccAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAACTATGAGTTAAAACAGGCGGGGATCCCCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCCGGAActg                                                                                                                                                                             >  2:331243/1‑249 (MQ=255)
                                                                        gTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGGTAAAGAACCTGGAACgggagg                                                                                                                                                                         >  1:631657/1‑247 (MQ=255)
                                                                               cccccccccGTGGTTCTGCCGAAGGGCAATTCTGTGAAATGGATGTGGTGCCAGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCACTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGt                                                                                                                                                                  <  2:700114/247‑1 (MQ=255)
                                                                                    cgtcggggTTCTGCCGAAGGGCCATTCTGCGAACCGGATGTGGGGCCAGGCGACAGGATCCGCTTTACGCGGACGGGATGCCGGACACAACGTACTGAGTCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGGCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTATTTGCCAGTAAACa                                                                                                                                                             <  2:294062/245‑1 (MQ=255)
                                                                                       cGTGGTTCGGCCGAAGCGCAATACTGCGAACTGGAGGTGGTGCCAGGCGACCGGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTc                                                                                                                                                          <  2:361594/251‑1 (MQ=255)
                                                                                                gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑CGGTTAACGAACCTGGAACGGAAGTTGCCCGTAAACAGTCGGCCCCGcc                                                                                                                                                 >  1:486933/1‑250 (MQ=255)
                                                                                                          aaTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGATATCACCTGGAGCGGAACACTTCTGCGCCCCACTC‑ACGGTAACGACCCTGGAACCGGAGTTTCCCGTAAACAATCCGcccccccgcccgacctg                                                                                                                                       >  2:206059/1‑238 (MQ=255)
                                                                                                             ccTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCtt                                                                                                                                    <  2:375828/250‑1 (MQ=255)
                                                                                                                gcgaacggtaggtggtcCCAGGCGGCCGGAACCGCTTTCCGCTGGCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑GGGAGCCAGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCATTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑aa                                                                                                                                 <  2:380229/239‑1 (MQ=255)
                                                                                                                 ggAACTGGATGTGGTGCCAGGCGACCTGACCCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGa                                                                                                                               <  2:333623/250‑1 (MQ=255)
                                                                                                                   aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAACTATGAGTTAAAACAGGCGGGTATCCCCTCGAGCGGAACACCGCTGCGCCCAACTCCAGGT‑AACGAACATGGAACGGTAAGTGCCCGTTAACAATCGCCCCCCCTCCACGACCTGCTT‑AAgttt                                                                                                                             >  2:109446/1‑248 (MQ=255)
                                                                                                                   aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGAtt                                                                                                                             <  1:245129/251‑1 (MQ=255)
                                                                                                                    aCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTa                                                                                                                            <  1:425626/251‑1 (MQ=255)
                                                                                                                     cTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCCGGCGCGGAACACTCCTTCGCCAGACAC‑CGGTTTACGAACCCGGGACCGTAGTTGCCAGTCAACAGTCCGCCCCCCTCCACGATCTGCCT‑TAGATTAt                                                                                                                           >  2:281082/1‑251 (MQ=255)
                                                                                                                      tGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCCGCACGATATGCTT‑AAAATTAtt                                                                                                                          >  1:465648/1‑250 (MQ=255)
                                                                                                                       ggttgtgGTGCCAGTCGCCCTGAACCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGc                                                                                                                         <  2:703081/248‑1 (MQ=255)
                                                                                                                        gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCt                                                                                                                        <  1:145828/251‑1 (MQ=255)
                                                                                                                        gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGGACAACGCTTCGCCCGAATT‑AGGGTTACGAACCTGGAACGGTAGTTGCCAGGAAACAGTCGGGCCCGCTGCACGATCTGCTT‑AAGGTTATGCt                                                                                                                        >  1:541488/1‑251 (MQ=255)
                                                                                                                           gtggtgCCATGCGGCCCGAACTGCTTTCCGCGGACGGGTTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTGGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTgca                                                                                                                     <  2:383258/251‑3 (MQ=255)
                                                                                                                           gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGGAGTTGCCCGTAAACAGTCGGCCCCGCTGCACGCTCTGCCT‑TAGATTATGCTGa                                                                                                                      >  1:769304/1‑250 (MQ=255)
                                                                                                                                 ccAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGACCC‑AGGGTAACGAACCTGGAACGGGGGTTGCCAGGAAACAGACGGCCCCGCGGCACGATCCGCTT‑TAGAATATGATGAAAAAgac                                                                                                               >  2:677476/1‑249 (MQ=255)
                                                                                                                                  cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGAACACCCGGGGCGGATCAACGCTGCGCCAGACTC‑AGGTTTACGAACCTGGAACGGTAAATGCCAGCCAACAGCGGGCCCCCCGTCACGATCTGCTT‑GAGATTTTGCTgttaaagtac                                                                                                              >  1:127883/1‑242 (MQ=255)
                                                                                                                                   aGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGAACACTGCTGCGCCAGACTC‑AGGTTAATGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCCCTGCACGATCTGCTT‑AAGATTTTGCTGGAAAAGTCgg                                                                                                             >  1:677241/1‑251 (MQ=255)
                                                                                                                                   aGGCGACCTGAACCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCgg                                                                                                             <  1:535826/251‑1 (MQ=255)
                                                                                                                                    ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGa                                                                                                            <  1:292075/251‑1 (MQ=255)
                                                                                                                                          cTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGCAAACTGCTGCGCCAGACTC‑AGGATAACGAACCTGGAAAGGTAGTTGCCAGTAAAACGTCGGCCCCGCTGCCCGATATGCTT‑TAGATTCTGCTGAAAAAATCGGACAACAt                                                                                                      >  1:726341/1‑251 (MQ=255)
                                                                                                                                           tGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGGTAACGAACCTGGAATGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTTAAAAAATCGGACAACATg                                                                                                     >  1:690719/1‑251 (MQ=255)
                                                                                                                                           ggAACCGCTTTGCGCGGAGGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCGATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAATGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATg                                                                                                     <  2:690719/250‑1 (MQ=255)
                                                                                                                                                 gCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATAACCAGGAGCGGGAAACTGCTGCGCCAGACTC‑AGGTTAACGAAACTGGAAAGGTAGTTGCCAGGAAACAGTAGGCCCCGGTGCACGCTTCGCTTTAAG‑TTATGCTGGAAAAGTCGGACAACATGATcccc                                                                                               >  1:481562/1‑248 (MQ=255)
                                                                                                                                                   tttATGCTGACGGGTGGCCCGCCACGACGTGCTGAGCCGCCCCCGTTGGCTGTTGCCGGGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCg                                                                                              <  2:91531/250‑1 (MQ=255)
                                                                                                                                                           gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑t                                                                                      <  1:361/249‑1 (MQ=255)
                                                                                                                                                              gggAAGCCGGCCACAAGGGTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGGCGGACAACATGATCGCCGATACGG‑TTTTcc                                                                                 <  2:811888/251‑1 (MQ=255)
                                                                                                                                                                gaggcctgccACAACGGTCTGGGCCGCTCCCTGTGGATGTGGCCGTGCAGGAGTGG‑GCCAGCTATGCCTGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTTGAGCGGAACACTGCGTCGCCAGACTC‑AGGTTATCGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTT‑AAGATTACGCTGTAAAAGGCGGACAACTTGATCGCCGATACGG‑TTTTCCGc                                                                               <  2:508413/248‑1 (MQ=255)
                                                                                                                                                                            aaCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCCTGATAGGCCAt                                                                   >  1:740042/1‑251 (MQ=255)
                                                                                                                                                                             aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGGACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTT‑AAGATTATGCTGAAAAAATCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATg                                                                  >  1:402824/1‑251 (MQ=255)
                                                                                                                                                                             aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATg                                                                  <  1:405599/251‑1 (MQ=255)
                                                                                                                                                                                    gAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGGCAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCtt                                                           >  1:20096/1‑251 (MQ=255)
                                                                                                                                                                                            cccGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                                                    <  1:805292/250‑1 (MQ=255)
                                                                                                                                                                                                 tGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTgg                                               <  1:613070/250‑1 (MQ=255)
                                                                                                                                                                                                    cTTTTGCCGTGCAGGA‑TGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCGGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGCGCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGag                                                                                            <  2:647225/203‑2 (MQ=255)
                                                                                                                                                                                                      tttGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACAt                                          <  1:109446/250‑1 (MQ=255)
                                                                                                                                                                                                        tGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTACTGCGCCAGACTC‑ACGGTAACGAACCTGGAACAGTAATTGCCAGTAAACACTCGGCCCCGCTGCACCGTACGCTT‑AAGATTATGCTGCAAAAGCAGGGCAACATGATCGCCGATACGGATTTTCC‑CCCGAATGGCCCTGAGCCGTTCAATGTGCCTTGAACATGGc                                       >  2:45558/1‑251 (MQ=255)
                                                                                                                                                                                                               cAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCgg                                <  1:677476/251‑1 (MQ=255)
                                                                                                                                                                                                                       gAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGGTCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCgg                           >  1:6923/1‑249 (MQ=255)
                                                                                                                                                                                                                         gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTAt                                                                                                                           <  2:651981/152‑1 (MQ=255)
                                                                                                                                                                                                                          ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGTAACGGAACACTGCTGCGCCAGACTC‑AGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGATGCACGATCTGCTT‑AAGATTATACTTAAAAAGTCGGAAAACATGCTCGCCGATACGG‑GTTTCCGCATGATAGGCCCTGCGCGCTTCCATATGGCTGGGAAATGGCGGGCCCGGTGGGACGcc                      >  1:298529/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAATGGTAGTTGCCAGGAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGt                <  1:184636/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGTGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCa              <  1:45558/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGGAGTTTCCAGATAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAAGTTATGCTGAAACAGTCGGACAAAATGATCGCCGATACGG‑TTTTCCGCCTGATAGGCCATGCGCGCTTCAAATTGCCTGGAACATGGCTGGCCgggctggacccccggcgtcc              >  1:222393/1‑234 (MQ=255)
                                                                                                                                                                                                                                       tGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCGGACGCCCTGCCTCAGATcc         >  1:403368/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cAATTCTGAAATATGAGTTAAACCAGGCGGGTATCACCTGGAGCGGAACACTGCGGCGCCAGACTC‑AGGTTAACGACCCTGGAACGGTAGTTGCCAGTAAACATTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTGCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTg       <  2:486933/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ttCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTTGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTACACGATCTGCTT‑AAGATTATGCTGAAAAAGACGGACAACATGATCGCCGATACGG‑GTTTCCGCATGGTAGGCCATGCGCGCTTCAATGTGCCTGGGACATGGCGGGACGGGGCGGGAGCCGTGAGTCAGATACTgcgg    >  1:124064/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCa  <  1:686758/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGA‑TGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTT‑AAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGG‑TTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCA  >  NC_000913/3329500‑3329987

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