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.9% 145.5 / ‑6.1 48D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (16/31);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGC  >  NC_000913/3329498‑3329971
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
aCAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACTTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAg                                                                                                                                                                                                                                    <  1:180888/251‑1 (MQ=255)
aCAGTGCCCCAAAGCCGGGGGATATGGGTTTTATCCGCGTGGCCTCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGAt                                                                                                                                                                                                                                       <  2:158243/248‑1 (MQ=255)
 cagggCCCCAAAGCCTGGTGATATGGCTTTAATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGGGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGt                                                                                                                                                                                                                                   <  2:259342/247‑1 (MQ=255)
      cgtccctcttccgatctatGGCTTTTATACGCGTGGCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGACCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaaa                                                                                                                                                                                                                              <  2:410700/235‑1 (MQ=255)
      ccccAAAGCCTGGTGCTATGGCTTTTATACGCGTGGCATCTTCTTACCCCGTTACGATGTTCAGCCCGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCCATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCCGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCACTTCTGAAATATGCGTTaaa                                                                                                                                                                                                                               >  2:215702/1‑250 (MQ=255)
          aaaGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAgg                                                                                                                                                                                                                          <  1:470360/251‑1 (MQ=255)
             gCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTAAGCTAGGTACGCACCCTCCACCGTGGTTCTGCCGAAGAGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCAAAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGAGGGATGGAGCCAGCTATGCCGGTGAAATTCTGAAATTTGTGTTAAAACAggagg                                                                                                                                                                                                                       >  2:227078/1‑248 (MQ=255)
               cTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGt                                                                                                                                                                                                                     >  1:4267/1‑251 (MQ=255)
                 gttGATATGGCTTTTATACGCGTGGCACCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTAt                                                                                                                                                                                                                   <  2:286894/249‑1 (MQ=255)
                 ggtggtagggCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTAt                                                                                                                                                                                                                   <  2:224978/243‑1 (MQ=255)
                   ggATATGGATTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGTCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCa                                                                                                                                                                                                                 <  1:320339/250‑1 (MQ=255)
                       atGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaaa                                                                                                                                                                                                                              >  1:410700/1‑234 (MQ=255)
                           cTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATCACCTggggg                                                                                                                                                                                                         >  1:466003/1‑248 (MQ=255)
                              gtaaccgcccggcagtttaggtcaccGTTACGATGTGAACCAAGGTACGCACCCACCCCCCAGGTTCTGCTGATGCCTAATTCTGCGACCGGGTTGTGGTGCAAGTCGACCTGATCCGCTGTCCGCTGACGGGAGAACTGCCACAACGTTCTGAGCCGATCCCGTTGGCTTTTGCCGTGCAGGATGGAGCTAGCTATGCCTGTGCAATTCCGAAATATGAGTTAAAACAGGTCGGTATCACCTGGAGCgg                                                                                                                                                                                                       <  2:266703/226‑1 (MQ=255)
                               tataCGCGTGGCATCTTATTACCCCGATACGATGTTCAGCCAGGTACGCACCCTCCCACGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGACGTGACGCCAGGCGACATAAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCAGAGCCGCACCCCTTGGCTTTTGCCGTGCAGGTAGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAccagggtcgga                                                                                                                                                                                                      >  1:457010/1‑241 (MQ=255)
                                 aaCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAaca                                                                                                                                                                                                   <  1:170079/250‑1 (MQ=255)
                                   cgcgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGTGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATTACCTGGGGCGGAgcacc                                                                                                                                                                                                 >  2:337668/1‑250 (MQ=255)
                                           aTCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTATAACAGGCGGGTATCTCCCGGAGAGGaaaaatgcgggcgc                                                                                                                                                                                         >  2:165641/1‑239 (MQ=255)
                                                            cGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTGTCACCTGTCGCGGAACCCTGCTGCGCCA‑GACTAAGGTTAAcat                                                                                                                                                                        >  1:293880/1‑248 (MQ=255)
                                                              aTGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGTGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAAcc                                                                                                                                                                     <  1:337668/251‑1 (MQ=255)
                                                                   cAGCCAGGCACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAGGCCTGCCACAACGTTCTGAGCCGCGCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGa‑a                                                                                                                                                               <  1:268990/251‑1 (MQ=255)
                                                                     gCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑Cg                                                                                                                                                            <  1:228310/251‑1 (MQ=255)
                                                                                             gTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGcc                                                                                                                                    <  1:150659/251‑1 (MQ=255)
                                                                                                 tGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGGACACTGCTGCGCCA‑GACTCCGGGTAACGAAACTGGA‑A‑CGGGTGGTGCCAGTTAACAGGCGGGCCCGc                                                                                                                                >  1:126330/1‑251 (MQ=255)
                                                                                                           caacaCTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGTTTTACGCTGGCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCt                                                                                                                      <  2:176611/247‑1 (MQ=255)
                                                                                                               aCTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                                                                                                   <  1:372706/250‑1 (MQ=255)
                                                                                                                    gaagtGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCATGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAt                                                                                                             <  1:191695/247‑1 (MQ=255)
                                                                                                                         ggATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGCGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATACCCTGGAGCGGAACACTGCTGCGCCC‑GAATCCGGGTAACGAACCTGGA‑A‑CAGGAGATGCCAGTAAACAGACGGCCCCGCTGCAGCAACTGCTTAAggtttagc                                                                                                        >  2:269098/1‑244 (MQ=255)
                                                                                                                             gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaa                                                                                                    <  1:429609/251‑1 (MQ=255)
                                                                                                                                gtgCCCGGAGACCTGAAACGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGCTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCCCCCGGGGCGGGACACTGGTGAGCCAAGA‑TAAGGTTAAAGCACCTGGA‑AAGGGGAGTT‑CCAGTGAACAGTCGGCCCAGCCGCACGaaccccctaagaatacacttcaaa                                                                                                  >  2:342331/1‑227 (MQ=255)
                                                                                                                                     aGGCGACCTGACCCGCTTTACGGTGAGGGAAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGGTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTTGGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCgg                                                                                            <  2:271943/251‑1 (MQ=255)
                                                                                                                                      ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGa                                                                                           <  1:86368/251‑1 (MQ=255)
                                                                                                                                              gaaccgtgtatctcTGACGGGATGCCTGTCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGa                                                                                   <  2:466003/238‑1 (MQ=255)
                                                                                                                                                 ccGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTGTTTCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTGAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCg                                                                                <  1:56598/251‑1 (MQ=255)
                                                                                                                                                       tACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTTGAGCGGGACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑AGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGGTCTGCTTTAGATTTTGCTGAAAAAGTCGGACAAAATTATAGCCGAtc                                                                          >  1:260735/1‑250 (MQ=255)
                                                                                                                                                         cGCTGCCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAACGGAACACTGCTGCGCCA‑GACTCCGGGTAAAGAACCTGGA‑A‑CGGGAGTTGCCAGTAAACAGTCGGCCCCCCTGCACGCTATGCTTAAGATTATGGTGAAAAAACCGGAAAACATGATCGCCGATACg                                                                        >  1:211865/1‑251 (MQ=255)
                                                                                                                                                          gCTGACGGGATGCCTGCCACAACTTTCTGAGCCGCTCCCGTTGGCTTTTTCCGTGCAGGATGGAGCCAGCTATGCCGTTGCAATTCTTAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCTCTGCACGATCTGCTTAAGATTATGCTGAAAAATTCGGACAACATGATCGCCGATACgg                                                                       <  1:315508/251‑1 (MQ=255)
                                                                                                                                                                  gATGCCTGCCACAACGTTCTGAGCCGCTCCCGGTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGc                                                               <  1:122336/251‑1 (MQ=255)
                                                                                                                                                                       ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGaca                                                                                         <  1:358298/220‑1 (MQ=255)
                                                                                                                                                                         gccACAACGTTCTGAGCCCCTCCCGTTGGCTTTTGCCTTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAg                                                        <  1:95628/251‑1 (MQ=255)
                                                                                                                                                                          ccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGAAATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTCTCACCCGGGGCGGCACCATGCTGCGCCA‑GACTCAGGTTAACGAACCTGGACA‑CG‑TAGTTGCCAGCAAAAAGTCTGCCCCCCTGCACGATCTGCTTAAGATTATACTGAAAAAGTCGGAAAACATTATCGCCGCGACGGTTTTccccaataaa                                                         >  1:63004/1‑241 (MQ=255)
                                                                                                                                                                                             ccccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                                    <  1:107000/250‑1 (MQ=255)
                                                                                                                                                                                              cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGc                                   <  1:408257/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGTGCAGGATTGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAATCAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGTTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCgg                     <  2:472554/251‑1 (MQ=255)
                                                                                                                                                                                                                  aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaa                                                                                        >  1:263741/1‑178 (MQ=255)
                                                                                                                                                                                                                  aGGATGGAGCCAGCTATGCCGGTGCAAGTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaa                                                                                        <  2:263741/178‑1 (MQ=255)
                                                                                                                                                                                                                   ggATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGt              <  1:431545/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGTAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGc  <  1:342331/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
ACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCA‑GACTCAGGTTAACGAACCTGGA‑A‑CGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGC  >  NC_000913/3329498‑3329971

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: ‑