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.7% 212.7 / NA 62intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (27/32);  minor base T (1/1);  total (29/33)
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

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCT  >  NC_000913/3560206‑3560697
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
tcttctTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCTAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCTGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCAGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGTGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cg                                                                                                                                                                                                                                                     >  1:19720/1‑250 (MQ=255)
     ttttACCATCAACGTGTTAGCGACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGGTTTAAC‑CGAACTGGAAGGCAGCGCCGATCCGCCCTGATGG‑cgcagg                                                                                                                                                                                                                                                 >  2:752095/1‑245 (MQ=255)
     ttttACCATCAACGCGTTAGCAACGTTGAGCTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGCTGGCCTTGCGATCGTGCCACGCCGCGTTAAC‑CGACCTGGACGGCAGCGCCGCCCCGCCATGATGG‑cgaagg                                                                                                                                                                                                                                                 >  2:100995/1‑245 (MQ=255)
        tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                              >  1:576838/1‑250 (MQ=255)
            aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtt                                                                                                                                                                                                                                         >  1:241338/1‑251 (MQ=255)
              cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCATGtttt                                                                                                                                                                                                                                       >  1:799037/1‑251 (MQ=255)
              cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtttt                                                                                                                                                                                                                                       >  1:801135/1‑251 (MQ=255)
                         cAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                                                                                                                            >  1:833035/1‑251 (MQ=255)
                            cGTTGAGATTATTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTACGGCGTGGTGCCGATATCGAAAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTACCACGGCGTGTTAAC‑CGCACTGGAAGGAAGCGCCGCACCGCCATGATGG‑CGCAGGATAAGGTTTTGATCTGCCAGCTcc                                                                                                                                                                                                                         >  2:679892/1‑250 (MQ=255)
                                      gttgGTGACAATGCGAAAATTGGTGTGATTGGGTGGTGCTGGCGCTCCCTCTCCGGGGGGGGTGCCTATATCGATAAACGGCGGCGAGCCATTGGGGATTTGCCCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                                                                                                                               <  1:249545/251‑1 (MQ=255)
                                        tgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                             >  1:129961/1‑251 (MQ=255)
                                              cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGACCATTGTGGTCGCGGCGAAt                                                                                                                                                                                                        >  2:770689/1‑250 (MQ=255)
                                                   cgcAAATTGCTGAGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCTAGGTTCCGATATCGATACACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGTGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCATTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCCGGTTTTGCTCCGCCAGCTCATTGAGGTCGAGGCGAATAGTCtc                                                                                                                                                                                                  >  1:677908/1‑250 (MQ=255)
                                                      aaaTTGTTGTGAGTGAGAAGGGCGTGCGTTACCGCTCCCGGGGTGTTGCCGATATAGATAAACAGCGTCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTAAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                               <  1:605774/251‑1 (MQ=255)
                                                          tGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                                           >  2:462706/1‑251 (MQ=255)
                                                           gctgtgtgtGAGCGGTGCGTGCGCGTCCGTGTCCGGGGGGGGGCCGATATCGATAAACAGGGTCGAGCCATGGGGGATTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑a                                                                                                                                                                                         <  1:466549/243‑1 (MQ=255)
                                                            cTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑Ac                                                                                                                                                                                        >  1:243882/1‑251 (MQ=255)
                                                                    gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAAt                                                                                                                                                                                >  1:89549/1‑251 (MQ=255)
                                                                                  cTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGctct                                                                                                                                                                  >  1:103102/1‑251 (MQ=255)
                                                                                     ccGCTTCGGGGGGGGTGCGGATATCAATAAACGGCGTGGGGCCATTGGGGTTTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGTGTGTTAAC‑CGAGCGGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCg                                                                                                                                                               <  1:100995/251‑1 (MQ=255)
                                                                                             ggCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGGTGCGCGCTATTTCTTCGGTCTGGGTGGCCCAGCGCTCGTACAACGGCGTGTTAAC‑CGAACTGGAAGGAAGCGCCGAACCGCCATGATGG‑CTCAGGATCACGTATTGCTCCGCCAGCCCTTTGAGGAATCGGCGAATAGTCTTCGG‑CAGTCCGGAGAAATGctcttcgggatctttggaactacc                                                                                                                                                       >  2:244001/1‑231 (MQ=255)
                                                                                                gtggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                     >  1:52915/1‑250 (MQ=255)
                                                                                                   gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGCAGGCAGCGCCGCCCCGCCATGAGGGAG‑CAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTc                                                                                                                                                                                                    >  2:711600/1‑201 (MQ=255)
                                                                                                               aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTAATGCCATAACCCTGCTGTTTAAc                                                                                                                                     >  1:339337/1‑251 (MQ=255)
                                                                                                                    cacctgcgatccATGGGGGTTTGGCTCCTCCACTTTGCGGTCGTTGCGCTCTTGTTTTTCGTTCTGGTTGGCTTTGCGTTCGTGCCCCGGCTTGTTACC‑CGAATTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTG‑ACTGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                                <  1:75527/241‑1 (MQ=255)
                                                                                                                       cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                              >  1:81365/1‑250 (MQ=255)
                                                                                                                           gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                               >  1:245677/1‑182 (MQ=255)
                                                                                                                            ggcctttGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                                                                                                                        <  1:671708/246‑1 (MQ=255)
                                                                                                                                  gggggATTTGCTCCGCCACTTTGCGGGCGGTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                                                                                                                   <  1:93446/249‑1 (MQ=255)
                                                                                                                                   ggggTTTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt                                                                                                                 <  1:838780/251‑1 (MQ=255)
                                                                                                                                    gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCACCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTTTACCAGCTCTTCGGTACTGACCTAACCCTGCTGTTTAACCAGTTCGATAATAcctctt                                                                                                                  >  2:222002/1‑245 (MQ=255)
                                                                                                                                                 ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTCGtgtttgt                                                                                                   >  2:401666/1‑251 (MQ=255)
                                                                                                                                                      tggcggggggtGCGCTCTTTTGCCGCGGGCGGGGGGGCCTTGCGATGGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCCCCCCACCGCCATGATGG‑CGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                                               <  2:432085/240‑1 (MQ=255)
                                                                                                                                                       tGCGGGCGTGGCGTTCTTTTTCTGCGGTCTGGGTGGCTGTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAGGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAt                                                                                              <  2:801135/250‑1 (MQ=255)
                                                                                                                                                        gggggCGATGCGCGTTTTTTTTTGGGTTGGGGTGGCCTGGCTATCGGGCCCGGGCGGGTTACC‑CGAAGTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                            <  2:243882/249‑1 (MQ=255)
                                                                                                                                                                        ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGGTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                            >  1:232623/1‑251 (MQ=255)
                                                                                                                                                                        ttttCTTCGGTCTGGGTGACCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCTAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCCGCTCTTCGGTACTGACCTAAACCTGCTGTTTAACCCGTTCGATAATACCATCGGGACGTTGGGTTTTTTTCATTTATAAATCCCTGGCatt                                                                            >  2:482278/1‑251 (MQ=255)
                                                                                                                                                                              ccGGTCTGGGTGGCCGTGCGATCGTGCCCCGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttc                                                                      <  1:835331/250‑1 (MQ=255)
                                                                                                                                                                                ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                    >  1:598334/1‑251 (MQ=255)
                                                                                                                                                                                 gTCTGGGTGGCCTGGCGATCGGGCCACGGGGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtt                                                                   <  1:770689/251‑1 (MQ=255)
                                                                                                                                                                                  tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ATGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                                  >  1:128799/1‑251 (MQ=255)
                                                                                                                                                                                     gggTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                                  <  1:542327/248‑1 (MQ=255)
                                                                                                                                                                                          gCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                          <  1:752095/251‑1 (MQ=255)
                                                                                                                                                                                             gtGCGATCGTGCCACGGGGGTTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTg                                                        <  2:241338/249‑1 (MQ=255)
                                                                                                                                                                                               gcaaTCGTGCATCGGTGTGTTAAC‑CGATCTGGAAGGCAGCGCCCCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                                     <  1:862867/248‑1 (MQ=255)
                                                                                                                                                                                                  aTCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGaa                                                  <  1:15727/251‑1 (MQ=255)
                                                                                                                                                                                                  aTCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGa                                                   >  1:56896/1‑250 (MQ=255)
                                                                                                                                                                                                   tCGTGCCACGGCTGGTTAAC‑CGACCGGTAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAAGACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGaa                                                  <  1:461458/250‑1 (MQ=255)
                                                                                                                                                                                                        cccgggCGTTTTAAC‑CGACGTGGAAGCGAGCGCCTCCACGACATGATGG‑CTCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                            <  2:667340/247‑1 (MQ=255)
                                                                                                                                                                                                        ccccGGCGTGTTAAC‑CAAACTGGAGGGCAGCGCAGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                            <  1:658636/248‑1 (MQ=255)
                                                                                                                                                                                                        ccACGGGGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                            <  1:462706/251‑1 (MQ=255)
                                                                                                                                                                                                             ggtggtTAAC‑CGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGc                                       <  2:81365/247‑1 (MQ=255)
                                                                                                                                                                                                             gggtgtTAACACGAAATGGAAGGCGACGCCCCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCGCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAg                                        <  2:755724/249‑1 (MQ=255)
                                                                                                                                                                                                                   tAAC‑CGAACTGGAAGGCCGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCTAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTCGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCg                                  <  1:482278/250‑1 (MQ=255)
                                                                                                                                                                                                                     aacctaccgggAAGGCAGCGCCCCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                                   <  2:232623/239‑1 (MQ=255)
                                                                                                                                                                                                                        ggaacgggagggcggccccgcTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaaa                             <  2:103102/234‑1 (MQ=255)
                                                                                                                                                                                                                         gAACTGGAAGGCCGCGCGGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCTAATCAACAAAAGCCATCGCTAAAc                            <  1:223044/251‑1 (MQ=255)
                                                                                                                                                                                                                          accggggggccAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTGGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaaa                             <  2:50902/239‑1 (MQ=255)
                                                                                                                                                                                                                           aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc                          >  1:307751/1‑251 (MQ=255)
                                                                                                                                                                                                                                aaGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGc                     >  2:822704/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cgccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg              <  1:847628/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ccgacccGCCATGATGG‑CGCAGGATCAGGTGTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGc            <  2:52915/246‑1 (MQ=255)
                                                                                                                                                                                                                                            cACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAt         >  1:449846/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtg      <  2:860001/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    tgatgG‑CGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCAGTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCt  <  2:716413/250‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAC‑CGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG‑ACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCT  >  NC_000913/3560206‑3560697

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