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.G96.6% 213.9 / ‑0.6 59intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (31/26);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑C‑C‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCG  >  NC_000913/3560209‑3560685
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       
tcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑cgc                                                                                                                                                                                                                                        >  2:42743/1‑247 (MQ=255)
     tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGagc                                                                                                                                                                                                                                 >  2:356295/1‑248 (MQ=255)
      aCCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATc                                                                                                                                                                                                                                 >  2:273205/1‑249 (MQ=255)
           catcGTGTTAGCAACGGGGAGTTTGTGGGTGACAATGCGCAAATTG‑CTGTGATGGAGCAGTGCGTGCGCGACC‑GCTTCGGGCGTGGTGCC‑AATATCGATAAACAGCGTCGAGCGTTTGGGGATTTGCTCC‑GCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGCCATCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCATGttt                                                                                                                                                                                                                           <  1:248821/247‑1 (MQ=255)
           cAACGTGTTAGCAACGTTGAGATTGTTGGTGAAAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGtttt                                                                                                                                                                                                                          >  1:162311/1‑251 (MQ=255)
                  ttAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGc                                                                                                                                                                                                                   >  1:233286/1‑251 (MQ=255)
                     gCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                  >  2:425104/1‑249 (MQ=255)
                      cAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGc                                                                                                                                                                                                               >  1:245770/1‑251 (MQ=255)
                          ggtGAGACTGTTGGAGCCAATGCCCAAATTG‑ATGCGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAGCAGCGTCGAGCCATTGGGGCTTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCCTTTTCTGCGGCCTGGGAGGCAGTGCGATCGTGCCACGGCGTGTTAACACAACTGGAAGGAAGCGACCG‑ACCGCAATGAAGG‑C‑CCAGCATCAGGTGTTGCgccgcccgcccat                                                                                                                                                                                                           >  2:174644/3‑244 (MQ=255)
                              agaTTGTGGGTGACATTGCG‑AAATTGGCTGTGATTGAGCAGTGCGTGCGGTACC‑GCTTCCGGCGTGGTGCG‑GATTTCGAAAAACGGCGTCGGGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                                                                                                                       <  1:18464/251‑1 (MQ=255)
                                      ggTGACAAGGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTCCC‑GCTGCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGACCCATTGGGGATTTGCTCC‑GCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                                                                                                               <  1:215393/251‑1 (MQ=255)
                                       gTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                              >  1:424159/1‑251 (MQ=255)
                                        gtgaaatggcgaaattggttgggATTGAGCAGGGCGTGCGCTACC‑GCTTCCGGGGTGTTGCG‑TATATCGAAAAACGGCGTCGTGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                              <  2:343891/229‑1 (MQ=255)
                                          aCAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                           >  1:226469/1‑251 (MQ=255)
                                                    aaTTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                 >  2:299362/1‑251 (MQ=255)
                                                    aaTTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                 >  2:191574/1‑251 (MQ=255)
                                                          gtgtgATTGAGAATGGCGTGCGCTACC‑GCTTCCGGCGGGGGGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                                                                                                            <  1:114072/250‑1 (MQ=255)
                                                              gATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                          >  1:397620/1‑249 (MQ=255)
                                                                  gAGCGGTGCGTGC‑CTACCCGCTTCCGGGGTGTTGCC‑AATATCAAAAAACAGCGTGGAGCCATGGGGGGTTTGCTCC‑GCCTCGTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCGTGCGATCGGGCCACGGCGTGTTAACGGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                    <  1:359324/251‑1 (MQ=255)
                                                                  gAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                  >  2:20873/1‑205 (MQ=255)
                                                                  gAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                    >  2:100109/1‑251 (MQ=255)
                                                                  gAGCAGTGCGTCCGCTCCC‑GCTCCCGGCGTGGTGCC‑GATATGGATAACCAGCGTCGAGCCTTTGGGGATTTCCTCC‑GCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                  <  1:20873/205‑1 (MQ=255)
                                                                      aGTGCGTGAGCTACA‑GCTTCCGGAGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCGGCCGGAGAAATGctc                                                                                                                                                                >  1:198109/1‑251 (MQ=255)
                                                                             gcgcgccc‑GCTTCCGGCGTGGTGCC‑GATATCGATAACCAGCGCCGAGCCATGGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                         <  1:358658/245‑1 (MQ=255)
                                                                                 tACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                     >  2:235373/1‑251 (MQ=255)
                                                                                      gttccgggggtggtgCC‑GATATCGATAAACAGCGTCGACCCTTGGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                 <  1:223856/242‑1 (MQ=255)
                                                                                      ccTTCCGGCGTG‑TGCCCGTTTTCGATAAACAGCGGCGACCCATGGGGGTTTTGCCCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTACCCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                 <  1:261198/250‑1 (MQ=255)
                                                                                       cttaaggCGTGATGCG‑GATATCGATAAACAGCGACGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGAGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTACCCGAACTGGCAGGCAGCG‑CCGCACCACCATGATGGAG‑C‑AGGATCAGCTTTTGCTCCGCCAGCTCTTTGAGGTCGCGGCGAATAAGCTGCGGGCCGACGGAGAAATACTCTACTAACTCTTCggc                                                                                                                                                 >  2:300391/6‑249 (MQ=255)
                                                                                              cGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                         >  1:404202/1‑251 (MQ=255)
                                                                                                     cc‑GATATCGATAAACGGCGTCGAGCCATGGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                                                                                                                  <  1:112352/251‑1 (MQ=255)
                                                                                                     cc‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                                                                                                                  >  1:87591/1‑251 (MQ=255)
                                                                                                        gATATCAATAACCAGCGTCGACCCATTGGGGTTTTGCTCC‑GCCACTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                <  1:313006/251‑1 (MQ=255)
                                                                                                               aTAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                                                                                                         >  1:114356/1‑251 (MQ=255)
                                                                                                                 aaaCAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                                                                                                        >  2:268489/1‑250 (MQ=255)
                                                                                                                    cAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                    >  1:411821/1‑251 (MQ=255)
                                                                                                                    cAGCGTCGACCCATGGGGGATTTGCCCC‑GCCCTTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCGGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                    <  2:92250/251‑1 (MQ=255)
                                                                                                                       cGTCGGGCCATTGGGTATTTGCCCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                                                                                                                 <  1:395490/251‑1 (MQ=255)
                                                                                                                             gCCATTGGGGTTTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                           <  1:19537/251‑1 (MQ=255)
                                                                                                                                aTTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGCTGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                        >  1:358868/1‑251 (MQ=255)
                                                                                                                                 gtggggtttggcccccgccgttgtGCGGGTGATGCGCTCTTTTTTTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                          <  2:405235/226‑1 (MQ=255)
                                                                                                                                  gtgggATTTTCTCCCGCCCCTTGGGGGGCGATGCGTTTTTTTTCGTGGGTCGGGGTGGCCTTGGGATCGGGCCACGGGGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                                                                                                       <  2:157973/249‑1 (MQ=255)
                                                                                                                                            cccc‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                            <  1:268489/249‑1 (MQ=255)
                                                                                                                                                     cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                    >  2:69811/1‑251 (MQ=255)
                                                                                                                                                      tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                    >  2:370094/1‑250 (MQ=255)
                                                                                                                                                         gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                >  1:262418/1‑251 (MQ=255)
                                                                                                                                                                 ggcgcTCTTTTTCTGCGGTCGGGGTGGCCTTGCGACCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCCCCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                        <  2:300692/250‑1 (MQ=255)
                                                                                                                                                                     ctctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTgg                                                                    >  1:215177/1‑251 (MQ=255)
                                                                                                                                                                                          ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                >  2:144936/1‑250 (MQ=255)
                                                                                                                                                                                                   aTCGTGCCACGGCGTGTTAACCGAACGGGAAGGCGGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGaa                                      <  2:222747/251‑1 (MQ=255)
                                                                                                                                                                                                    tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                     >  1:204712/1‑251 (MQ=255)
                                                                                                                                                                                                      gTGCCACGGCGCGTTAACCGAACTGTAAGGCAGCG‑CCCCACCGCCATGATGG‑CGC‑AGGATCAGGTCTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTCACCAGTTCGATAATACCGTTGTGACGTTGTGCTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                    >  1:114326/1‑250 (MQ=255)
                                                                                                                                                                                                            gggCGTTTTAACAGAACGGGAGGGCGGCG‑CCGCTCCGCCGTGTTGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                             <  2:389143/250‑1 (MQ=255)
                                                                                                                                                                                                                  gtTAACCGAACTGGAAGGCAGCG‑ACGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                       <  1:31578/251‑1 (MQ=255)
                                                                                                                                                                                                                     acccGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGc                     <  1:146296/248‑1 (MQ=255)
                                                                                                                                                                                                                          aaCTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATACCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcc               >  1:416777/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggAAGGCGGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg           <  2:245770/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggAAGGCAGCG‑CCGCACCGCCATGCTGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg           >  1:398725/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gggCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt        <  1:34296/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGCAGCG‑CCGCACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt        <  1:345789/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cg‑ccgcACCGCCATGATGG‑CGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg  >  2:408918/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       
TCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACC‑GCTTCCGGCGTGGTGCC‑GATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCGCACCGCCATGATGG‑C‑C‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCG  >  NC_000913/3560209‑3560685

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: ‑