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

CTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560207‑3560686
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       
cttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑cgc                                                                                                                                                                                                                                        >  1:292615/1‑249 (MQ=255)
  tcttTTACCATCAACGTGTAAGCAAGTTTGAAATTGTTGGTGAAAATGCGCAAATTGCGGTGAGTGAGCAGTGCGTGCGCTACCGCTTCCGTCTTTGGCCCAATATGAAAAAACGTCGTCGAGCCA‑TTGGGGATGTG‑CTCCGCCAGTTTCCGGGGGATGCGCTCTTTTTCTTCGGTCGGTGGGGCCTTGCGATCGGGCCACGGTGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑cg                                                                                                                                                                                                                                         <  1:459546/248‑2 (MQ=255)
                             ttGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCAtt                                                                                                                                                                                                           >  2:344260/1‑251 (MQ=255)
                               gggaTTGTGGGGGACAATGCGCAAATTGCTGTGATTGGGAAGTGGGTGCGCTACCGCTCCGGGCGGGGGGCGGATATCGAGAAACGGGGGCGCGCCA‑GGGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCTCCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                                                                                                                         <  1:216218/249‑1 (MQ=255)
                               gagaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                          >  2:354053/1‑250 (MQ=255)
                                ggaTTTTTGGTGACAATGCGCAAATTGCTGGGATGGAGCAGGGCGTGCGCTCCCGTTCCCGGGGGGGTGCCGTTGTCGTTAACCAGCGGCGAGCCA‑TGGGGGATTGG‑CCCCGCCACTGTTCGGGCGATGCGCTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                                                                                                                        <  1:89264/250‑1 (MQ=255)
                                       gggTGACAATGCGCAAATTGCTGTGATTGAGCAGGGCGGGCGCTACCGCTTCCGGCGTGGTTCCGATATCGATAACCGGCGTCGGGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                 <  1:372308/250‑1 (MQ=255)
                                       gggTGAAAATGGGCAAATTGCTGGTTTTGAGCGGGGCGTGCCCTCCCGTTTCGGGCGGGGTGCCGATATCGAAAAACAGCGTCG‑GCCAATTGGGGATTTC‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                                                                                                                  <  1:163200/249‑1 (MQ=255)
                                             cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAAt                                                                                                                                                                                            >  2:240834/1‑250 (MQ=255)
                                              agggcgcAAATTGTTGTGATTGAGCGGGGCGTGCGCTCCCCCTTCCGGCGTGGTCCCGATACCGAAAACCAGCTTCGGGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                                                                                                          <  1:394924/248‑1 (MQ=255)
                                               aTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                                                                                                         >  1:111271/1‑251 (MQ=255)
                                                   gcAAATTGCTGTGATGGAGCAGTGCGTGCGCTCCCGCTCCGGGCGTGTGGCCGATATCGATAACCAGCGTGGAGCCT‑TTGGGGTTTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTGGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGc                                                                                                                                                                                     <  1:164770/251‑1 (MQ=255)
                                                          gCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                               >  1:239983/1‑250 (MQ=255)
                                                           gtgggATTGAGCAGTGCGTGCGCTACCGCTCCGGGCGTGGTGCGGAAATCGATAAACAGCGTCGAGCCT‑TGGGGGTTTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGAc                                                                                                                                                                             <  1:167351/247‑1 (MQ=255)
                                                           cTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                                                                                                             >  2:249139/1‑251 (MQ=255)
                                                                   gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATG‑ATGGAG‑CAGGATCAGGTTTTGCTCCGCAAGCCCATTGCGGTAGCGGCGACTAATCTGCGGGATGACGGAGAAAt                                                                                                                                                                     >  1:180248/1‑251 (MQ=255)
                                                                            gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                             >  1:105972/1‑250 (MQ=255)
                                                                                        ttCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                                                >  2:246094/1‑251 (MQ=255)
                                                                                             gggtggtgCCGATATCGAAAAACAGGGGCGAGCCT‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                                                                                                                           <  1:382377/249‑1 (MQ=255)
                                                                                              cGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                          >  1:470042/1‑251 (MQ=255)
                                                                                                gggggcctataccGATAACCGGCGTGGAGCCA‑TTGGGGTTTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                                                                                                                        <  1:164474/239‑1 (MQ=255)
                                                                                                    gccgaaagagaaaaaccggcgtgAGCCA‑GTGGGGGTTTGTCTCCTCCATTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATGGCGCCCCGGCGTGTTAACCGAACGGGAAGGCAGCCCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                     <  2:75793/229‑1 (MQ=255)
                                                                                                                   cagcgccggtccA‑TTGGG‑ATTTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                     <  1:354053/241‑1 (MQ=255)
                                                                                                                    cggggcggggcatttgggtttttG‑CCCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                     <  1:428756/231‑1 (MQ=255)
                                                                                                                      cgtggAGCCA‑TTGGGTATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                                                                                                                  <  1:17449/247‑1 (MQ=255)
                                                                                                                      cGTCGAGCAA‑TTGGGGATTTG‑ATCAGCCCCTTTGCGGGCGATGCGCTCTTTTTATTCGGTATGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATTCGGTCCTGACATACCCCTGCTGTTTAACCAGTTCg                                                                                                                   >  2:455280/1‑250 (MQ=255)
                                                                                                                       gcgtagagtcattgggtgttgccaccGCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCGTGCGAGCGTGCCCCGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGATAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                    <  2:350468/226‑1 (MQ=255)
                                                                                                                         cGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTATGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGGAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTCCAAGCTCTTCGGTACTGACCTCACCCTGCTGTTTAACCAGTTCGata                                                                                                                >  2:217834/1‑250 (MQ=255)
                                                                                                                          cgggcattgggggTTTTG‑CTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                                                                                                               <  1:199146/242‑1 (MQ=255)
                                                                                                                                              cccGCCACTTTGCGGGCGTTGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCCCGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt                                                                                            <  1:132616/250‑1 (MQ=255)
                                                                                                                                                  ccACGTTGCGGGCGATGCGGTCTTTTTCGTCGGTCTGGGTGGCCTTGCGTGCGTGCCACGTCGTGTAAACCTCACTGGAGGGCAGCCCCGCACCGCCATG‑ATGG‑CGCTGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGATCTGACGGAGAAATGCTCTGCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgt                                                                                        <  1:304441/251‑1 (MQ=255)
                                                                                                                                                   cccTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGGGCCACGGCGTGTTACCCGAACGGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                       <  1:367355/249‑1 (MQ=255)
                                                                                                                                                     cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                     >  1:369367/1‑251 (MQ=255)
                                                                                                                                                      tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGTCGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATAAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGGTCGATAATACCGTTGTGACGTTGTGTTTGTTTct                                                                                    >  1:305102/1‑250 (MQ=255)
                                                                                                                                                       tggcgtgcgagccgcTCTTTTCCGTCGGGCGGGGGGCCCTTGCGAGCTTGCCACGTCGTGTTAACCGAACGGGAACGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                     <  1:468942/237‑1 (MQ=255)
                                                                                                                                                           gggCGATGCGCTATTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGCCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                          >  2:252476/1‑224 (MQ=255)
                                                                                                                                                            ggCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTGGCGTTCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                              <  1:205839/251‑1 (MQ=255)
                                                                                                                                                              cGATGCGCTTGTTTTCTTGGGTGTGGGGGTCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCTCCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                            <  2:221312/251‑1 (MQ=255)
                                                                                                                                                                         ttttCTTCGGTCTGGGTGGACTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCCTG‑CTGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCCTTTATAAATCCCTGGAatt                                                                 >  2:71906/1‑251 (MQ=255)
                                                                                                                                                                            tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                               <  1:452552/250‑1 (MQ=255)
                                                                                                                                                                             cggcggtgggggTGGCCTTGCGATCGTGCCCGGGCGGGTTAACCGACCTGGAGGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattatt                                                              <  2:239983/241‑1 (MQ=255)
                                                                                                                                                                               tcGGTCTGGGTGGCCTGGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAtttt                                                            <  1:240834/250‑1 (MQ=255)
                                                                                                                                                                                  gTCTGGGTGGCCTTGAGATCGTGCCACGGCGTGTTAACCGCACTAAAAGCCCGAGCCGCACCGACGTGCAGGA‑C‑CAGGATCCGGTTTTCCTACGCCAGCTCCTTGAGGTCGCGGCCACTAGCCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCATCGGACCTGACCTAAACCTGCTGCTTAACCAGTTCGATAACACCGTTCTGACGTTGTGTTTGTTTCATTTATAAATCGCTGGAATTAttttcgtt                                                        >  1:254753/1‑251 (MQ=255)
                                                                                                                                                                                    cggggTGGCCGTGCGATGGTGCCACGGCGTGGTAACCGACCTGGATGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                       <  1:315344/248‑1 (MQ=255)
                                                                                                                                                                                           gCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                               >  2:402679/1‑251 (MQ=255)
                                                                                                                                                                                            ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                              <  1:228067/251‑1 (MQ=255)
                                                                                                                                                                                                      gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                     >  2:261560/1‑250 (MQ=255)
                                                                                                                                                                                                             ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                               >  1:110437/1‑249 (MQ=255)
                                                                                                                                                                                                              gCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAAGTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGc                            <  1:319369/251‑1 (MQ=255)
                                                                                                                                                                                                                  gttacccgacctgtacgccAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTGTTGCTCCGCCAGCTCATGGAGGTCGCGGCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTACCAGCGCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATGCCGTTTTGACTTTGTGTTTTTTTCAGTTATAAATCCCGGTAATTATTTTCGTTGTCGCGCATTGATCTAATCAACAAAAGCCATc                        <  2:438090/233‑1 (MQ=255)
                                                                                                                                                                                                                     acccaaacgggacggcgccgcccagccGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                     <  2:120929/226‑1 (MQ=255)
                                                                                                                                                                                                                      aCCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                     >  2:150641/1‑250 (MQ=255)
                                                                                                                                                                                                                              ggaatgcggccccggaCCGCCAGG‑ATGG‑CGCAGGAACAGTTTGTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGGCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAACCCTGCGGTTTAACGAGTTCGATAATACCGGTGTGACGTTGTGTTGGTTTCATGTATAAGTCCCTGGGATTATATTCGTTTTCGCGGATTAAGGGAAGCAACGAAGGCCAGCGCGAACCCCACg            <  2:233128/236‑1 (MQ=255)
                                                                                                                                                                                                                              ggAAGGCAGCCCCGCCCCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGCACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc             <  1:308541/250‑1 (MQ=255)
                                                                                                                                                                                                                               gAAGGCAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGATCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACgg           >  1:80132/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gccggcccgccccgccATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGGCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAACCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTaa       <  2:292615/245‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGCGCCGCACCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAc      <  1:246094/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gccgcTCCGCCATG‑ATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGCCCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg  >  2:457297/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       
CTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCA‑TTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATG‑ATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560207‑3560686

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