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.G98.5% 235.9 / ‑4.6 68intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (29/38);  total (29/39)
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.89e-01

TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTG  >  NC_000913/3560208‑3560693
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
ttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                                                                                                               >  1:369455/1‑248 (MQ=255)
ttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                                                                                                               >  1:457707/1‑248 (MQ=255)
  cttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgcagg                                                                                                                                                                                                                                            >  2:569655/1‑246 (MQ=255)
    ttgaccttCAACGTGTTAGCAACGGTGAGATGGTTGGTAACAAAGCGCAATTGGCTGTGATTGAGCAGGGCGGCCGCTCCCGCTTCGGGGGTGGGGCCGATATCGAAAAACAGCGTCGAGCCATGGGGGATTTGGTCCGCCACTTGGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAGGGcgcag                                                                                                                                                                                                                                             <  1:100671/241‑5 (MQ=255)
         cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGt                                                                                                                                                                                                                                     >  2:519049/1‑251 (MQ=255)
         cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                                                                                                                                                                                                     >  1:386656/1‑251 (MQ=255)
             aacgtgtaagaaaggtgggtatttttgGTGAAATGGCGCAATTTGCTGTGATTGGCCAGTGCGGGCTCTCCCTCTCCCGGCGGGTTCCCAATATCGATAAACAGCGGCGCGCCATGGGGGATTTGTTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGTCAGCGCCGCACCGCCATGTTGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                                 <  2:15814/227‑1 (MQ=255)
              aCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGc                                                                                                                                                                                                                                >  2:563827/1‑251 (MQ=255)
                  gttagacaCGTTGAGATAGATGGTGACAATGGGCAAATGGCTGAGCTTGAGCAGTGCGCGCGCTAACGCTTACGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGCTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTAGGACTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCAGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                                                                                                                            >  2:640313/8‑251 (MQ=255)
                        aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                                      >  1:413120/1‑251 (MQ=255)
                                        gggACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTGCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                 <  1:470203/238‑1 (MQ=255)
                                        gTGACAATGCGCAAATTGCTGTGATTGCGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATCTCGATAAAACGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                 >  2:470203/1‑240 (MQ=255)
                                        gTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                                      >  1:187619/1‑251 (MQ=255)
                                                gcgcAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                                                                                                                              >  2:220159/1‑251 (MQ=255)
                                                  gacaatttGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGGGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                                                                                                                                                                             <  1:533876/245‑1 (MQ=255)
                                                   gcAAATTGCTTTGATTGAGCAGGGCGTGCGTTACCGCTTCGGGCGTGGTGCGGATACCGATAAACAGCGTCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                                                                                                                            <  1:610261/251‑1 (MQ=255)
                                                     aaaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGGCGGCCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                          <  1:159614/251‑1 (MQ=255)
                                                               gATTGAGCAGCACGCGCGCCACAGCTTACTGCGTGGAGCCGACCTCGATAACCCGCGTCGAGCCATTGGTGATTTGCTCCGCCACTTTGAAGGCGATGAGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGCAGGCAGCGCCGCACAGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGACGCGGCGAATAGACTGCGGGCTgcaggcg                                                                                                                                                                                >  2:472699/1‑245 (MQ=255)
                                                                 ttGAGCAGTGCGTGCGCTACCGCTTCCGGCGAGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                                                                                                                               >  2:175527/1‑250 (MQ=255)
                                                                       gtggctccccctcccgttacgggtgtggtgCTGATATCTAAAAACAGCGCGGCGCTATTGGGGTTTTTCCCCGCCACTTTGCGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGGGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCGGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                        <  2:60090/227‑1 (MQ=255)
                                                                                 cgcccccTTCCGGCTTGGTGCCAATACCGATAAACGGCGCGGAGCCATGGGGGATTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTCCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctct                                                                                                                                                              <  1:584166/245‑1 (MQ=255)
                                                                                   cccGCTTCCGGGGGGGTGCCAAATTGGATAAACAGGTGCGAGCCTGTGGGTATTTGCTCGGCCGCTTTGCGGGCGATGCGCTCGTTTCCTTCGGTTGGGGGGGCTTTGGTACCGTTCCAGGGCGTGTTAACCGAACGGGAATGCGGGGCCGCTCCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCTAATAGTCTGCGGGCGGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                            <  2:303105/250‑1 (MQ=255)
                                                                                   aCCGCTTCCGGCGTGGGGCCAATATCGATAAACCGCGTCGAGCCAGTGGGGATTTGCTCCGCCACTTTGCGGGCGTTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctct                                                                                                                                                              <  1:5987/249‑1 (MQ=255)
                                                                                    ccgctccgggtgtggtccggATATCGAAAAACAGCGTCGACCCTTTGGGGTTTTGCTCCGCCACTTGGCGGGCGATGCTCTCTTTTTCTGCGGTCTGGGTGGCCTTGCTATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                                                                                                                           <  2:541543/232‑1 (MQ=255)
                                                                                      gCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                         >  1:542315/1‑251 (MQ=255)
                                                                                                          atCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                                     >  2:580081/1‑251 (MQ=255)
                                                                                                             gATAAACAGCGGCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                                                                                                                                  <  1:397975/251‑1 (MQ=255)
                                                                                                                  aacggcgccgggCCATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                              <  1:308595/240‑1 (MQ=255)
                                                                                                                    caccggggaccaagtGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTGTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTAAACCGAACGGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                            <  1:43843/237‑1 (MQ=255)
                                                                                                                      cGTCGGGCCATGGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                                                                                                                         <  1:563827/251‑1 (MQ=255)
                                                                                                                          cGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAAGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACAAGTTCGata                                                                                                                       >  1:213150/1‑250 (MQ=255)
                                                                                                                                     ggtttgcgtccgcccgtttGCGGGGGATGCGGTCTTTTTCTTCGGACGGGGGGGCCGTGCGATCGTGCCACGGCGTGTGAACCGAACGGGAAGGCACCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa                                                                                                           <  2:437168/235‑1 (MQ=255)
                                                                                                                                       aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                     <  1:379631/143‑1 (MQ=255)
                                                                                                                                       aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                     >  2:379631/1‑143 (MQ=255)
                                                                                                                                                gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                <  1:568472/251‑1 (MQ=255)
                                                                                                                                                    cgttgcgggggtggcgcgttttttCTTCGGTCGGGTGGGCGTTGCTCGCGTGCCGCGGCGTTTTATCCGAACTGGAAGGCGTGGCCGTACTGCTATGTTGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATTGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACTCTGGTGTTTAACCAGTTCGATAATACCGTTGGGACGTTGtgtttgtt                                                                                              <  2:373952/230‑1 (MQ=255)
                                                                                                                                                     tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGACACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGATCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGTTCGATAATAACGTTGTGACGTTGTGTTTGTTTcc                                                                                           >  2:474046/1‑250 (MQ=255)
                                                                                                                                                      ttGCGGGGGATGCGCTCTTTTTCTTCGGCCTGGGGGGCCTGGCGACCGGGCCACGGGGTGTTAACCGAACGGGAAGGCAGGGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                            <  1:519049/249‑1 (MQ=255)
                                                                                                                                                      ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAt                                                                                          >  1:433086/1‑251 (MQ=255)
                                                                                                                                                        gggggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                        <  1:460002/249‑1 (MQ=255)
                                                                                                                                                        gctggcgtggctctgtttttCTCGGGTCTGGGTGTCCTGGCGATCGTCCCACGGCGTGGTAAACGACCGGTAAGGCAGGGCGGCCCCGCCATGAGGGCGCAGGATCAGGTGTTGCTCCGCCAGCTCATTGATGTCGCGGCGATTAGTCTGCGGGCTGACGGATAAATGCGCTTCCAGCTCTTCGGTACTGACATTACCCTGCTGTTTAACCATTTCGATAATACCGTTGTGACGTTGTGTTTTTTTCAtg                                                                                         <  2:96541/235‑2 (MQ=255)
                                                                                                                                                        gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCCttt                                                                                        >  2:586128/1‑251 (MQ=255)
                                                                                                                                                        gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt                                                                                         <  1:647085/250‑1 (MQ=255)
                                                                                                                                                          gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                      >  2:69112/1‑251 (MQ=255)
                                                                                                                                                           ggCGATGCGCTTTTTTTCTGCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                     <  2:369455/251‑1 (MQ=255)
                                                                                                                                                            gCGATGCGCTTTTTTTTGTCGGTCTGGTGGGCCTTGCGAGCGGGCCACGGTTGTTTAACCGAATGGGAAGTCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                    <  1:428176/251‑1 (MQ=255)
                                                                                                                                                                tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                >  1:288066/1‑251 (MQ=255)
                                                                                                                                                                   gttttttttCTTGGGTCGGGGTGGCCGTGTGAGCGTGCCACGGCGGGTGACCCTAACGGAAGGGCAGCGCGCCACCGCCATGATGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTg                                                                             <  2:424125/247‑1 (MQ=255)
                                                                                                                                                                    ctctTTTTCTTCGGCCTGGGTGGCCTTGCGATCTTGCCACGGCGGGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                              <  1:64532/249‑1 (MQ=255)
                                                                                                                                                                     tctTTTTATTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGAAGGAGAAATGCTATACACGATCTTCGGTACTGACATAACACTGCTGTTTAAACAGTTCGATAATAACGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTgg                                                                            >  2:90167/1‑250 (MQ=255)
                                                                                                                                                                        ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAACCTGCTGTTTAACCAGTTCGATAAAACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                        >  1:296176/1‑251 (MQ=255)
                                                                                                                                                                          ttcttcGGTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATGGGGGACGCGGCGAATAGTCTGCGGGCTGACGGAGCACTGCTCGACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTAGATAATACAGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGActtat                                                                      >  2:668812/1‑251 (MQ=255)
                                                                                                                                                                           tcttcGGGGTGGGTGGCCGTGCGACCGTGCCACGGCGTGGTAACGTAACGGGTAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                      <  2:177411/250‑1 (MQ=255)
                                                                                                                                                                            cttcGGTCTGGGTGGCCTTGCGATCGGGCAACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTGATAAATCCCTGGAATTAttt                                                                    >  2:571541/1‑251 (MQ=255)
                                                                                                                                                                                ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                <  1:638849/251‑1 (MQ=255)
                                                                                                                                                                                 gTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAGCCCTGGAATTAttttcgtt                                                               <  1:352347/251‑1 (MQ=255)
                                                                                                                                                                                       gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGAAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                         <  1:181021/251‑1 (MQ=255)
                                                                                                                                                                                         ggCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCGGCACCGCCATGATGGCGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCGGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAGCCCTGCTGTGTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                        <  1:472699/250‑1 (MQ=255)
                                                                                                                                                                                          gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                      <  1:66846/251‑1 (MQ=255)
                                                                                                                                                                                                gCGATCGGGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                                 <  1:613722/251‑1 (MQ=255)
                                                                                                                                                                                                   aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGaa                                              <  1:447804/251‑1 (MQ=255)
                                                                                                                                                                                                     cggcccACGGCGTGTTACCCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTGCGTTTTCGCGCATTGAGCGaa                                              <  2:163562/245‑1 (MQ=255)
                                                                                                                                                                                                          cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                        >  1:571380/1‑250 (MQ=255)
                                                                                                                                                                                                             ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                                                                                                                                   >  1:275876/1‑156 (MQ=255)
                                                                                                                                                                                                                   gtgaaCCGAGCTGGATGTCCGCGCCGCTCCGCCATGACTGAGCAGGATCAGGTTTTGCGCCTCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGATAAAAGCTCTTCCAGCTCTTCGGTACTGACATAACCCTTCTGTTTAACCAGTTCGATAATACCGTTTTTACTTTGTGTTTGTTTCATTTTTAAATCCCTGGAATTATTTTCGTTATCGCGCAGTGATCGACTcaacataatccttc                               <  2:406911/248‑10 (MQ=255)
                                                                                                                                                                                                                      atccaaacgggacggCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCGTCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCTCCTTGATCGAATCAACAAAATCCATCGc                             <  2:116707/237‑1 (MQ=255)
                                                                                                                                                                                                                          gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGACTCAACAAAAGCCATCGCTAAAc                        >  1:335619/1‑251 (MQ=255)
                                                                                                                                                                                                                             cgggAAGGCAGCGCAGCACCGCCATGATGGCGCAGGATCAGGTTTTGATCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAGTGAGCGAATCAACAAAAGCCATCGCTAAAcccc                     <  2:168625/249‑2 (MQ=255)
                                                                                                                                                                                                                                 aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTCTTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGc                 >  1:588125/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc            >  2:318905/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAt     <  1:474046/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtg  >  1:128658/1‑251 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTG  >  NC_000913/3560208‑3560693

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