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,7430GT97.7% 146.5 / ‑3.2 44D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (22/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCA  >  NC_000913/3329500‑3329987
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
aGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                >  1:423313/1‑251 (MQ=255)
           gCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCCCCCGCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGAGCTGAACCGCTTTACGCTGACTGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTTGCTTTTGCCGTGCAGGATGTAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                                                                                                                     <  1:495962/251‑1 (MQ=255)
            ccTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGCTGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCggg                                                                                                                                                                                                                                    >  2:49445/1‑251 (MQ=255)
                  gATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAAGTCTGAAATATGAGTTAAAACAGGCGGGTATCAc                                                                                                                                                                                                                              <  1:407337/251‑1 (MQ=255)
                       ggCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                                          <  1:203993/250‑1 (MQ=255)
                       ggCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                                          >  1:504964/1‑250 (MQ=255)
                         cTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGCAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGaac                                                                                                                                                                                                                       >  2:179694/1‑249 (MQ=255)
                           tttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGACGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCACTTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCgg                                                                                                                                                                                                                     >  1:324063/1‑251 (MQ=255)
                             tataCGCGTGGCATCTTATTACCCCGTTCCGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGa                                                                                                                                                                                                                    <  1:212165/250‑1 (MQ=255)
                                   cgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                                                                                                                             <  1:346801/251‑1 (MQ=255)
                                              attaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACt                                                                                                                                                                                                  <  1:253468/251‑1 (MQ=255)
                                                    ccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAAACTCAgggt                                                                                                                                                                                            >  1:375787/1‑249 (MQ=255)
                                                         aCGATGTTCAGCCAGGTACGCCCCCCCCTCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAATCGCTTTACGCTGACGGGATGCCTGCCACAACTTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAGCAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACg                                                                                                                                                                                        <  1:478471/250‑1 (MQ=255)
                                                           gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAAAACTGCTGCGCCAGACACAGGTTAACAAAc                                                                                                                                                                                     >  2:353627/1‑251 (MQ=255)
                                                                        gTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCCGACTCAGGTTAACGAACCTGGAACGGTAGt                                                                                                                                                                        >  1:331804/1‑251 (MQ=255)
                                                                              aCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGcc                                                                                                                                                                    >  1:36477/1‑249 (MQ=255)
                                                                                   cccccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTaa                                                                                                                                                               <  1:248402/249‑1 (MQ=255)
                                                                                    ccccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACa                                                                                                                                                            <  1:188473/251‑1 (MQ=255)
                                                                                                gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAAAGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGcc                                                                                                                                                >  1:131536/1‑250 (MQ=255)
                                                                                                 ccGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGACATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACAATGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCACCGCTg                                                                                                                                               >  1:310972/1‑251 (MQ=255)
                                                                                                   gAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCa                                                                                                                                             <  1:263966/251‑1 (MQ=255)
                                                                                                         aaaTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGATTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCt                                                                                                                                       <  2:173027/250‑1 (MQ=255)
                                                                                                                gCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTGGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTaa                                                                                                                                 <  2:89769/250‑1 (MQ=255)
                                                                                                                 cGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTa                                                                                                                                  >  1:177023/1‑248 (MQ=255)
                                                                                                                  gAACTGTATGTGGTGCCAGGCGACCTGATCCGTTTTACGCTGACTGGATGCCTGTCACAACGGTCTGAGCCGCTCCCGTTGGCTGTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATAGGAGTTAAAACAGGCGGGTATCACCTGGAGAGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAt                                                                                                                              <  2:314848/251‑1 (MQ=255)
                                                                                                                  gAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAAAACTGCAGCGCCAGAATCAGGTTAACGAAACTGGGACGGGAATTGCCAGTAAAAAGTCGGCCCCGCTGCCCGATCTgtgta                                                                                                                                  >  1:493276/1‑243 (MQ=255)
                                                                                                                  gAACGGGATGTTGTGCCTGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTGTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACATGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAt                                                                                                                              <  1:337333/251‑1 (MQ=255)
                                                                                                                      gggATGTTGTCCCCGGCGACCTGACCCGCTTTACGCTGACGTGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATg                                                                                                                          <  2:430752/250‑1 (MQ=255)
                                                                                                                       tgATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGATGCGCCAGACTCAGGTTACCAAAACTGGACAGGAAGTTGCCCGTAACAAGTCAGCCCCGCCGCCCGAACTGCTTATGAtttagc                                                                                                                         >  2:126143/2‑247 (MQ=255)
                                                                                                                                     gggACCGGAACCGCGTTGCGCGGGCGGGATGCCTGCCAGAGCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCTAGACTCAGGTTAACGAACCTGGAACTGTAGTTGCCAGGAAACATTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGac                                                                                                           <  2:324063/249‑1 (MQ=255)
                                                                                                                                            gAACCGCTTTACGCTGACGGGATGCCGGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGTTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGTACAACATGa                                                                                                    <  1:379303/251‑1 (MQ=255)
                                                                                                                                               ccGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCCGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAATCGGGCCCGCCGGCCGATCTGCTTTAGATTTTGCTGAAAAAGTTGGACAACATTATcc                                                                                                 >  2:182902/1‑250 (MQ=255)
                                                                                                                                                     tccGCTGACGGGATGCCTGCCACAGCGTTTTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCTAGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTGAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTTAAGATTATGCGGAAAAAGTCGGACAACAGGATCGCCGAt                                                                                            <  2:230507/248‑1 (MQ=255)
                                                                                                                                                      aCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATAc                                                                                          >  1:384838/1‑251 (MQ=255)
                                                                                                                                                                           caACGTTCTGAGCCGCTCCCGTTGGCTTGTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCCATTCTGAAATATGAGATAAAACCGGCGGGTAGCCCCTGGGGCGGAACACTGATGTGACCGAATCAGGTTCACGAACCTGGAACGGGAGGTGCCAGTAAACACTCGGCCACGCTGAACGAACTGGCTAAGGTTCTGTTGAAAAAGTCGGACAACATGATCGCCTACACGGTTTTCTGCATAATAgggc                                                                      >  2:118951/1‑248 (MQ=255)
                                                                                                                                                                                 tCTGAGCCGCTCCCGTTCGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGACCACCTGGAGCGGGACCATGATGCGCCAGACTCCGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGGCAACATGATCGCCGATACGGGTTTCCGCATGATAGGCCATgtgcg                                                               >  1:174885/1‑251 (MQ=255)
                                                                                                                                                                                       ccGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAAAACTGCTGCGCCAGACTCAGGTTAACGAACCTGGGAAGGTAGTTGCCAGTAAACAGGAGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAAACTGATCGCCGATACGGTTTTCCGCATGATAGGCCAAGCGCGCTTcc                                                          >  1:392968/1‑249 (MQ=255)
                                                                                                                                                                                            cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCCGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGGTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTGTTCCGCATGATAGGCCATGCGCGCTTCAATGTGc                                                   <  2:350051/251‑1 (MQ=255)
                                                                                                                                                                                               gTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGAATCAGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTg                                                >  1:60468/1‑251 (MQ=255)
                                                                                                                                                                                                            ggtcAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGc                                   <  2:4220/248‑1 (MQ=255)
                                                                                                                                                                                                                   aTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGTTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCg                            <  1:110341/251‑1 (MQ=255)
                                                                                                                                                                                                                                 accGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGGCAGACTCACGTTACCGAACCTGGAAAGGTAAGTGCCAGTAAACAGTCGGCCCCGCTGCACGAACTGCTTTAGAGTATGCTGAAAAAGTCGGACAACATGATCGCCGATACCGGTTTCCGCATTCTAGGCCCTGCGCGCTTCTATTGTCCTGGAACATGGCGGGCCggggccggcgcccagcgtcc              >  1:224786/2‑234 (MQ=255)
                                                                                                                                                                                                                                        cAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACCGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAAAATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCCGGTCGGACGCCGCGCGGCAGATCCTg       >  1:217494/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCGGGAACGTTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGTTAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCa  <  2:174066/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: