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.5% 221.4 / NA 66D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (39/25);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGC  >  NC_000913/3329497‑3329989
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
tACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATg                                                                                                                                                                                                                                                         >  1:219231/1‑250 (MQ=255)
 aCAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATg                                                                                                                                                                                                                                                         >  1:498688/1‑249 (MQ=255)
 aCAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAg                                                                                                                                                                                                                                                       <  1:208526/251‑1 (MQ=255)
   aGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                     >  2:16623/1‑251 (MQ=255)
   aGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGAGGGAGCCAGCTATGCCGGTGCAATTCTGTAATATGAGtt                                                                                                                                                                                                                                                     >  2:164183/1‑251 (MQ=255)
      gCCCCAAAGCCTGGTTATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCGCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTaa                                                                                                                                                                                                                                                   <  2:456432/250‑1 (MQ=255)
        cccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAAc                                                                                                                                                                                                                                                <  1:30638/251‑1 (MQ=255)
                    tGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGTATGTGGTGCCCGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGTTATGCCGGTGCAATTCTGAAATATGAGTTAAACAAGGCGGGGATc                                                                                                                                                                                                                                     >  2:401278/1‑250 (MQ=255)
                     gATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGCTGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTATGAAATCTGAGTTAAAACAGGCGGGTATCAc                                                                                                                                                                                                                                   >  2:397448/1‑251 (MQ=255)
                                    cgcgTGGCATCTTATTACCCCGTTACGATGTTAAGACAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCCATTCTGACATATGAGTTCAAACAGGCGGTTATAAccatgaaggggaaaat                                                                                                                                                                                                                    >  1:361178/1‑237 (MQ=255)
                                                tattaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCa                                                                                                                                                                                                           <  1:286621/248‑1 (MQ=255)
                                                tattaCACCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCCGGCGACCTGAACCGCTTTACGCTGACGGGATGACTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGAACAATGCTGGGCCAGAc                                                                                                                                                                                                        >  1:180009/1‑251 (MQ=255)
                                                       ccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGtt                                                                                                                                                                                                 <  1:22922/251‑1 (MQ=255)
                                                           tACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACg                                                                                                                                                                                             >  1:132777/1‑251 (MQ=255)
                                                              gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAAc                                                                                                                                                                                          >  2:104405/1‑251 (MQ=255)
                                                                      gCCAGGTACGCACCCTCCCCCTTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACg                                                                                                                                                                                  <  2:219231/251‑1 (MQ=255)
                                                                            tACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTATGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACATGGAACGGGAGtt                                                                                                                                                                            >  1:484081/1‑251 (MQ=255)
                                                                            accgcgcgcTCCCCCGTGGTGCTGCCGAATCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTTAACCGCTGTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCGATGCCGGTGCAATTGTGAAATATGAGTTAAAACAGGCGGGTAGCACCGGGAGCGGAACACTGCTGCGCCAGACTGAGGTTAACGAACCTGGAACGGTAGtt                                                                                                                                                                            <  2:211213/243‑1 (MQ=255)
                                                                             aCGCACCCTCCCCCGTGGTGCTGCCGAAGCGCAATACTGCGAACGGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGtt                                                                                                                                                                            <  2:171559/250‑1 (MQ=255)
                                                                              cGCACTCTCCCCCGTGGTTCTGCCGAAGCGCACTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGACGGAGAAAATCTGAAATATGAGTTAAAAAAGGCGGGAATCACCTGGAGCGGAACAATGCTGCGACAGACTCCGGTTCACGCACCTGGAACggcaggtgg                                                                                                                                                                          >  2:41117/1‑244 (MQ=255)
                                                                                          cGTGGTTCTGCCGAAGCGCAATGCTGCGAACTGGATGTGGGGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTTCCACAACTTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTTAAATATGAGTTAACACAGGCTGGTAGCACCTTGAGCGGAACACTGCTGTTCCATACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTc                                                                                                                                                             <  1:466201/251‑1 (MQ=255)
                                                                                             ggTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAAGCGACCTGCACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGGGCGGAACACAGCCGCGCCACACTCAGGTTAACACAACTGGGACGGGAGTTGGCAG‑TAAACAGCCGGc                                                                                                                                                          >  1:269160/1‑251 (MQ=255)
                                                                                              gTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACACTTCTGCGCCAGAATCAGGTTAACGAACCTGGGACGGTAGTTGCCAG‑TAAACAGTCGGcc                                                                                                                                                         >  1:169778/1‑251 (MQ=255)
                                                                                                 cTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGAATCAGGGTAACGAACCTGGGACGGTAGTTGCCAG‑TAAACAGTCGGccccc                                                                                                                                                      >  1:416242/1‑250 (MQ=255)
                                                                                                   gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCt                                                                                                                                                    >  1:268212/1‑251 (MQ=255)
                                                                                                       aaGCGCCATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCAACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCCCCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAG‑GAAACAGTCGGCCCCCCTGCAc                                                                                                                                                >  2:72402/1‑251 (MQ=255)
                                                                                                            cAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACt                                                                                                                                                                                                       <  2:470542/192‑1 (MQ=255)
                                                                                                            cAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACACTGCTGCGCCAGACt                                                                                                                                                                                                       >  1:470542/1‑192 (MQ=255)
                                                                                                             aaTACTGCGAACAGGATGTGGTGCCAGGCGACTTGAACCTCTTTACGCTAACGGGATGCCTGCAACAACGGTCTGAGCCACTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAATACAGGCGGGTATTAACTTGATCGGAACACTTCTGAGACAGAATCAGGTTAACGATCCCGTAACGGGAACTGCCAG‑TAAACAATCGGCGCCGCTGCACgggctg                                                                                                                                          >  2:391178/1‑246 (MQ=255)
                                                                                                               tACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGCACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCt                                                                                                                                        >  2:247355/1‑251 (MQ=255)
                                                                                                                aCTGCGACCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATAACCAGGAGCGGAACACTGCTGCGCCAGACACAGGGTAACTAACACGGGACGGTCGTTGCCAGCGAAAC‑GTCGGCCCCGCTGAACGATcgcgtt                                                                                                                                       >  2:355667/1‑246 (MQ=255)
                                                                                                                    cGAACTGGCTGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGACACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGAGCCCGACTCAGGTTAACGAACCTGGAACGGCAATTGCCAG‑GAAACAAACGGCGCCGCCGCAccaacttctttac                                                                                                                                    >  2:249750/1‑238 (MQ=255)
                                                                                                                       aCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                                               <  2:193243/157‑1 (MQ=255)
                                                                                                                       aCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTgg                                                                                                                                                                                                                               >  1:193243/1‑157 (MQ=255)
                                                                                                                       aCTGGATGTGGGGCCAGGCGACCTGAACCGCTTGTCGCTGACGGGATGCCTGCCACAACGTGCTGAGCCGCGCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTa                                                                                                                                <  1:249750/251‑1 (MQ=255)
                                                                                                                            aTGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑Tg                                                                                                                          <  1:291581/251‑1 (MQ=255)
                                                                                                                            aTGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑Tg                                                                                                                          <  1:178854/251‑1 (MQ=255)
                                                                                                                                    ccAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTTACGAACCCGGGACGGGAGTTGGCAG‑GAAACAGTCGGCCCCGCCGCACGAGCTGTTTTAGATTCTTC‑TGAAAAAGTc                                                                                                                  >  2:84595/1‑251 (MQ=255)
                                                                                                                                      aGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAATTGCCAG‑TAAACAGTCGGCCCCGCTGCCCGATCTGCTTTAAATTATTC‑TGGAAAAGTCgg                                                                                                                >  1:437053/1‑251 (MQ=255)
                                                                                                                                         cGACCTGAACTGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATAACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGCTCTGCTTCAGATTATGC‑TTAAAAAGTAGGac                                                                                                              >  2:392450/1‑250 (MQ=255)
                                                                                                                                         cGACCTGAACCGCTTTCCGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACTGGCGGGTATCACCTGGAGCGGAACACTGCGGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGaca                                                                                                             <  1:349630/251‑1 (MQ=255)
                                                                                                                                            ccTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGCGCCCGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATCAACCGGAGCCGCCCTCTCATGCGCCAGAATCAAGTTAACGCACCTGGGACGGTGGTTGCCAA‑TAAAAAGGTGGCCCCGCGGCCAGCACTGCTTAAGATTAACG‑TGAAAAAGtacgaccac                                                                                                           >  1:205219/1‑242 (MQ=255)
                                                                                                                                              tGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACACTGCTGCGCCAGAATCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATg                                                                                                        >  1:30471/1‑251 (MQ=255)
                                                                                                                                                    gCTTTACGCTGCCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATTCCGGTGCAATTCGGAAATATGAGGTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGAGCTGCTTAAGATTATGC‑GGAAAAAGTCGGACAACAGGATCGcc                                                                                                  <  2:66499/251‑1 (MQ=255)
                                                                                                                                                         aCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATAc                                                                                             <  1:467933/251‑1 (MQ=255)
                                                                                                                                                                cGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTc                                                                                      <  1:450528/251‑1 (MQ=255)
                                                                                                                                                                cGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTc                                                                                      >  1:51258/1‑251 (MQ=255)
                                                                                                                                                                        ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGCTCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGAt                                                                              >  1:277748/1‑251 (MQ=255)
                                                                                                                                                                             acaACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCTCATGATAGGcc                                                                         <  2:197522/251‑1 (MQ=255)
                                                                                                                                                                               agcGTTCTGGGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGTCGGGTATCACCTGGGGCGTAACAGTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCAt                                                                       <  1:41117/249‑1 (MQ=255)
                                                                                                                                                                                     cTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACAATGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGCGGTGACAG‑GAAACAGTCGGCCCCGCCGCACCATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACACGTTTGCCGATACGGGTTTCCGCCTGAAAGGCCATgcgcgc                                                                 >  2:361190/1‑251 (MQ=255)
                                                                                                                                                                                         gCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAAAAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAAAATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGATTc                                                              >  1:370583/1‑250 (MQ=255)
                                                                                                                                                                                         gCCGCTCCCGTTGGCGTTTGCCGTGCAGGATGGAGCCATCTATGCCGGTTCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCTCCTGGAGCTGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCGGCTTAAGAGTATGG‑TGAAAAAGTGGGACATCATGATCGCCGATACGGTTTTCCGCATGAGAGGCCATACGCTCTTCa                                                             <  2:498688/251‑1 (MQ=255)
                                                                                                                                                                                               cccGTTGCCTTTTGCCTTGCAGGAGGGCGCCAGCTATGCCGGTGCATGTCTGAAATATGAGTTAAAGCAGGCGGGTATCCCCTGGAGCTGAACACTGCTGCGCCACATGCAGGTTAACGATCCTTGAACGGTAGTTGCCAG‑TAAACTGTCGGCCCCGCTCCACGTTCTTCTTAAGATTATTC‑TGAAAAAGTCGGACATCATGATCGCCTATATTGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGc                                                       <  2:184565/251‑1 (MQ=255)
                                                                                                                                                                                                  gTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAG‑GAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACCTGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTg                                                    >  1:463235/1‑251 (MQ=255)
                                                                                                                                                                                                           tGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACCCTGCTGCGCCAAAATCAGGGTAACGAACCTGGAAAGGTAGTTGCCCG‑TAAACAGTCGGCCCCGCTGCACGGTCTGCTTAAGATTATGCCTGAAAA‑GGCGGACAACATGAACCCCGATACGGGTTTTCGCATGATAGGCCATGCGCGCTTTAATGTGCCTGGAACAGGGc                                           >  2:372067/1‑251 (MQ=255)
                                                                                                                                                                                                              cGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGGCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGGAGGACAACATGATCGCCGGTACGGGTTTCCGCATGAAAGGCCATGCGCGCTTCAATATGCCTGGAACATGGC‑ggg                                       >  2:145231/1‑251 (MQ=255)
                                                                                                                                                                                                                     gATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTc                               >  1:198895/1‑251 (MQ=255)
                                                                                                                                                                                                                             cAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGGCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGACGCCg                       <  2:463235/251‑1 (MQ=255)
                                                                                                                                                                                                                              aGCTATGCCGGTGCAATTCTGAAAAAGGAGTTAAAACAGGCGGTTATCACCTTGCGCGGAACACTTCTGCTCCATACACAGGTTAACGAACCTGGAACGGGAGTTTCCAG‑TAAACATCCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATTATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGACGCCGt                      <  1:465774/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATTGC‑GGGC‑CGGGTCGGACGCCGTGCGTc                 <  2:51258/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGAGGAACACTGCTGCGCCAGACTCAGGTTTACGAACCTGGAACGGTAGTTGCCCG‑TAAACAGTCGGCCCCGCTGCACGATCTGATTAAGATTATGC‑TGAAAACATCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCCTCAATGGGCCTGGAACATGGC‑GGGCGCGGG‑CGGACGCCGTGCGGCAGaa             >  1:55036/1‑249 (MQ=255)
                                                                                                                                                                                                                                           cAATTCTGAAATATGAGTTAAAACAGGCGGGTCTCACCTGGAACGGAACACTGCTGCGCCAAACTCAGGTTAAAGAACCTGGAACGGGAGTTGCCCG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGA‑TGAAAAAATCGGACAACAAGATCGCCGATACGGGTTTTCGCATGATAAGCCATGCGCGCTTCAATGTGCCTGGAACATGGccgggc‑cgggacgacgccctgagtgcgaacctg          >  2:58407/1‑228 (MQ=255)
                                                                                                                                                                                                                                           cAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGACGCCGTGCGTCAGATCCTg         <  1:188673/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGc  <  1:22482/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAAAATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGGCGCCGAGCGTCAGATCCTGCGCCAGc  >  1:159827/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
TACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAG‑TAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGC‑TGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGC‑GGGC‑CGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGC  >  NC_000913/3329497‑3329989

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