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

ACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGC  >  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)
                                                            cGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTGTCACCTGTCGCGGAACCCTGCTGCGCCAGACTAAGGTTAAcat                                                                                                                                                                        >  1:293880/1‑248 (MQ=255)
                                                              aTGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGTGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAcc                                                                                                                                                                     <  1:337668/251‑1 (MQ=255)
                                                                   cAGCCAGGCACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAGGCCTGCCACAACGTTCTGAGCCGCGCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGaa                                                                                                                                                               <  1:268990/251‑1 (MQ=255)
                                                                     gCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACg                                                                                                                                                            <  1:228310/251‑1 (MQ=255)
                                                                                             gTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGcc                                                                                                                                    <  1:150659/251‑1 (MQ=255)
                                                                                                 tGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGGACACTGCTGCGCCAGACTCCGGGTAACGAAACTGGAACGGGTGGTGCCAGTTAACAGGCGGGCCCGc                                                                                                                                >  1:126330/1‑251 (MQ=255)
                                                                                                           caacaCTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGTTTTACGCTGGCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCt                                                                                                                      <  2:176611/247‑1 (MQ=255)
                                                                                                               aCTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                                                                                                   <  1:372706/250‑1 (MQ=255)
                                                                                                                    gaagtGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCATGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAt                                                                                                             <  1:191695/247‑1 (MQ=255)
                                                                                                                         ggATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGCGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATACCCTGGAGCGGAACACTGCTGCGCCCGAATCCGGGTAACGAACCTGGAACAGGAGATGCCAGTAAACAGACGGCCCCGCTGCAGCAACTGCTTAAggtttagc                                                                                                        >  2:269098/1‑244 (MQ=255)
                                                                                                                             gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaa                                                                                                    <  1:429609/251‑1 (MQ=255)
                                                                                                                                gtgCCCGGAGACCTGAAACGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGCTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCCCCCGGGGCGGGACACTGGTGAGCCAAGATAAGGTTAAAGCACCTGGAAAGGGGAGTTCCAGTGAACAGTCGGCCCAGCCGCACGaaccccctaagaatacacttcaaa                                                                                                  >  2:342331/1‑227 (MQ=255)
                                                                                                                                     aGGCGACCTGACCCGCTTTACGGTGAGGGAAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGGTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTTGGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCgg                                                                                            <  2:271943/251‑1 (MQ=255)
                                                                                                                                      ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGa                                                                                           <  1:86368/251‑1 (MQ=255)
                                                                                                                                              gaaccgtgtatctcTGACGGGATGCCTGTCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGa                                                                                   <  2:466003/238‑1 (MQ=255)
                                                                                                                                                 ccGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTGTTTCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTGAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCg                                                                                <  1:56598/251‑1 (MQ=255)
                                                                                                                                                       tACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTTGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAAGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGGTCTGCTTTAGATTTTGCTGAAAAAGTCGGACAAAATTATAGCCGAtc                                                                          >  1:260735/1‑250 (MQ=255)
                                                                                                                                                         cGCTGCCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAACGGAACACTGCTGCGCCAGACTCCGGGTAAAGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCCCTGCACGCTATGCTTAAGATTATGGTGAAAAAACCGGAAAACATGATCGCCGATACg                                                                        >  1:211865/1‑251 (MQ=255)
                                                                                                                                                          gCTGACGGGATGCCTGCCACAACTTTCTGAGCCGCTCCCGTTGGCTTTTTCCGTGCAGGATGGAGCCAGCTATGCCGTTGCAATTCTTAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCTCTGCACGATCTGCTTAAGATTATGCTGAAAAATTCGGACAACATGATCGCCGATACgg                                                                       <  1:315508/251‑1 (MQ=255)
                                                                                                                                                                  gATGCCTGCCACAACGTTCTGAGCCGCTCCCGGTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGc                                                               <  1:122336/251‑1 (MQ=255)
                                                                                                                                                                       ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGaca                                                                                         <  1:358298/220‑1 (MQ=255)
                                                                                                                                                                         gccACAACGTTCTGAGCCCCTCCCGTTGGCTTTTGCCTTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAg                                                        <  1:95628/251‑1 (MQ=255)
                                                                                                                                                                          ccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGAAATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTCTCACCCGGGGCGGCACCATGCTGCGCCAGACTCAGGTTAACGAACCTGGACACGTAGTTGCCAGCAAAAAGTCTGCCCCCCTGCACGATCTGCTTAAGATTATACTGAAAAAGTCGGAAAACATTATCGCCGCGACGGTTTTccccaataaa                                                         >  1:63004/1‑241 (MQ=255)
                                                                                                                                                                                             ccccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                                    <  1:107000/250‑1 (MQ=255)
                                                                                                                                                                                              cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGc                                   <  1:408257/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGTGCAGGATTGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAATCAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGTTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCgg                     <  2:472554/251‑1 (MQ=255)
                                                                                                                                                                                                                  aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaa                                                                                        >  1:263741/1‑178 (MQ=255)
                                                                                                                                                                                                                  aGGATGGAGCCAGCTATGCCGGTGCAAGTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaa                                                                                        <  2:263741/178‑1 (MQ=255)
                                                                                                                                                                                                                   ggATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGt              <  1:431545/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGTAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGc  <  1:342331/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
ACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGC  >  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: