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

AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCA  >  NC_000913/3329500‑3329987
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
aGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                  >  1:304228/1‑251 (MQ=255)
      ccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCACCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTCGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAAc                                                                                                                                                                                                                                             >  1:527254/1‑250 (MQ=255)
           gCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTATGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                                                                                                                       <  1:8955/251‑1 (MQ=255)
                       ggCTTTTATACGCGTGGCATCTTATCACCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                                            <  1:174740/250‑1 (MQ=255)
                           ttgataCGCGTGGCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACGGGATGGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGc                                                                                                                                                                                                                         <  2:413833/246‑1 (MQ=255)
                                    gTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCCTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGtgacacctcggggcgcaccaccgc                                                                                                                                                                                                              >  2:348734/1‑228 (MQ=255)
                                     gggCATCTTAGTCCCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGAGGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgctg                                                                                                                                                                                                             <  2:769123/250‑1 (MQ=255)
                                               ttaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCAACTTGAGCGGGAAACTGCCGCGGCAGaatc                                                                                                                                                                                                   >  2:32447/1‑248 (MQ=255)
                                                          cGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACCGGCGGGGCTCCCCTGGAGCGGGACACTGCTGCGCCACACTCACGGTAACgc                                                                                                                                                                                        >  2:613070/1‑249 (MQ=255)
                                                          cGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCAACTCAGGTTAACGa                                                                                                                                                                                        >  1:160272/1‑250 (MQ=255)
                                                           gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAc                                                                                                                                                                                      <  1:316467/251‑1 (MQ=255)
                                                                    ccAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAACTATGAGTTAAAACAGGCGGGGATCCCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCCGGAActg                                                                                                                                                                             >  2:331243/1‑249 (MQ=255)
                                                                        gTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAAAGAACCTGGAACgggagg                                                                                                                                                                         >  1:631657/1‑247 (MQ=255)
                                                                               cccccccccGTGGTTCTGCCGAAGGGCAATTCTGTGAAATGGATGTGGTGCCAGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCACTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGt                                                                                                                                                                  <  2:700114/247‑1 (MQ=255)
                                                                                    cgtcggggTTCTGCCGAAGGGCCATTCTGCGAACCGGATGTGGGGCCAGGCGACAGGATCCGCTTTACGCGGACGGGATGCCGGACACAACGTACTGAGTCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGGCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTATTTGCCAGTAAACa                                                                                                                                                             <  2:294062/245‑1 (MQ=255)
                                                                                       cGTGGTTCGGCCGAAGCGCAATACTGCGAACTGGAGGTGGTGCCAGGCGACCGGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTc                                                                                                                                                          <  2:361594/251‑1 (MQ=255)
                                                                                                gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCCGGTTAACGAACCTGGAACGGAAGTTGCCCGTAAACAGTCGGCCCCGcc                                                                                                                                                 >  1:486933/1‑250 (MQ=255)
                                                                                                          aaTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGATATCACCTGGAGCGGAACACTTCTGCGCCCCACTCACGGTAACGACCCTGGAACCGGAGTTTCCCGTAAACAATCCGcccccccgcccgacctg                                                                                                                                       >  2:206059/1‑238 (MQ=255)
                                                                                                             ccTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCtt                                                                                                                                    <  2:375828/250‑1 (MQ=255)
                                                                                                                gcgaacggtaggtggtcCCAGGCGGCCGGAACCGCTTTCCGCTGGCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGAGGGAGCCAGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCATTAAACAGTCGGCCCCGCTGCACGATCTGCTTaa                                                                                                                                 <  2:380229/239‑1 (MQ=255)
                                                                                                                 ggAACTGGATGTGGTGCCAGGCGACCTGACCCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGa                                                                                                                               <  2:333623/250‑1 (MQ=255)
                                                                                                                   aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAACTATGAGTTAAAACAGGCGGGTATCCCCTCGAGCGGAACACCGCTGCGCCCAACTCCAGGTAACGAACATGGAACGGTAAGTGCCCGTTAACAATCGCCCCCCCTCCACGACCTGCTTAAgttt                                                                                                                             >  2:109446/1‑248 (MQ=255)
                                                                                                                   aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAtt                                                                                                                             <  1:245129/251‑1 (MQ=255)
                                                                                                                    aCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTa                                                                                                                            <  1:425626/251‑1 (MQ=255)
                                                                                                                     cTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCCGGCGCGGAACACTCCTTCGCCAGACACCGGTTTACGAACCCGGGACCGTAGTTGCCAGTCAACAGTCCGCCCCCCTCCACGATCTGCCTTAGATTAt                                                                                                                           >  2:281082/1‑251 (MQ=255)
                                                                                                                      tGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCCGCACGATATGCTTAAAATTAtt                                                                                                                          >  1:465648/1‑250 (MQ=255)
                                                                                                                       ggttgtgGTGCCAGTCGCCCTGAACCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGc                                                                                                                         <  2:703081/248‑1 (MQ=255)
                                                                                                                        gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCt                                                                                                                        <  1:145828/251‑1 (MQ=255)
                                                                                                                        gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGGACAACGCTTCGCCCGAATTAGGGTTACGAACCTGGAACGGTAGTTGCCAGGAAACAGTCGGGCCCGCTGCACGATCTGCTTAAGGTTATGCt                                                                                                                        >  1:541488/1‑251 (MQ=255)
                                                                                                                           gtggtgCCATGCGGCCCGAACTGCTTTCCGCGGACGGGTTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTGGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTgca                                                                                                                     <  2:383258/251‑3 (MQ=255)
                                                                                                                           gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCCGTAAACAGTCGGCCCCGCTGCACGCTCTGCCTTAGATTATGCTGa                                                                                                                      >  1:769304/1‑250 (MQ=255)
                                                                                                                                 ccAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGACCCAGGGTAACGAACCTGGAACGGGGGTTGCCAGGAAACAGACGGCCCCGCGGCACGATCCGCTTTAGAATATGATGAAAAAgac                                                                                                               >  2:677476/1‑249 (MQ=255)
                                                                                                                                  cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGAACACCCGGGGCGGATCAACGCTGCGCCAGACTCAGGTTTACGAACCTGGAACGGTAAATGCCAGCCAACAGCGGGCCCCCCGTCACGATCTGCTTGAGATTTTGCTgttaaagtac                                                                                                              >  1:127883/1‑242 (MQ=255)
                                                                                                                                   aGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGAACACTGCTGCGCCAGACTCAGGTTAATGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCCCTGCACGATCTGCTTAAGATTTTGCTGGAAAAGTCgg                                                                                                             >  1:677241/1‑251 (MQ=255)
                                                                                                                                   aGGCGACCTGAACCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCgg                                                                                                             <  1:535826/251‑1 (MQ=255)
                                                                                                                                    ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGa                                                                                                            <  1:292075/251‑1 (MQ=255)
                                                                                                                                          cTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGCAAACTGCTGCGCCAGACTCAGGATAACGAACCTGGAAAGGTAGTTGCCAGTAAAACGTCGGCCCCGCTGCCCGATATGCTTTAGATTCTGCTGAAAAAATCGGACAACAt                                                                                                      >  1:726341/1‑251 (MQ=255)
                                                                                                                                           tGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAATGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTTAAAAAATCGGACAACATg                                                                                                     >  1:690719/1‑251 (MQ=255)
                                                                                                                                           ggAACCGCTTTGCGCGGAGGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCGATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAATGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATg                                                                                                     <  2:690719/250‑1 (MQ=255)
                                                                                                                                                 gCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATAACCAGGAGCGGGAAACTGCTGCGCCAGACTCAGGTTAACGAAACTGGAAAGGTAGTTGCCAGGAAACAGTAGGCCCCGGTGCACGCTTCGCTTTAAGTTATGCTGGAAAAGTCGGACAACATGATcccc                                                                                               >  1:481562/1‑248 (MQ=255)
                                                                                                                                                   tttATGCTGACGGGTGGCCCGCCACGACGTGCTGAGCCGCCCCCGTTGGCTGTTGCCGGGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCg                                                                                              <  2:91531/250‑1 (MQ=255)
                                                                                                                                                           gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGt                                                                                      <  1:361/249‑1 (MQ=255)
                                                                                                                                                              gggAAGCCGGCCACAAGGGTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGATCGCCGATACGGTTTTcc                                                                                 <  2:811888/251‑1 (MQ=255)
                                                                                                                                                                gaggcctgccACAACGGTCTGGGCCGCTCCCTGTGGATGTGGCCGTGCAGGAGTGGGCCAGCTATGCCTGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTTGAGCGGAACACTGCGTCGCCAGACTCAGGTTATCGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTTAAGATTACGCTGTAAAAGGCGGACAACTTGATCGCCGATACGGTTTTCCGc                                                                               <  2:508413/248‑1 (MQ=255)
                                                                                                                                                                            aaCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCCTGATAGGCCAt                                                                   >  1:740042/1‑251 (MQ=255)
                                                                                                                                                                             aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGGACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTTAAGATTATGCTGAAAAAATCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATg                                                                  >  1:402824/1‑251 (MQ=255)
                                                                                                                                                                             aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATg                                                                  <  1:405599/251‑1 (MQ=255)
                                                                                                                                                                                    gAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGGCAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCtt                                                           >  1:20096/1‑251 (MQ=255)
                                                                                                                                                                                            cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                                                    <  1:805292/250‑1 (MQ=255)
                                                                                                                                                                                                 tGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTgg                                               <  1:613070/250‑1 (MQ=255)
                                                                                                                                                                                                    cTTTTGCCGTGCAGGATGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCGGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGCGCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGag                                                                                            <  2:647225/203‑2 (MQ=255)
                                                                                                                                                                                                      tttGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACAt                                          <  1:109446/250‑1 (MQ=255)
                                                                                                                                                                                                        tGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTACTGCGCCAGACTCACGGTAACGAACCTGGAACAGTAATTGCCAGTAAACACTCGGCCCCGCTGCACCGTACGCTTAAGATTATGCTGCAAAAGCAGGGCAACATGATCGCCGATACGGATTTTCCCCCGAATGGCCCTGAGCCGTTCAATGTGCCTTGAACATGGc                                       >  2:45558/1‑251 (MQ=255)
                                                                                                                                                                                                               cAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCgg                                <  1:677476/251‑1 (MQ=255)
                                                                                                                                                                                                                       gAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGGTCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCgg                           >  1:6923/1‑249 (MQ=255)
                                                                                                                                                                                                                         gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTAt                                                                                                                           <  2:651981/152‑1 (MQ=255)
                                                                                                                                                                                                                          ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGTAACGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGATGCACGATCTGCTTAAGATTATACTTAAAAAGTCGGAAAACATGCTCGCCGATACGGGTTTCCGCATGATAGGCCCTGCGCGCTTCCATATGGCTGGGAAATGGCGGGCCCGGTGGGACGcc                      >  1:298529/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAATGGTAGTTGCCAGGAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGt                <  1:184636/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGTGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCa              <  1:45558/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTTCCAGATAACAGTCGGCCCCGCTGCACGATCTGCTTAAAGTTATGCTGAAACAGTCGGACAAAATGATCGCCGATACGGTTTTCCGCCTGATAGGCCATGCGCGCTTCAAATTGCCTGGAACATGGCTGGCCgggctggacccccggcgtcc              >  1:222393/1‑234 (MQ=255)
                                                                                                                                                                                                                                       tGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCGGACGCCCTGCCTCAGATcc         >  1:403368/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cAATTCTGAAATATGAGTTAAACCAGGCGGGTATCACCTGGAGCGGAACACTGCGGCGCCAGACTCAGGTTAACGACCCTGGAACGGTAGTTGCCAGTAAACATTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTGCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTg       <  2:486933/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ttCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTTGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTACACGATCTGCTTAAGATTATGCTGAAAAAGACGGACAACATGATCGCCGATACGGGTTTCCGCATGGTAGGCCATGCGCGCTTCAATGTGCCTGGGACATGGCGGGACGGGGCGGGAGCCGTGAGTCAGATACTgcgg    >  1:124064/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCa  <  1:686758/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCA  >  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: