Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,923 C→T G277G (GGG→GGA nlpD ← murein hydrolase activator NlpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,9230CT94.4% 92.8 / NA 37G277G (GGG→GGAnlpDmurein hydrolase activator NlpD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (17/17);  minor base G (0/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATT  >  NC_000913/2867688‑2868168
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
aaCTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCTCAGAAGCGCCa                                                                                                                                                                                                                                          <  1:273971/251‑1 (MQ=255)
                ccccATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACtt                                                                                                                                                                                                                          <  1:325578/250‑1 (MQ=255)
                 cccATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCACTCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACttt                                                                                                                                                                                                                         >  2:221705/1‑251 (MQ=255)
                     tGGTCGCTAATTTTTGCCCCGCCTTAACTTATTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTCTGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGCTAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCATCGGGTTCGATCACTTTGccc                                                                                                                                                                                                                     >  2:56233/1‑251 (MQ=255)
                            tATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCg                                                                                                                                                                                                              >  1:102915/1‑251 (MQ=255)
                                  tttccccTCCTTAACTTCTTGTTGTTCCCGGACCAGCTTTGTGTCGTTATGGGCGTAGGCGCTCAGGTAATCATCAGTATGTTTGATGATAATCAGATTACCTTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTTCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCGTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTAGCCCTCATTCGGCCaa                                                                                                                                                                                                         <  2:311352/247‑2 (MQ=255)
                                                  tcttGTTGTTCTCGGACCAGCGTTGTGTCGTTATGTTCTGATGCACTCAGGTATTTATCATTATTTTTGATGATTATCTGATTACCGTAGCCGCGCACCTCTTTTCCAGCATATACAACGCGGCCAGCTGCGTTCTCGATACTTGCCTTTTCTTTGCTGCCTGCGATATCAATCCCCTTGTTGTCTCCCTCTGAGGCGCCGAAGTTTTCTATCACGTTGCCCTGAGTCGGCCAGCGCCATGTCGAGATat                                                                                                                                                                                         <  2:278586/250‑2 (MQ=255)
                                                  tcttGTTGGTCCCGGACCAGCATTGTGTCGTTATGGGCTTAGGCACTCAGGAAATCATCATAATGTTTGCTGATACTCAGCTTTCCGTAGCCGAGCAGCTCGTTACCAGCATAAACAACGCGGCCATCTACGTTCGCGATAATTGCCTGGCCTTTGCTTCCTGCGATATCACTCCCCTTGTTCCCTCCCTCAGAAGCGCAAAAGGATTCGTTCACTTTGCCCTCGGTCGGCCAGCGCCGGGAGGAGAtgt                                                                                                                                                                                         >  2:236464/1‑248 (MQ=255)
                                                   cttGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGc                                                                                                                                                                                      >  1:7442/1‑251 (MQ=255)
                                                             ccGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTg                                                                                                                                                                           <  1:70315/251‑1 (MQ=255)
                                                              cGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGCTTg                                                                                                                                                                           >  1:344950/1‑250 (MQ=255)
                                                                          ggggcGTTATGGGCGTAGGCACTCAGGTAGTCATCATTATGTTTGATGATAATCTGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCGGCCTGCGATATCAAGCCCCTTGTTGCCGCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGGCAGCGCCAGGGGGAGATAGGCGTACGGGTTTATGTACTGCTGAc                                                                                                                                                              <  2:205280/247‑1 (MQ=255)
                                                                                  aTGGGCGTAGGCACACAGGTAATCCTCATTATGTTTGATGCTAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCACCGCCAAGAGCAGAAAGACGTACTGGGTGATCTACTGCTGACAATCGGCt                                                                                                                                                      >  2:377821/1‑251 (MQ=255)
                                                                                        gTAGGCTCTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTTCTGACAGTCGGCTCGGTtg                                                                                                                                                <  2:236538/251‑1 (MQ=255)
                                                                                          aGGCACTAAGGTAATCATCATTATGTTTGCTGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGGGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTtgtg                                                                                                                                              >  2:451673/1‑251 (MQ=255)
                                                                                                       atcatcATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGt                                                                                                                                  >  1:19440/1‑250 (MQ=255)
                                                                                                        acatcatTATGTTTGATGATAATCATATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCTATAATTGCCGGTCCTTTGCGGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCTATCACTTGGCCCTCAGGCGGCCAGCGCCAGGTGGAGATATGCGTATTGGTTTATGTACTACTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTAc                                                                                                                                <  1:292300/250‑1 (MQ=255)
                                                                                                         catcatTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTAcc                                                                                                                               >  2:125686/1‑251 (MQ=255)
                                                                                                            catTATGTTTGATGATAATCAGATTACCGTCGCCGCGCAGTGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGGCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAATCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGGCAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGt                                                                                                                             <  2:386445/250‑1 (MQ=255)
                                                                                                               tATGTTTGATGATAATCAGATTACCGTAGCCGCGCCGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTAc                                                                                                                          >  1:394828/1‑250 (MQ=255)
                                                                                                                  gTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGgc                                                                                                                      <  1:71841/251‑1 (MQ=255)
                                                                                                                                 gATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCGGCGCCAGGTGGAGAGAGGCGTACTGGTTGAGGTACTGCTGACATTCGGCTCGGTTGTGTTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTc                                                                                                       <  1:64544/251‑1 (MQ=255)
                                                                                                                                 gATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTc                                                                                                       >  1:396010/1‑251 (MQ=255)
                                                                                                                                     cccGTAGTCGCTCAGCGCGTTACCAGCATAACCAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTTCCTCCCTCAGAAGCGCCAAAGGTTTCGATCGCTTTGCCCTCTGTCTGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGCCAGTCGGCTCGGTTGTGCTTGCTGTTCGTACCGTTCCAGGCGCTGTGACCGTGGTCGCAg                                                                                                   <  2:373431/250‑1 (MQ=255)
                                                                                                                                          aGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTgg                                                                                               <  2:102915/250‑1 (MQ=255)
                                                                                                                                                gcaacgcgtaaccgtcATAAACCACGCGGCCATCTGCGGTCGCGCTACTTGCCTGTACTTTGCAGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGACGCGCCAAAGGGTTCGATCACTTTGCCCTCAGTCGCCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGCTGTACTGCTGACAGTCGGCTCGGTTGCGCTTTCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCttgttg                                                                                        >  2:76662/16‑251 (MQ=255)
                                                                                                                                                 cAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTtgttgt                                                                                       <  1:321808/251‑1 (MQ=255)
                                                                                                                                                             aGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGATTCGATCCCTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACAttt                                                                           >  1:310796/1‑251 (MQ=255)
                                                                                                                                                                           ggCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCGCTCTGt                                                             <  1:126279/251‑1 (MQ=255)
                                                                                                                                                                                   gCGGTCGCGATAAGTGCCTGTCCTTGGCTTCCTGCGATATCAATCCCCTTGTTGCCGCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCGCTGTCGGCCAGCGCCTGGTGGAGTGGTGCTTACGGTGTGATGTACGGCTGACATTCGGCTCGGTTGTGCTTGCTTTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGTGTGTTGTTCTGCAACAGTTTGTTAGCACTCTTGTCACCCg                                                      <  2:168646/250‑1 (MQ=255)
                                                                                                                                                                                     ggtagcgATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCCTGTTGCCTCCCTCAGAAGCGACCAAGGTGTCGATCACTTTGCCCTCAGGCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGCTGTACTGCTGCCAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGGTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGa                                                     >  2:434121/5‑249 (MQ=255)
                                                                                                                                                                                                  gCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACg                                       <  2:19440/250‑1 (MQ=255)
                                                                                                                                                                                                     tGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTAGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAAt                                   >  1:439227/1‑251 (MQ=255)
                                                                                                                                                                                                                   tgcGATATCAATCCCCTTGTTGCCTCCCTCCGAAGCGCCAACGGTTTCGCTCCCTTTGCCCTCAGTCGGCCAGCGCCAGGAGCAGATAGGCGTACTGGTTGATGTACTGCTGCCAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACg                     >  1:427038/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ccTTGTTGCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTg       <  1:1020/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gttgCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGTCTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCATCGGTGGaa    <  1:197189/250‑1 (MQ=255)
                                                                                                                                                                                                                                       tgCCTCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAAtt  <  2:384972/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
AACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGATAATCAGATTACCGTAGCCGCGCAGCGCGTTACCAGCATAAACAACGCGGCCATCTGCGGTCGCGATAATTGCCTGTCCTTTGCTGCCTGCGATATCAATCCCCTTGTTGCCCCCCTCAGAAGCGCCAAAGGTTTCGATCACTTTGCCCTCAGTCGGCCAGCGCCAGGTGGAGATAGGCGTACTGGTTGATGTACTGCTGACAGTCGGCTCGGTTGTGCTTGCTGTTGGTACCGTTACAGGCGCTGTGACCGTGGTCGCAGTTGGCTTGTTGTTCGGCAACATTTTGTTAGCACTCTGTTCACCCGAAGACTCAGAATACGTAATTGTCGGTTGCGACGCAACAGCAACGGTGGAATT  >  NC_000913/2867688‑2868168

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: