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% 217.5 / NA 67D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (26/39);  total (27/40)
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.88e-01

CAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAG  >  NC_000913/3329499‑3329988
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
cAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGt                                                                                                                                                                                                                                                 >  1:353462/1‑251 (MQ=255)
 aGTGCCCCAAATCCTGGTGATTTGGCTTTTATACGCGTGGCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCGCCGTGGGTCTGCCGAAGCGCAATACTGCGAACGGGATGTGTTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                <  2:359260/251‑1 (MQ=255)
 aGTGCCCCAAAGCCTGGTGATATGGCTTTTATTCGCGTGGCATCTTATTACCCCGTTACTATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                <  1:34257/251‑1 (MQ=255)
 aGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                >  2:351997/1‑251 (MQ=255)
  atgccccaaagaacggagATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTAACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGATCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                                                >  2:267733/18‑250 (MQ=255)
     ccccAAAGCCTGGTGATATGGCTTTTATACGCTTGGCACCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTTCCCCGTGTATCTGCCGAAGCGCAATTCTGCGAACTGGATGTGGTGCCAGGCGACTATAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCTCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGATTTaaaa                                                                                                                                                                                                                                            <  1:425935/251‑1 (MQ=255)
     ccccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaaa                                                                                                                                                                                                                                            <  1:445683/251‑1 (MQ=255)
       ccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAAc                                                                                                                                                                                                                                           >  1:368504/1‑250 (MQ=255)
             ccTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGCTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCggg                                                                                                                                                                                                                                    >  2:363034/1‑251 (MQ=255)
             ccTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCggg                                                                                                                                                                                                                                    <  1:187870/251‑1 (MQ=255)
                      atGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTg                                                                                                                                                                                                                           >  1:3848/1‑251 (MQ=255)
                           ttttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCg                                                                                                                                                                                                                      >  1:326249/1‑251 (MQ=255)
                               ataCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGCGTTAAAACAGGCGGGGATCAACTGGGGCGGAac                                                                                                                                                                                                                  >  1:312276/1‑251 (MQ=255)
                                taCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCgggac                                                                                                                                                                                                                  >  1:112692/1‑247 (MQ=255)
                                 aCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAacac                                                                                                                                                                                                                <  1:73704/251‑1 (MQ=255)
                                 aCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCGGAAATATGAGTTAAAACAGGCGTGTATCACCTGGAGCGGAacac                                                                                                                                                                                                                <  2:238652/251‑1 (MQ=255)
                                     gTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgct                                                                                                                                                                                                            <  1:10321/251‑1 (MQ=255)
                                             ttattaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCa                                                                                                                                                                                                                               <  2:243185/224‑1 (MQ=255)
                                             ttattaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCa                                                                                                                                                                                                                               >  1:243185/1‑224 (MQ=255)
                                                    cccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCCGGAGCGGAAAACTGCTGCGCCAGACTCCGGt                                                                                                                                                                                             >  1:316788/1‑251 (MQ=255)
                                                       gTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTa                                                                                                                                                                                           >  1:409275/1‑250 (MQ=255)
                                                          aCGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGa                                                                                                                                                                                       <  1:295113/251‑1 (MQ=255)
                                                          aCGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACGGGATGGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTAGCACCTGGAGCGGAACACTGCTGCGCCAGTCGCAGGTTAACGa                                                                                                                                                                                       <  2:312276/251‑1 (MQ=255)
                                                          aCGAGGTTCAGCCAGGTACGCACCCTCCCCCGTGGTGCTGCCGAAGCGCAATACTGCGAACGGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTGAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGa                                                                                                                                                                                       <  2:415737/251‑1 (MQ=255)
                                                             aTGTTCAGCCAGGTACGCTCCCGCCCCCGTGGTTCGGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCGGCCACAGCTTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTaa                                                                                                                                                                                                                                              <  1:303317/193‑1 (MQ=255)
                                                                         gTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGGCAGACTCAGGGTAACGAAACTGGAACGGGAGt                                                                                                                                                                        >  1:149105/1‑251 (MQ=255)
                                                                            cGCACCCTCCCCCGTGGCTCTGCCGAAGCGCAATCCTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGACGTGCAGGATGGAGCCAGCTATGCCGTTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGCACACTGCTGCGCCCGAATCAGGTTaaaaaaccttggcagggagatgg                                                                                                                                                                     >  2:315700/1‑230 (MQ=255)
                                                                                cccTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGGGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGt                                                                                                                                                                 <  2:299487/251‑1 (MQ=255)
                                                                                  cTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTa                                                                                                                                                                >  1:340253/1‑250 (MQ=255)
                                                                                       ccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCCGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACACTCCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAAc                                                                                                                                                             >  1:432010/1‑248 (MQ=255)
                                                                                            gTTCTGCCGAAGCGCAATACGGCGAACTGGATGTGGTGCCAGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTTCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACGGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGTTAGTGGCCAGTAAACATTCGGcc                                                                                                                                                     <  2:293089/251‑1 (MQ=255)
                                                                                              tCTGCCGAAGCGCAATTCTGCGAACTGGATGTGGGGCCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCggtcc                                                                                                                                                    <  2:249101/250‑4 (MQ=255)
                                                                                              gcTGCCGAAGCGCAATACTGCGAACTGGGTGTGGGGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGccc                                                                                                                                                    <  2:353462/249‑1 (MQ=255)
                                                                                                 gCCGAAGCGCAATACTGCGAACTGGATGTTGGGCCAGGCGACCTGAGCCGCTTTACGCGGACGGGATGCCTGCCCCAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACACGGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGTCCCCGCt                                                                                                                                                <  2:453172/251‑1 (MQ=255)
                                                                                                       gcgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGa                                                                                                                                          <  1:2025/251‑1 (MQ=255)
                                                                                                          cAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCt                                                                                                                                       <  1:180471/251‑1 (MQ=255)
                                                                                                              aCTGCGAACTGGATGTGGTGCCAGGCGACCTTAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCGGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                                                                                                                    <  2:326249/250‑1 (MQ=255)
                                                                                                                 gCGAACTGGTGGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTa                                                                                                                                  <  2:3848/249‑1 (MQ=255)
                                                                                                                                   cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCg                                                                                                              <  1:285155/251‑1 (MQ=255)
                                                                                                                                             gAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACAt                                                                                                      <  1:95713/249‑1 (MQ=255)
                                                                                                                                             gAACCGCTTTACGCTGACGGGATGCCGGCCACAACGTTCTGAGCCGCGCCCGTTGGCTTGTGCCGTGCAGGAGGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTGAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGt                                                                                                                                                                           <  2:313789/180‑1 (MQ=255)
                                                                                                                                                       aCGCTGCCGGGATGCCTGGCCCAACGTTCTGAGCCGCTCCCGCTGGCTTTTGCCGTGCCGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAACATGAGTTAAAACAGGCGGGGATTCCCTGGGGCGGGACACTGCTGCGCCAGACTCAGGGTAACGAACCTGTAACGGTAGATGCCAGTACAAAATCAGCCCCGCTGCCCGATCTCCTTTAGAATATGCTGCAAAACGCCGACAACATGTTCGcctctac                                                                                          >  2:213216/1‑246 (MQ=255)
                                                                                                                                                          cTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTTGAGCGGAACACTGCTGCGCCAGACTAAGGTTAAAGAACCTGGGACGGTAGTTGCCAGTAAACAGACGGCCCCGCTGCACGATCTGCTTAAAGTTATGCTGAAAAAGTCGGACAAACTGATCGCCGGTACGGt                                                                                       >  2:124538/1‑251 (MQ=255)
                                                                                                                                                           tGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTTAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGtt                                                                                      >  1:318980/1‑251 (MQ=255)
                                                                                                                                                            tacgtgATGCATGCCACAACGTGCTGAGCCGCTCCCGTTGGCTTTTTCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACATGCGTGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACATGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGtt                                                                                      <  2:346462/245‑1 (MQ=255)
                                                                                                                                                            gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGttt                                                                                     <  2:112692/251‑1 (MQ=255)
                                                                                                                                                                        gccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGAt                                                                           <  1:244563/249‑1 (MQ=255)
                                                                                                                                                                        gccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATa                                                                          <  2:340253/250‑1 (MQ=255)
                                                                                                                                                                        gccACAAAGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGTCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAg                                                                         <  2:152855/251‑1 (MQ=255)
                                                                                                                                                                           acaACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGcc                                                                      <  2:292822/251‑1 (MQ=255)
                                                                                                                                                                                      aGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGGGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGGACGGGAGTTGCCAGTCAACAGTAGGCCCAGCTGCACGATCTGCTTAAGATTATGATGAAAAAGTCGGACAACCTGATCGCCGATACGGTTTTTCGCATGATAGGCCATGCGCGCTTc                                                           >  1:301383/1‑251 (MQ=255)
                                                                                                                                                                                       gCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGGTGCCAGCAAACAGTCGGCCCCCCCGCACGGACTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGAAACGGTTTTCCCCATGATAGGCCATACGCGCTTc                                                           >  2:253485/1‑250 (MQ=255)
                                                                                                                                                                                          acTCCCGTTGGCTTTTGACGTGCCGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAAAATCAGGGTAACGAAACTGGAACGGGAGGTGCCAGTAAACAGTCGGCCCCGCTGCACACACTGCTTAACGTTATGCTGAAAAAGTCGGACCACCTGATCGACGATACGGGTTTCAGCCAGATAGGCCCAGCGCGCTTCAAtg                                                       >  2:344567/2‑251 (MQ=255)
                                                                                                                                                                                                                  ggATGGGGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGt                               <  1:293272/251‑1 (MQ=255)
                                                                                                                                                                                                                   gATGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACTCTGCTGCGTCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCTCGCTGCTGGAGCTGCGTAAGAGTATGCTGAAAAAGTCGGTCAACAGGATCTCCGATACGGTTGTCCGCATGATAGGCCAGGCGCGCGTCAATGTGACTGGAACATGGCGGGCCGGGTc                              <  2:292372/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCg                             <  1:363034/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCTGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCg                             <  1:267733/251‑1 (MQ=255)
                                                                                                                                                                                                                     gtgAGCTAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAACCAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGGCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACATTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCAGGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCgg                            <  2:261533/249‑1 (MQ=255)
                                                                                                                                                                                                                       gAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGGTGCCAGTAAACAGACGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATACGCGCTTCAATGTGCCTGGAACAAGGCGGGGCGGGTCgggg                          >  2:402232/1‑249 (MQ=255)
                                                                                                                                                                                                                               tATGCCGGTGCAATTCTGCAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTACGCCAGAATCAGGGTAACGAACCTGGAAAAGTAGTTGCCAGGAAAACGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATAAACGCCGATACGGGTTTCCGCATGATAGGCCCTGCGCGGTTCAATATGCCCGGAACCTGGGGGGGCGGGTCGGAGGCCGTGCg                  >  2:450948/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCa               <  1:213216/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGtttt                                                                                    <  1:255837/179‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGtttt                                                                                    >  2:255837/1‑179 (MQ=255)
                                                                                                                                                                                                                                      gTGCAATTCTGAACTATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGACTCAGGTTAACGAACCTGGAACGGTAGTTTCCAGTAAACAGTCGGCCCCg                                                                                                                                                  >  1:90407/1‑116 (MQ=255)
                                                                                                                                                                                                                                      gTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCg                                                                                                                                                  <  2:90407/116‑1 (MQ=255)
                                                                                                                                                                                                                                          aaTTCTGAAATAGGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGGCTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTgc       <  2:9955/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGGACACTGCTGCGCCAGAATCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTTACCAGACGGCCCCGCTGCACGATTTGCTTAAGATTATGCCGCAAAAGTCGGACAACCTGATCGCCGATACGGTTTTCCGGATGCAAAGCCCTGGGCGGCTCAATGTGCCTGGAACATGGCGGGCCGGGCCGGCCGCCCGGCGTCCGATccctcgcccg  >  1:208010/1‑243 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
CAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAG  >  NC_000913/3329499‑3329988

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