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

AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCG  >  NC_000913/3329500‑3329956
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
aGTGCCCCAAAGCCTGGTGATATGGCTTTTCTACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGtt                                                                                                                                                                                                                   >  2:275856/1‑251 (MQ=255)
    ccccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCCGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGACATATGAGTTaaa                                                                                                                                                                                                                >  2:370745/1‑250 (MQ=255)
            ccTGGTTATATGGCCTTTATACGAGTGGCATCTTATCACCCCCTTACGATGTTCAGCCAAGTACGCACTCTCCCCCGTGATTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGGGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGACGGTGCACTACTGAAATCTGGGTTAAAACAGGCggg                                                                                                                                                                                                       >  2:129364/1‑251 (MQ=255)
               ggggATATGGCTTTTATACCCGGTGCCTCTTTTTTCCCCGTTACGATTTTCAGCCCGGTCCCGCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGATCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGTATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGt                                                                                                                                                                                                      <  2:274079/246‑1 (MQ=255)
                   atatGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAcc                                                                                                                                                                                                >  2:174288/1‑251 (MQ=255)
                    tatcgCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAccc                                                                                                                                                                                               >  2:461464/5‑250 (MQ=255)
                     atGGCTTTTATACGCGTGGCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTg                                                                                                                                                                                              <  1:275856/251‑1 (MQ=255)
                        gCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAg                                                                                                                                                                                           <  1:96470/251‑1 (MQ=255)
                          ttttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGAGCg                                                                                                                                                                                         >  1:375350/1‑251 (MQ=255)
                           tttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCg                                                                                                                                                                                         <  1:255501/250‑1 (MQ=255)
                             tataCGTGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTACCAACTGGAGCggg                                                                                                                                                                                       >  1:134014/1‑249 (MQ=255)
                                aCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAaca                                                                                                                                                                                    <  2:285040/250‑1 (MQ=255)
                                    gTGGCATCCTATTAACACGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGAGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCCGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGCTCCCCTGGAGCGGGACAccgct                                                                                                                                                                               >  2:174234/1‑251 (MQ=255)
                                          tCTTATTACCCCGTTACGTTGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCAGACGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGCTCTGAGCCGCTCCCGCTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGCGTTAAAACAGGCGGGTGTCACCCGGAGCGGCACACTGTTGCGcc                                                                                                                                                                          >  2:406809/1‑250 (MQ=255)
                                             tatttCCCCTTTGCGATGTTCAGCAAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCTGCCACAACGTTCGGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTTCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGTAACACTGCTGCGCCa                                                                                                                                                                         <  2:402099/248‑1 (MQ=255)
                                              attaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGGGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACt                                                                                                                                                                    <  2:376581/251‑1 (MQ=255)
                                                           gATGTTAAGCCCGATACGCACCATCCCCCGTGGTTCTGCCGCAGCGCAATACTGCGAACTGGATGTGGTGCAAGGCGACCTGCCCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGACTCAGGTTAACGAAc                                                                                                                                                      >  2:191375/1‑251 (MQ=255)
                                                                  aGCCAGGTACGCACCCTCCCCCGTGGGTCTGCCCAGGCGCAATACGGCGAACTGGATGGGTTGCCAGGCGACCTGATCCGCTTTACGCGGACGGGTTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATTAGTTAAAACAGGCGGGTATCACCGGGAGCGGAACACTGCTGCGCCAGACTCAGGTTATCGATCCTTGAAc                                                                                                                                               <  2:180658/251‑1 (MQ=255)
                                                                    gccagttacggtcccTCCCCCGGGGTGCTGCCGAAGCGCAATACTGCGAACTGGATGTTGTGCCAGGCGCCCTGAACCGCTTTACGCTGGCGGGATGCCTGCCACAACGTGCTGGTCCGCTCCCGTTGGCTTTTGCCGTTCAGGAGGGAGCCAGCTATGCCGGTGCAATTCTGAAATAGGAGTTAAAAGAGGCGGGTATCTCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGATCCTGGAACg                                                                                                                                              <  2:451888/239‑1 (MQ=255)
                                                                              cACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGAACACCGCTGCGCCAGACTCCGGTTAACGGACCCGGAACGGGAGTTGccc                                                                                                                                    >  2:76149/1‑250 (MQ=255)
                                                                                   ccccccGTGGTTCAGCCGAAGCGCAATACTGCGAACTGGATGGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTaaa                                                                                                                               <  2:253387/250‑1 (MQ=255)
                                                                                             ttCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGAATCAGGTTAACGAACCTGGAAAGGTAGTTGCCAGGAAACAGGCGGccc                                                                                                                     >  2:131875/1‑251 (MQ=255)
                                                                                                     aaGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACGGCTGCGCGAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCAc                                                                                                             <  1:191375/251‑1 (MQ=255)
                                                                                                      aGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACg                                                                                                            <  2:235135/251‑1 (MQ=255)
                                                                                                      aGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAAAGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACg                                                                                                            >  1:415152/1‑251 (MQ=255)
                                                                                                        cgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                          <  2:248588/251‑1 (MQ=255)
                                                                                                         gcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATc                                                                                                         <  1:461464/251‑1 (MQ=255)
                                                                                                         gcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAAAACTGCTGCGCCAGAACCAAGGTTACGAAACCGGAACGGGAGTTGCCAGTAAAAAGGCGGCCCAGCTGCACGATc                                                                                                         >  2:260098/1‑251 (MQ=255)
                                                                                                               cTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACGCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACCGTCGGCCCCGCTGCACGATCTGCtt                                                                                                    >  1:90388/1‑250 (MQ=255)
                                                                                                               cTGCGAACTGGATGTGGAGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCCGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTTGAGGGGGACCCTGCGCCGCCCGACCCAGGTTAACGAACCTGCAAAGGGAGGTGGCAGGAAACAAGCGGCACCCCTGCACGAtagtcttt                                                                                                   >  2:67403/1‑244 (MQ=255)
                                                                                                                                tgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTTAAGATTATGCTGaaaaa                                                                                  >  1:424259/1‑250 (MQ=255)
                                                                                                                                ggCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAACATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAg                                                                                 <  2:427479/250‑1 (MQ=255)
                                                                                                                                           ctgaaaagcattACGCTGACGGGATGCCTGCCACAACGCTCCGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGAGGGCATCACCTGGGGCGGGACCCTGCTGCGCCAGAATCAGGGTCACGCAACTGGAAAGGGCGTTGCCCGTAAACAGTCGGCCCCCCTGCACGCTCTTCTTAGGGTTATTCTGAAAAAGACGGACCACAt                                                                      >  2:61001/11‑251 (MQ=255)
                                                                                                                                               aCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATc                                                                  <  1:370745/251‑1 (MQ=255)
                                                                                                                                                     ttACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCg                                                              <  2:330778/249‑1 (MQ=255)
                                                                                                                                                        cGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGGACGATCTGGTTAAGATTATGCTGGAAAAAGCGGACAACATTATCGCCGATACg                                                         >  2:128451/1‑251 (MQ=255)
                                                                                                                                                         gCTGACGGGATGCCTGCCACACCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTGTCACCTGGAGCGGCACCCTGCTGCGCCAGACTCAGGTTAAAGAACCTGGAACGGTAGTTGCCCGTAAAAAATCCGCACCGCTGCACGATCTGCTTAAGATTATGCTGGAAAAGGCGGACAAAATGAACGCCGaatagg                                                        >  2:444517/1‑246 (MQ=255)
                                                                                                                                                         gCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACTATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACgg                                                        <  1:316310/251‑1 (MQ=255)
                                                                                                                                                                ggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGAATCAGGTTAACGAACCTGGAAAGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGCTACGGTTTTCCg                                                 >  1:399485/1‑251 (MQ=255)
                                                                                                                                                                   tgcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGACTCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTTAGATTATGCTGAAAAAGTCGGACAACATTATCGCCGCTCCGGGTTTCCGCAt                                              >  1:303085/1‑251 (MQ=255)
                                                                                                                                                                       tgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCCGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGATCGCCGAAACGGGTTTCCGCATGAt                                           >  1:343047/1‑250 (MQ=255)
                                                                                                                                                                       ggccACAACGTTCTGCGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATa                                          <  2:332744/250‑1 (MQ=255)
                                                                                                                                                                              aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCCGCGCCCGACTCCAGGTAACGAACCTGGGACGGGAGTTGCCAGGAAACCGTCGGACCCGCTGCACGCTCTGCTTAAGGTTATGCTGACAAAATAGGACAACATGAGCGCCGATACGGTTTTCCGCATTATAGGCCAtt                                   >  2:450241/1‑250 (MQ=255)
                                                                                                                                                                                     gAGCCGCTCCAGCTGGCTCTTGCCGTGCAGGATGGAGCCAGCTATGACGGCGCAATTCTGAAATATGAGTTCAAACAGGCGGGCATTTCCTGGAGAGGGACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGCAGTTGCCCGTAAACAGTCGGCCCCGCTGCCCGCTCCGCTTAAGATTATGCTGAAAAAGACGGACAACATGCTCGCCGCTACGGTTTTCCGCATGATCGGCCATGCGCCCtt                            >  2:41369/1‑251 (MQ=255)
                                                                                                                                                                                      aGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTc                           <  1:262052/251‑1 (MQ=255)
                                                                                                                                                                                         cGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTTGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGCTACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAt                        >  2:346851/1‑251 (MQ=255)
                                                                                                                                                                                            tCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGGGGAACACTGCTGCGCCCGAATCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCCCTGCCCGATCTGCTTAAGATTTTGCTGAAAAAGTCGGACAACATTATCGCCGAAACGGTTTTTCGCCTTGTAAGCCCAGCGCGCTTCAAtgtg                     >  2:327494/1‑251 (MQ=255)
                                                                                                                                                                                             cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAAtgtg                     <  2:365097/250‑1 (MQ=255)
                                                                                                                                                                                               cGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCt                  <  1:346851/251‑1 (MQ=255)
                                                                                                                                                                                                         tGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGGACATgg         <  1:131875/250‑1 (MQ=255)
                                                                                                                                                                                                               gCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATTCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTGCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCg  <  1:444517/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGGCAACATGATCGCCGATACGGTTTTCCGCATGATa                                          >  2:86806/1‑200 (MQ=255)
                                                                                                                                                                                                                          ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATa                                          <  1:86806/200‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
AGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCG  >  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: