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.1% 193.5 / ‑6.2 53D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (28/23);  minor base C (0/1);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCG  >  NC_000913/3329500‑3329956
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
aGTGCCCCAAAGCCTGGTGATATGGCTTTTCTACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                   >  2:275856/1‑251 (MQ=255)
    ccccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCCGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGACATATGAGTTaaa                                                                                                                                                                                                                >  2:370745/1‑250 (MQ=255)
            ccTGGTTATATGGCCTTTATACGAGTGGCATCTTATCACCCCCTTACGATGTTCAGCCAAGTAC‑GCACTCTCCCCCGTGATTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGGGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGACGGTGCACTACTGAAATCTGGGTTAAAACAGGCggg                                                                                                                                                                                                       >  2:129364/1‑251 (MQ=255)
               ggggATATGGCTTTTATACCCGGTGCCTCTTTTTTCCCCGTTACGATTTTCAGCCCGGTCCCGC‑CCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGATCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGTATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGt                                                                                                                                                                                                      <  2:274079/246‑1 (MQ=255)
                   atatGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAcc                                                                                                                                                                                                >  2:174288/1‑251 (MQ=255)
                    tatcgCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAccc                                                                                                                                                                                               >  2:461464/5‑250 (MQ=255)
                     atGGCTTTTATACGCGTGGCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTg                                                                                                                                                                                              <  1:275856/251‑1 (MQ=255)
                        gCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAg                                                                                                                                                                                           <  1:96470/251‑1 (MQ=255)
                          ttttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGAGCg                                                                                                                                                                                         >  1:375350/1‑251 (MQ=255)
                           tttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCg                                                                                                                                                                                         <  1:255501/250‑1 (MQ=255)
                             tataCGTGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTACCAACTGGAGCggg                                                                                                                                                                                       >  1:134014/1‑249 (MQ=255)
                                aCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAaca                                                                                                                                                                                    <  2:285040/250‑1 (MQ=255)
                                    gTGGCATCCTATTAACACGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGAGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCCGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGCTCCCCTGGAGCGGGACAccgct                                                                                                                                                                               >  2:174234/1‑251 (MQ=255)
                                          tCTTATTACCCCGTTACGTTGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCAGACGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGCTCTGAGCCGCTCCCGCTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGCGTTAAAACAGGCGGGTGTCACCCGGAGCGGCACACTGTTGCGcc                                                                                                                                                                          >  2:406809/1‑250 (MQ=255)
                                             tatttCCCCTTTGCGATGTTCAGCAAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCTGCCACAACGTTCGGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTTCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGTAACACTGCTGCGCCa                                                                                                                                                                         <  2:402099/248‑1 (MQ=255)
                                              attaCCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGGGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑Ct                                                                                                                                                                    <  2:376581/251‑1 (MQ=255)
                                                           gATGTTAAGCCCGATAC‑GCACCATCCCCCGTGGTTCTGCCGCAGCGCAATACTGCGAACTGGATGTGGTGCAAGGCGACCTGCCCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGA‑CTC‑AGGTTAACGAAc                                                                                                                                                      >  2:191375/1‑251 (MQ=255)
                                                                  aGCCAGGTAC‑GCACCCTCCCCCGTGGGTCTGCCCAGGCGCAATACGGCGAACTGGATGGGTTGCCAGGCGACCTGATCCGCTTTACGCGGACGGGTTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATTAGTTAAAACAGGCGGGTATCACCGGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTATCGATCCTTGAAc                                                                                                                                               <  2:180658/251‑1 (MQ=255)
                                                                    gccagttacggtcccTCCCCCGGGGTGCTGCCGAAGCGCAATACTGCGAACTGGATGTTGTGCCAGGCGCCCTGAACCGCTTTACGCTGGCGGGATGCCTGCCACAACGTGCTGGTCCGCTCCCGTTGGCTTTTGCCGTTCAGGAGGGAGCCAGCTATGCCGGTGCAATTCTGAAATAGGAGTTAAAAGAGGCGGGTATCTCCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGATCCTGGAACg                                                                                                                                              <  2:451888/239‑1 (MQ=255)
                                                                              cACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGAACACCGCTGCGCCAGA‑CTC‑CGGTTAACGGACCCGGAACGGGAGTTGccc                                                                                                                                    >  2:76149/1‑250 (MQ=255)
                                                                                   ccccccGTGGTTCAGCCGAAGCGCAATACTGCGAACTGGATGGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTaaa                                                                                                                               <  2:253387/250‑1 (MQ=255)
                                                                                             ttCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGA‑ATC‑AGGTTAACGAACCTGGAAAGGTAGTTGCCAGGAAACAGGCGGccc                                                                                                                     >  2:131875/1‑251 (MQ=255)
                                                                                                     aaGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACGGCTGCGCGAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCAc                                                                                                             <  1:191375/251‑1 (MQ=255)
                                                                                                      aGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACg                                                                                                            <  2:235135/251‑1 (MQ=255)
                                                                                                      aGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGGTAACGAACCTGGAAAGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACg                                                                                                            >  1:415152/1‑251 (MQ=255)
                                                                                                        cgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                          <  2:248588/251‑1 (MQ=255)
                                                                                                         gcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATc                                                                                                         <  1:461464/251‑1 (MQ=255)
                                                                                                         gcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAAAACTGCTGCGCCAGAACCA‑AGGTTA‑CGAAACCGGAACGGGAGTTGCCAGTAAAAAGGCGGCCCAGCTGCACGATc                                                                                                         >  2:260098/1‑251 (MQ=255)
                                                                                                               cTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CGC‑AGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACCGTCGGCCCCGCTGCACGATCTGCtt                                                                                                    >  1:90388/1‑250 (MQ=255)
                                                                                                               cTGCGAACTGGATGTGGAGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCCGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTTGAGGGGGACCCTGCGCCGCCCGA‑CCC‑AGGTTAACGAACCTGCAAAGGGAGGTGGCAGGAAACAAGCGGCACCCCTGCACGAtagtcttt                                                                                                   >  2:67403/1‑244 (MQ=255)
                                                                                                                                tgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTTAAGATTATGCTG‑aaaaa                                                                                  >  1:424259/1‑250 (MQ=255)
                                                                                                                                ggCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAACATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAg                                                                                 <  2:427479/250‑1 (MQ=255)
                                                                                                                                           ctgaaaagcattACGCTGACGGGATGCCTGCCACAACGCTCCGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGAGGGCATCACCTGGGGCGGGACCCTGCTGCGCCAGA‑ATC‑AGGGTCACGCAACTGGAAAGGGCGTTGCCCGTAAACAGTCGGCCCCCCTGCACGCTCTTCTTAGGGTTATTCTG‑AAAAAGACGGACCACAt                                                                      >  2:61001/11‑251 (MQ=255)
                                                                                                                                               aCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATc                                                                  <  1:370745/251‑1 (MQ=255)
                                                                                                                                                     ttACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCg                                                              <  2:330778/249‑1 (MQ=255)
                                                                                                                                                        cGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGGACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGGACGATCTGGTTAAGATTATGCTGGAAAAAG‑CGGACAACATTATCGCCGATACg                                                         >  2:128451/1‑251 (MQ=255)
                                                                                                                                                         gCTGACGGGATGCCTGCCACACCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTGTCACCTGGAGCGGCACCCTGCTGCGCCAGA‑CTC‑AGGTTAAAGAACCTGGAACGGTAGTTGCCCGTAAAAAATCCGCACCGCTGCACGATCTGCTTAAGATTATGCTG‑GAAAAGGCGGACAAAATGAACGCCGaatagg                                                        >  2:444517/1‑246 (MQ=255)
                                                                                                                                                         gCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACTATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACgg                                                        <  1:316310/251‑1 (MQ=255)
                                                                                                                                                                ggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGA‑ATC‑AGGTTAACGAACCTGGAAAGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGCTACGGTTTTCCg                                                 >  1:399485/1‑251 (MQ=255)
                                                                                                                                                                   tgcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGA‑CTC‑AGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTTAGATTATGCTG‑AAAAAGTCGGACAACATTATCGCCGCTCCGGGTTTCCGCAt                                              >  1:303085/1‑251 (MQ=255)
                                                                                                                                                                       tgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCCGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGGCGGACAACATGATCGCCGAAACGGGTTTCCGCATGAt                                           >  1:343047/1‑250 (MQ=255)
                                                                                                                                                                       ggccACAACGTTCTGCGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATa                                          <  2:332744/250‑1 (MQ=255)
                                                                                                                                                                              aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCCGCGCCCGA‑CTCCAGGT‑AACGAACCTGGGACGGGAGTTGCCAGGAAACCGTCGGACCCGCTGCACGCTCTGCTTAAGGTTATGCTG‑ACAAAATAGGACAACATGAGCGCCGATACGGTTTTCCGCATTATAGGCCAtt                                   >  2:450241/1‑250 (MQ=255)
                                                                                                                                                                                     gAGCCGCTCCAGCTGGCTCTTGCCGTGCAGGATGGAGCCAGCTATGACGGCGCAATTCTGAAATATGAGTTCAAACAGGCGGGCATTTCCTGGAGAGGGACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGCAGTTGCCCGTAAACAGTCGGCCCCGCTGCCCGCTCCGCTTAAGATTATGCTG‑AAAAAGACGGACAACATGCTCGCCGCTACGGTTTTCCGCATGATCGGCCATGCGCCCtt                            >  2:41369/1‑251 (MQ=255)
                                                                                                                                                                                      aGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTc                           <  1:262052/251‑1 (MQ=255)
                                                                                                                                                                                         cGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTTGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGCTACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAt                        >  2:346851/1‑251 (MQ=255)
                                                                                                                                                                                            tCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGGGGAACACTGCTGCGCCCGA‑ATC‑AGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCCCTGCCCGATCTGCTTAAGATTTTGCTG‑AAAAAGTCGGACAACATTATCGCCGAAACGGTTTTTCGCCTTGTAAGCCCAGCGCGCTTCAAtgtg                     >  2:327494/1‑251 (MQ=255)
                                                                                                                                                                                             cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                     <  2:365097/250‑1 (MQ=255)
                                                                                                                                                                                               cGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCt                  <  1:346851/251‑1 (MQ=255)
                                                                                                                                                                                                         tGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGGACATgg         <  1:131875/250‑1 (MQ=255)
                                                                                                                                                                                                               gCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATTCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTGCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCg  <  1:444517/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGGCAACATGATCGCCGATACGGTTTTCCGCATGATa                                          >  2:86806/1‑200 (MQ=255)
                                                                                                                                                                                                                          ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATa                                          <  1:86806/200‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTAC‑GCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGA‑CTC‑AGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTG‑AAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCG  >  NC_000913/3329500‑3329956

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