Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,455 +G intergenic (+589/+167) rtcR → / ← glpG DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G94.4% 102.8 / 0.0 36intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (14/20);  total (15/21)
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.89e-01

ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAC  >  NC_000913/3560205‑3560682
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
aTCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                                                                                                                                                                                      >  1:170102/1‑251 (MQ=255)
                          cAACGTTGAGATTGTTGGTGACAATGCGCCCATTGCTGAGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGAGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGATCCGccggc                                                                                                                                                                                                            >  2:53129/1‑248 (MQ=255)
                                  agatggtgggggACAATGCGAAAATTGCGGTGATTGAGCAGTGGGTGCGCTACCGCTTCCGGGGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                                                                                                                    <  1:242105/240‑1 (MQ=255)
                                    aTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAggg                                                                                                                                                                                                  >  2:22876/1‑250 (MQ=255)
                                       gtgggtaaCAATGCGAAAATTGCTGTGATTGAGCATTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGTGATTTGCTCCGCCACTTTGCGGGCGATTCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGTCGTGTTAACCGAACTGTAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                                                                                                                <  1:326121/243‑1 (MQ=255)
                                             gACAATGCGCAATTTGCTGTGAGTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGTTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                                                                                                         <  1:12489/251‑1 (MQ=255)
                                                aaTGCGGAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGCGCACTTGGGGATTTGCTCCGCCACTTCGCGGGCGATGCGCTCTTTTTATTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCGCCCAGCTCAATGAGGTCGCGGCGAAt                                                                                                                                                                                        >  2:286670/1‑249 (MQ=255)
                                                           tGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATAGATAAACTGCGGCGAGCCATTGGGGCGTTGCTCCGCCACTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTAGCCTTGCGATCGTGTCACGGCGTGTTCACCGGACTGGAAGGCCGCGCCGTACCACCATGGTGGCGCAGGATCAGGTTTCGCTCCGCCAGCTCATTCAGGACGCGGCGGTTAGTCTGCTGGCtt                                                                                                                                                                           >  2:55272/1‑250 (MQ=255)
                                                                 gatggaggtgtgcggcctctcccGCTTCCGCGGTGTTGCCGTTATCGATGAACGGGGTCTAGCCTTTGGGGATTTTCCCCGCTACTTGGGTGGCGATGTGCTCTTTTTCTTCGGTCTGGGGTGCCTTGCGATCGGGCCTCGGCGTGTAACTCGAACTGCAAGGCATTGCCTCACCGCCATGATGGTGCAGGATCAGGTGTTGCTCCGCCAGCTCATTGATGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                     <  1:360051/230‑1 (MQ=255)
                                                                   ttGAGCAGTGCGTACGCTACTCCTTCAGGCGAAGTCCCGCGATCACTAAACAGCGTCGAGCCATTCGGGCTCTGCTCCGCCACTTTGCGAGCGATGCGCTCTTTTTCTTCGGTCTGGGTGACCCTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGTGAGCCCCCTACCGCCATGATAGCCCAGCCTCAGTTGTTGCTCCGCCAGCTCACTGAGGTCGCAGCGAATCGACTACGGGCTGCCggggc                                                                                                                                                                    >  2:364576/1‑247 (MQ=255)
                                                                     gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGAGAAAt                                                                                                                                                                 >  1:122587/1‑251 (MQ=255)
                                                                              gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                         >  1:213597/1‑250 (MQ=255)
                                                                                         cTTCCGGAGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGGTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTAGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGt                                                                                                                                              >  1:22445/1‑250 (MQ=255)
                                                                                                   ggtgCCGATATAGATAAACAGCGTCGAGCCATTGGGGATTCGATCAGCCACTTTGCGGGCGCAGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGAAGAAATGCTCTACCAACTCTTCGGTCCTGCCATAAc                                                                                                                                   >  2:28062/1‑251 (MQ=255)
                                                                                                      gccgatatgaaaaaccggCGTCGGGCAATTGGGGATTGCCTCCGCCACTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                <  1:83574/234‑1 (MQ=255)
                                                                                                          atatCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                            >  1:139521/1‑251 (MQ=255)
                                                                                                          ataccGATAACCAGCGGCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                            <  1:94499/247‑1 (MQ=255)
                                                                                                                        gggcggcccATGGGGGTTTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                                                                                                              <  2:272270/244‑1 (MQ=255)
                                                                                                                           cGAGCCATTGGGGATTTGCTCCGCCACTTTGCGAGCGCTGCGCTCTTTTTCTTCGGTGTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGAAGCGCCGCAACGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                            >  1:209260/1‑250 (MQ=255)
                                                                                                                               ccATTGGGGATTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                       <  1:120964/251‑1 (MQ=255)
                                                                                                                                ccgtgtggATTTGCCCCGGCCCCTTTGCGGGCGATGGGCTCTTTTTCTCCGGGCTGGGTGGCCGTGCGATCTGGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGGTCGATAATAcc                                                                                                       <  2:47657/245‑1 (MQ=255)
                                                                                                                                cATTGGGTTTTTGCCCCGCGCCTTTGGGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCGTGCGCTCGTGCCACGGCGTGTTACCCGAACGGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                      <  2:80219/251‑1 (MQ=255)
                                                                                                                                       gATTTGCTCCGCCAATTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCCGCTCTTCGGTTCTGCCCTACACCTGCTGTTTCACCAGTTCGAACATACCATTGTGAc                                                                                               >  2:292106/1‑251 (MQ=255)
                                                                                                                                         tgtGCTCCGCCACTTTGGGGGTGTGGCGCTCTTTTGCTGCGGGCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGACCGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg                                                                                              <  2:37812/248‑1 (MQ=255)
                                                                                                                                               ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTAGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCACCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCCGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                         <  1:55272/249‑1 (MQ=255)
                                                                                                                                               ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                         <  1:22876/249‑1 (MQ=255)
                                                                                                                                                      gttggCGGGCGATGCGTTTTTTTTATTCGGTGTGGGTGGCCGTGCGAGTGTGCAACGGCGTGTTAACCGAACGGGAAGCCACCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTc                                                                                 <  2:241930/247‑1 (MQ=255)
                                                                                                                                                       ttggggggctgtGCGCTTTTTGTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACGTGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                 <  2:252146/239‑1 (MQ=255)
                                                                                                                                                        ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGTCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGTCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                                <  1:225998/250‑1 (MQ=255)
                                                                                                                                                         ggCGGGCGCTGCGCTCTTTTCCTTCGTTCTGGGTGGCCTTGCTGTCGTGCCCCGGCGGTTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATTCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt                                                                              <  2:335066/250‑1 (MQ=255)
                                                                                                                                                                   gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTGGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                    >  2:249166/1‑251 (MQ=255)
                                                                                                                                                                                                      cGTGCCACGGCGTGTTACCCGAAGTGGAAGGCAGGGCGGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                 <  2:51103/251‑1 (MQ=255)
                                                                                                                                                                                                        gcccACGGCGTGTTACCCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGAAATTATTTGCGTTTTCGCGCATTGAGCGAATcaa                               <  1:55589/249‑1 (MQ=255)
                                                                                                                                                                                                                   gtTAACCGAGCGGAAAGCCAGCGCCGCTCCGCCATGAGGGCGCAGTATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCGGCGGGCTGACGGAGAAATGCGCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCTCGCATTGAGCGAAGCAACAAAAGCCATc                    <  2:224009/251‑1 (MQ=255)
                                                                                                                                                                                                                    ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCg                   >  1:85547/1‑251 (MQ=255)
                                                                                                                                                                                                                          cGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTTCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGATTTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAGCGCTaaa              <  2:261017/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAc  >  1:129370/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAC  >  NC_000913/3560205‑3560682

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