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.G97.7% 151.7 / ‑4.3 43intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (21/21);  total (22/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 = 9.82e-01

TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTC  >  NC_000913/3560214‑3560698
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                               >  1:335770/1‑250 (MQ=255)
      cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCGTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGATGG‑CGCAGGCCCAGGtttt                                                                                                                                                                                                                                        >  2:75920/1‑251 (MQ=255)
      aaaggttttAGCAAGGTGGAGATTGTTGGGGACAATGCGAAAATTGCTGGGATTGAGCGGGGGGTGCGCTCCCGCTTCCGGGGGGGTGCCGATATCGATAACCAGCGTCGGGCCATGGGGGTTTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtttt                                                                                                                                                                                                                                        <  1:32418/244‑1 (MQ=255)
        aCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGc                                                                                                                                                                                                                                      >  1:362385/1‑251 (MQ=255)
                    cGTTGAGATTGTTTGTGCCAATGCGCACATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTCGTGCCGATATCGATAAACAGCGTCGAGCCATTGGAGATCTGCTCC‑GCCCCTATGCGGGCGCTGCGCTCTTTTACTTCGATCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTATCCGAACTGTCAGGCAGCGCCCCACCGCCATGATGG‑CGCAGGACCAGATTTTTCTCCGcccacacc                                                                                                                                                                                                                          >  2:82502/1‑245 (MQ=255)
                          gaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAACAGCGCGTGCGCTACCGCTTCCGGCGTGGAACCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCA‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑Agg                                                                                                                                                                                                                   >  1:42767/1‑251 (MQ=255)
                                      cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGTG‑CAGGATCAGGTTTTGCTCCGCCAGCT‑ACTGTAGGTCGCGGCGAAt                                                                                                                                                                                                        >  2:100962/1‑250 (MQ=255)
                                                   gCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTT‑AGGTCGCGGCGAATAGTCTGCGTGCTgg                                                                                                                                                                                          >  2:4686/1‑250 (MQ=255)
                                                       tgATTGAGCAGTGCGCGCGCTACCGCTTCCGGCGTGGTGACGATATCGCTACACAGCGTCGAGACATTGGGGATTTGCTCC‑GCAACTTTGCGGGCGATGCGCACTTTTTCTTCGGTATGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAAACGAACTGGAACGCAGCGCCGCACCGCAACGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                       >  1:289564/1‑250 (MQ=255)
                                                          ttGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                                                                                                                                    >  1:192985/1‑250 (MQ=255)
                                                            aagcatggcgtccgctccggCTTCGGGCGTGGTGCGGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCCCC‑GCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                                                                                                                                  <  1:104567/231‑1 (MQ=255)
                                                                         gcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCCGGTTTTGCTCCGCCAGCTCCTTT‑AGGTCGCGGCGAATAGTCTGCGGGATGACGGCGCAATGCTCTACCCGctc                                                                                                                                                                    >  2:266430/1‑251 (MQ=255)
                                                                                               cGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                                                                                                                              >  1:31638/1‑251 (MQ=255)
                                                                                                 aaatggATAAACAGCGTCGAGCCATTGGGTATTTGCTCC‑GCCTCTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                                            <  1:247486/246‑1 (MQ=255)
                                                                                                     cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAACTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                         >  1:197547/1‑250 (MQ=255)
                                                                                                             ggcggcGAGCCATTGGGGTGTTGCTCC‑GCCTCTTGGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                                <  1:266430/246‑1 (MQ=255)
                                                                                                                   gAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTcgc                                                                                                                                                                                                               >  1:381666/1‑166 (MQ=255)
                                                                                                                    gcgcgttGGGGGTTTGCCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCGTCGGTCTGGG‑TGGCCGTGCGACCGTGCCACGGCGTGTTACCCAACCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                          <  2:31638/246‑1 (MQ=255)
                                                                                                                      cctttGGGGTTTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                                       <  1:307613/248‑1 (MQ=255)
                                                                                                                              gattgccTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                                               <  1:215249/245‑1 (MQ=255)
                                                                                                                                          ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgt                                                                                                    >  2:140382/1‑251 (MQ=255)
                                                                                                                                            actttcagggCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCTATGATGG‑CGCAGGATCAGGTTTTGCTCCGCTAGCTCATTG‑AGGTGGCGGCGATTAGGCTGCGGGTTGTCGGAGAAATGCTCTACCAGCTCTCCGTGACTGACATAACCCTGCTGGTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                                  >  2:78238/8‑251 (MQ=255)
                                                                                                                                             cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGGTAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                                 >  2:36874/1‑251 (MQ=255)
                                                                                                                                              ttggCGGGCGAGGGGTTCTTTTCCTTCGGCCGGGG‑GGGCCTTGCGACCGTGCCCGGGGGGGTTAACGGACCGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                                 <  2:462302/247‑1 (MQ=255)
                                                                                                                                                    ggggATGCCCTTTTTTTCTCCGGTCTGGG‑TGGCCTTGCGACCGTGCCGCGGCGTGTTACCCGACCTGGAAGGCAGCGCCCCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                          <  2:358126/248‑1 (MQ=255)
                                                                                                                                                      ggATGCCCTCTTTTTCTTCGGTCGGGG‑TGGCCGTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                                        <  1:179091/250‑1 (MQ=255)
                                                                                                                                                          gcgcTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                     <  1:419568/250‑1 (MQ=255)
                                                                                                                                                          gcgcTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                                    >  1:12697/1‑251 (MQ=255)
                                                                                                                                                                gttttCTTCGGT‑TGGGCGGGCCTTGCGGTCGTGCCACGTGGTTTTAACCGAACGGGAAGGCAGCCCCTCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCATGTCATTG‑AGGTCGCGGCGAACAGGCTGCGGGCTGACGGAGAAATGCTGTACCAGCTCTTCGGTACTGACATATCCCTGCGGTTTAACCAGTTCGATAATACCGGGGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGaa                                                                               <  1:412884/249‑1 (MQ=255)
                                                                                                                                                                  tttCTTCGGTCGGGG‑GGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                             <  1:140382/250‑1 (MQ=255)
                                                                                                                                                                         ggTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                     >  1:103580/1‑251 (MQ=255)
                                                                                                                                                                         ggTCTGGG‑TGCCCGTTCTATCGGGCCACGGCGTGTTATCCGAACTGGAGGGCAGCGCCTCACCGCCATGATTG‑CGCATGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGATTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                     <  1:232137/251‑1 (MQ=255)
                                                                                                                                                                                    gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                           <  1:36874/251‑1 (MQ=255)
                                                                                                                                                                                        tGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTg                                                         <  1:117397/249‑1 (MQ=255)
                                                                                                                                                                                          ggatggggccgaggCGTGTGAGCCGATCTGGAGGGCGGCGCCGCACCGCCATGAGGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGGCGGAGAAATGCGCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGGTGTGTTGGTTTCATTTATAAATCCCTGGAATGATTTTCGTTTTCGCGCATTGAGCg                                                     <  2:42767/239‑1 (MQ=255)
                                                                                                                                                                                           gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                                       >  2:344259/1‑216 (MQ=255)
                                                                                                                                                                                           gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                                       <  1:344259/216‑1 (MQ=255)
                                                                                                                                                                                                 gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaac                                              >  2:367963/1‑251 (MQ=255)
                                                                                                                                                                                                   cggggCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaacaa                                            <  2:362766/248‑1 (MQ=255)
                                                                                                                                                                                                                           ggCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa                    <  1:91161/251‑1 (MQ=255)
                                                                                                                                                                                                                              aGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑CGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAACCCATCGCTAAACCCACGGCTAAcc                 >  2:413299/1‑251 (MQ=255)
                                                                                                                                                                                                                               ccgccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAccc                <  1:72881/250‑1 (MQ=255)
                                                                                                                                                                                                                                             tgatgG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTc  >  1:237860/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGG‑TGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTG‑AGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTC  >  NC_000913/3560214‑3560698

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