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.2% 196.4 / ‑4.5 54intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (28/25);  total (29/25)
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.83e-01

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560206‑3560686
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
tcttctTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑c‑c                                                                                                                                                                                                                                        >  2:75905/1‑251 (MQ=255)
   ttttttCCCATCAGCGTTTTAGCAACGTTAAAATTGTTGGTGAAAATGCGAAATTTGTTGTGATGGAGCAGGGCTTGCGCTCCCG‑CTTCCGGCGTGGTGCC‑GATAC‑CGAAAAACAGCGTGGGGCCATTGGG‑GTTTTGC‑TCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cgca                                                                                                                                                                                                                                       <  1:246323/248‑4 (MQ=255)
   tcttTTCCCATCAGCGTGTTGGCATCGTTGAGATTGTGGGGGACAGTGCGCAAATTGCTGGGATTGAGCAGTGCGGGCGTTCCCG‑CTCCCGGGGGGGGGCC‑GATAT‑GGAAAAACGGCGTCGAGCCTTTGGG‑GTTTTGC‑CCGGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCGGCACCGCCATGATGG‑cg                                                                                                                                                                                                                                         <  1:148365/248‑2 (MQ=255)
       ttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                  >  1:420120/1‑251 (MQ=255)
        acccaaatttgtgttgcaaaCGTTGAGATTGTTGGTGACAATGCATAAATTGTTTAGATTGAGCAGGGCGTGCGCTACC‑‑CTTTCGGGGGTGTGCCATTTAT‑CTAAAAACAGCGCCG‑GCCCTTGGGTGGTTTGT‑CCCGCCACTTTGCGGGCGATGCGCCCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTACCCGAACTGGAATGCAGCGCCGCACCGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                  <  2:92545/232‑1 (MQ=255)
         cccATCAACGTGTGAGCAATGTTAAGATTTTTGTTCAAAATGCGCAATTTGCTGGAATTGAGCGGTGCGTGCCCACCCG‑CTCCCGGCGTGGTGCT‑GTTAC‑CAATAACCGGCTCCGACCCATGGGG‑TTTTTGC‑TCCGCCATTTTGCGGGCGATGCGCTCTTTTTCTTCGTTCTGGTTGGCCTTGCGATCGTGCCACGGCTTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAg                                                                                                                                                                                                                                <  1:89841/250‑1 (MQ=255)
              cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGAGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGAAAAACAGCGACGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGAGGGCGATGCGCTATTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGACGGCAGAGCCGCACCGCCATGATGG‑CGCAGGATAAGGgttt                                                                                                                                                                                                                           >  1:57183/1‑251 (MQ=255)
                   tgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                      >  1:235528/1‑251 (MQ=255)
                     taacaaacgttggttttgttgtggACAATGCGCAATTTGCTTTGATGGAGCGTTGTGTGCG‑TCCCGGCTTCCGGCGTGTTCCA‑GATAT‑CAATAAACGGCGTCAGGCCTTTGGG‑AATTTCC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCg                                                                                                                                                                                                                     <  1:352594/227‑1 (MQ=255)
                        gCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                   >  1:65180/1‑249 (MQ=255)
                          aaggTTGAGATTGTGGGTGAAAATGCGAAATTTGGGGTGATTGAGCAGTGGGGGGCCTCCCG‑CTTCGGGCGGGGTGCC‑GATTT‑CGAAAAACAGCGCGGAGCCTTTGGG‑GATTTGC‑CCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                               <  1:20260/248‑1 (MQ=255)
                          aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGAG‑CAGGATCCGGGTTTGCTCCGCCAgac                                                                                                                                                                                                               >  2:482165/1‑249 (MQ=255)
                          aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                               >  1:212909/1‑251 (MQ=255)
                            ggTTGAGATTTTTGGTGAAAATGCGAAATTTGTTGTGTTTGAGCAGTGCGGGCCCTCCCC‑TTTCCGGGGTGGTGCC‑GATAT‑CGATAAACAGCGTGGACCCTTTGGG‑GTTTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                                                                                                             <  1:120302/250‑1 (MQ=255)
                                    ttgttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                     >  1:51393/1‑251 (MQ=255)
                                            gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                                                                                                              >  1:159503/1‑250 (MQ=255)
                                            gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                                                                                                              >  1:417455/1‑250 (MQ=255)
                                                    gaaatttgtggtgATTGGGCAGGGCGTGCGGTACCG‑CTTCGGGCGGGGTGCC‑GATAC‑CGATAAAAAGCGTCGAGCAATTGGG‑GATTTGC‑TCCGCCATTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                                                                                                                     <  1:40656/241‑1 (MQ=255)
                                                      aaaTTGCTGTGATTGAGCAGTGCGTGCGCTATCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCCTTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                   >  2:128025/1‑251 (MQ=255)
                                                             tgtgATGGAGCAGGGGGTGCGCTCCCT‑TTTCCGGGGGGGGGCC‑GATTT‑CGATAACCAGCGTCG‑GCCAATGGG‑GGTTTGCCTCCGCCACTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTG‑TTAACCGAACTGGAGGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCTCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                               <  1:376145/248‑1 (MQ=255)
                                                                  attgggcgttggcggcgtcaccctttCCGGGGTGGTGCC‑AATAC‑CGAAAAACGCCGTCGAGCCATTGGG‑GTTTTTC‑TCCGCCACTTTGCGGGCGATGCCCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                         <  1:226803/226‑1 (MQ=255)
                                                                    gcgcatggcggcggctcccG‑CTTCCGGCGTGGGGTC‑GA‑ATGCAATAATCAGCGTCTAGCCTTTGGG‑GATTTGC‑CCCGCCCTTTTGGTGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTT‑TTAACCGAACTGGAAGGCAGCGCTGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                                                                                                                      <  1:457722/233‑1 (MQ=255)
                                                                       cagttctggcgcTCCCG‑CTTCCGGCGGGGGCCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑CCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                  <  1:330353/243‑1 (MQ=255)
                                                                            cgtgcgCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCAGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                             >  1:443935/1‑251 (MQ=255)
                                                                                   tACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                      >  1:222099/1‑251 (MQ=255)
                                                                                    acccgctgctgggtgggtccc‑GATAT‑TGATAACCAGCGTCGAGCCATGGGG‑GATTTGC‑TCCCCCACTTTGCGGGCGAGGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCTTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                      <  1:60412/232‑1 (MQ=255)
                                                                                           tCCGGCGGGGTGCC‑GATAT‑CGATAAACAGCGTCGGGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                                                                                                 <  1:75905/249‑1 (MQ=255)
                                                                                              ggCGTGGTGCC‑GATAT‑AGATAAACAGCGTCGAGCCATTGGG‑AATTTGC‑TCCACCACTTTGCGGGCGACGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAAGTGGAAGGCAGCGCCGCACCGCCCTGATGG‑CGCAGGAACAGGTTTTGCTCCGGCAGCTCATTTAGGTCGAGGCGAATAGACTGCGGGCTGACGTCGAAATGCTCTACCAGCTCTTCGGTAATGAc                                                                                                                                            >  2:286730/1‑251 (MQ=255)
                                                                                                              aacGATAAACAGCGTAGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCCGGTTTTGCTCCGCCACCGCATTGCGGTCACTGCGAATAGTCCGCGGGCTGACGGACAACGGATCTACCAGCACTTCGGCACTGCCATACACCTGCTGtt                                                                                                                              >  2:32506/3‑251 (MQ=255)
                                                                                                                      aCAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                      >  1:57834/1‑251 (MQ=255)
                                                                                                                              gggagcaggtgtgtttttcttCCGCCACTTGGGGGGCGAGGCGTTTTTTTTCTTCGGTGTGGGTGCCCTGGCGACGGTGCCACGGGGTG‑TTAACCGACCTGGAGGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGGTGAACCAGTTCGata                                                                                                                <  2:317986/231‑1 (MQ=255)
                                                                                                                                   aTTGGG‑GATTTGC‑TCCGCGACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                         >  2:403101/1‑251 (MQ=255)
                                                                                                                                      ggg‑gGTTTGCCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTGGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                       <  1:110635/251‑1 (MQ=255)
                                                                                                                                             ttGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                 <  1:380951/250‑1 (MQ=255)
                                                                                                                                                   cccccACTTTGCGGGCGTGGCGCTTTTTTTCTTGGGTCGGGGGGGCCTTGCGACGGTGCCACGGCGTG‑TTACCCGAACGGGAGGCCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGTGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                           <  2:417455/248‑1 (MQ=255)
                                                                                                                                                         cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                     >  1:29697/1‑251 (MQ=255)
                                                                                                                                                          tttGCGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGGGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                     <  1:54907/250‑1 (MQ=255)
                                                                                                                                                            tGCGGGCGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGACCCGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAt                                                                                   <  1:422752/250‑1 (MQ=255)
                                                                                                                                                               gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                               >  1:217023/1‑251 (MQ=255)
                                                                                                                                                               gggCGATGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAATGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                               >  1:334162/1‑251 (MQ=255)
                                                                                                                                                                       cgcTCTTTTTATTCGGTCTGGTGGGCGTTGGGAGTGTGCCCCGG‑GTGGTTAACCGAACTGGGGGGCAGCCCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                       <  2:193186/251‑1 (MQ=255)
                                                                                                                                                                                     ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGATTTAttttcgt                                                         >  2:159592/1‑251 (MQ=255)
                                                                                                                                                                                              ggCCTTGCGCTGGGGCCAGGGCGTG‑TCAGTCGCACTGGAAGGCTTCGCCGCCCGGCCATGATGG‑CGCAGGTTCAGTTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATTATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATGTTCGTTTTAGCGCa                                                <  2:153739/251‑1 (MQ=255)
                                                                                                                                                                                                ccTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                              >  1:294732/1‑251 (MQ=255)
                                                                                                                                                                                                         cGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                     <  1:144120/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                                >  1:76628/1‑184 (MQ=255)
                                                                                                                                                                                                                       gttaccggaccTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                         <  1:2193/240‑1 (MQ=255)
                                                                                                                                                                                                                        ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                        <  1:368178/250‑1 (MQ=255)
                                                                                                                                                                                                                            ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                               >  1:465200/1‑191 (MQ=255)
                                                                                                                                                                                                                                 cTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGATTTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc               >  2:125539/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggcgcgcccgcACCGCCATGATGG‑CGCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCGCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACACAAGCCATCGCTAAACCCACGGCTa        <  2:39590/244‑1 (MQ=255)
                                                                                                                                                                                                                                        gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                              >  2:345157/1‑180 (MQ=255)
                                                                                                                                                                                                                                        gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                              <  1:345157/180‑1 (MQ=255)
                                                                                                                                                                                                                                        gCAGCGCCGAACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGACGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATCACCCTGCTGTTTAACCAGCTCGATAATACCGGTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTACGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa        >  2:280061/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aGCGCCGCAACGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTCTTTTCTTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc     >  2:105829/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gccgcACCGCCATGCTGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg  >  1:394085/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGa                                                                    >  1:44779/1‑181 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCG‑CTTCCGGCGTGGTGCC‑GATAT‑CGATAAACAGCGTCGAGCCATTGGG‑GATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTG‑TTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560206‑3560686

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