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

ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAC  >  NC_000913/3560205‑3560682
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
aTCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑c                                                                                                                                                                                                                                      >  1:170102/1‑251 (MQ=255)
                          cAACGTTGAGATTGTTGGTGACAATGCGCCCATTGCTGAGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGAGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCCGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGATCCGccggc                                                                                                                                                                                                            >  2:53129/1‑248 (MQ=255)
                                  agatggtgggggACAATGCGAAAATTGCGGTGATTGAGCAGTGGGTGCGCTACCGCTTCCGGGGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                                                                                                                    <  1:242105/240‑1 (MQ=255)
                                    aTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAggg                                                                                                                                                                                                  >  2:22876/1‑250 (MQ=255)
                                       gtgggtaaCAATGCGAAAATTGCTGTGATTGAGCATTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGTGATTTGCTCCG‑CCACTTTGCGGGCGATTCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGTCGTGTTAAC‑CGAACTGTAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                                                                                                                <  1:326121/243‑1 (MQ=255)
                                             gACAATGCGCAATTTGCTGTGAGTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGTTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                                                                                                         <  1:12489/251‑1 (MQ=255)
                                                aaTGCGGAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGCGCACTTGGGGATTTGCTCCG‑CCACTTCGCGGGCGATGCGCTCTTTTTATTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCGCCCAGCTCAATGAGGTCGCGGCGAAt                                                                                                                                                                                        >  2:286670/1‑249 (MQ=255)
                                                           tGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATAGATAAACTGCGGCGAGCCATTGGGGCGTTGCTCCG‑CCACTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTAGCCTTGCGATCGTGTCACGGCGTGTTCAC‑CGGACTGGAAGGCCGCGCCGTACCACCATGGT‑GG‑CGCAGGATCAGGTTTCGCTCCGCCAGCTCATTCAGGACGCGGCGGTTAGTCTGCTGGCtt                                                                                                                                                                           >  2:55272/1‑250 (MQ=255)
                                                                 gatggaggtgtgcggcctctcccGCTTCCGCGGTGTTGCCGTTATCGATGAACGGGGTCTAGCCTTTGGGGATTTTCCCCG‑CTACTTGGGTGGCGATGTGCTCTTTTTCTTCGGTCTGGGGTGCCTTGCGATCGGGCCTCGGCGTGT‑AACTCGAACTGCAAGGCATTGCCTCACCGCCATGAT‑GGTG‑CAGGATCAGGTGTTGCTCCGCCAGCTCATTGATGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                     <  1:360051/230‑1 (MQ=255)
                                                                   ttGAGCAGTGCGTACGCTACTCCTTCAGGCGAAGTCCCGCGATCACTAAACAGCGTCGAGCCATTCGGGCTCTGCTCCG‑CCACTTTGCGAGCGATGCGCTCTTTTTCTTCGGTCTGGGTGACCCTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGACGTGAGCCCCCTACCGCCATGATAGC‑C‑CAGCCTCAGTTGTTGCTCCGCCAGCTCACTGAGGTCGCAGCGAATCGACTACGGGCTGCCggggc                                                                                                                                                                    >  2:364576/1‑247 (MQ=255)
                                                                     gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGAGAAAt                                                                                                                                                                 >  1:122587/1‑251 (MQ=255)
                                                                              gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                         >  1:213597/1‑250 (MQ=255)
                                                                                         cTTCCGGAGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGGTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTAGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGt                                                                                                                                              >  1:22445/1‑250 (MQ=255)
                                                                                                   ggtgCCGATATAGATAAACAGCGTCGAGCCATTGGGGATTCGATCAG‑CCACTTTGCGGGCGCAGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGACCTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGAAGAAATGCTCTACCAACTCTTCGGTCCTGCCATAAc                                                                                                                                   >  2:28062/1‑251 (MQ=255)
                                                                                                      gccgatatgaaaaaccggCGTCGGGCAATTGGGGATTGCCTCCG‑CCACTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                <  1:83574/234‑1 (MQ=255)
                                                                                                          atatCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                            >  1:139521/1‑251 (MQ=255)
                                                                                                          ataccGATAACCAGCGGCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                            <  1:94499/247‑1 (MQ=255)
                                                                                                                        gggcggcccATGGGGGTTTTGCTCCG‑CCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                                                                                                              <  2:272270/244‑1 (MQ=255)
                                                                                                                           cGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGAGCGCTGCGCTCTTTTTCTTCGGTGTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGCAGGAAGCGCCGCAACGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                            >  1:209260/1‑250 (MQ=255)
                                                                                                                               ccATTGGGGATTTGCTCCG‑CCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                       <  1:120964/251‑1 (MQ=255)
                                                                                                                                ccgtgtggATTTGCCCCGGCCCCTTTGCGGGCGATGGGCTCTTTTTCTCCGGGCTGGGTGGCCGTGCGATCTGGCCACGGCGTGTTAAC‑CGAACGGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGGTCGATAATAcc                                                                                                       <  2:47657/245‑1 (MQ=255)
                                                                                                                                cATTGGGTTTTTGCCCCG‑CGCCTTTGGGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCGTGCGCTCGTGCCACGGCGTGTTACC‑CGAACGGGAGGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                      <  2:80219/251‑1 (MQ=255)
                                                                                                                                       gATTTGCTCCG‑CCAATTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGCAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCCGCTCTTCGGTTCTGCCCTACACCTGCTGTTTCACCAGTTCGAACATACCATTGTGAc                                                                                               >  2:292106/1‑251 (MQ=255)
                                                                                                                                         tgtGCTCCG‑CCACTTTGGGGGTGTGGCGCTCTTTTGCTGCGGGCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAAC‑CGACCGGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg                                                                                              <  2:37812/248‑1 (MQ=255)
                                                                                                                                               ccG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTAGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCACCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCCGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                         <  1:55272/249‑1 (MQ=255)
                                                                                                                                               ccG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                         <  1:22876/249‑1 (MQ=255)
                                                                                                                                                      gttggCGGGCGATGCGTTTTTTTTATTCGGTGTGGGTGGCCGTGCGAGTGTGCAACGGCGTGTTAAC‑CGAACGGGAAGCCACCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTc                                                                                 <  2:241930/247‑1 (MQ=255)
                                                                                                                                                       ttggggggctgtGCGCTTTTTGTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAAC‑CGAACGTGAAGGCGGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                 <  2:252146/239‑1 (MQ=255)
                                                                                                                                                        ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGTCGATCGTGCCACGGCGGGTTAAC‑CGAACTGGAAGTCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                                <  1:225998/250‑1 (MQ=255)
                                                                                                                                                         ggCGGGCGCTGCGCTCTTTTCCTTCGTTCTGGGTGGCCTTGCTGTCGTGCCCCGGCGGTTTAAC‑CGAACTGGAAGGCAGCGCCGCCCCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATTCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt                                                                              <  2:335066/250‑1 (MQ=255)
                                                                                                                                                                   gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTGGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                    >  2:249166/1‑251 (MQ=255)
                                                                                                                                                                                                      cGTGCCACGGCGTGTTACC‑CGAAGTGGAAGGCAGGGCGGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                 <  2:51103/251‑1 (MQ=255)
                                                                                                                                                                                                        gcccACGGCGTGTTACC‑CGACCTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGAAATTATTTGCGTTTTCGCGCATTGAGCGAATcaa                               <  1:55589/249‑1 (MQ=255)
                                                                                                                                                                                                                   gtTAAC‑CGAGCGGAAAGCCAGCGCCGCTCCGCCATGAG‑GG‑CGCAGTATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCGGCGGGCTGACGGAGAAATGCGCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCTCGCATTGAGCGAAGCAACAAAAGCCATc                    <  2:224009/251‑1 (MQ=255)
                                                                                                                                                                                                                    ttAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCg                   >  1:85547/1‑251 (MQ=255)
                                                                                                                                                                                                                          cGAACTGGAGGGCAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTTCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGATTTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAGCGCTaaa              <  2:261017/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGCGCCGCACCGCCATGAT‑GG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAc  >  1:129370/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAC  >  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: ‑