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.1% 179.8 / ‑4.2 53intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (25/27);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560206‑3560699
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
tcttctTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCAGCACCGCCATGATGG‑cg                                                                                                                                                                                                                                                       >  1:102889/1‑250 (MQ=255)
 cttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGCTGG‑cgc                                                                                                                                                                                                                                                      >  1:338675/1‑249 (MQ=255)
        tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTGACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                                >  1:167169/1‑250 (MQ=255)
         aCCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGCAGGCAGCGCCGCACCGCCATGATGG‑cgctggatc                                                                                                                                                                                                                                                >  2:335836/1‑241 (MQ=255)
          ccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAgg                                                                                                                                                                                                                                             >  1:277978/1‑251 (MQ=255)
            aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtt                                                                                                                                                                                                                                           >  2:362245/1‑251 (MQ=255)
              cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGGGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGAAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtttt                                                                                                                                                                                                                                         >  1:119910/1‑251 (MQ=255)
                     ttAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGAG‑CAGGATCAGGTTTTGCTCCGc                                                                                                                                                                                                                                  >  1:88049/1‑251 (MQ=255)
                       atcactgTTGAGATTGTTGGTGAAAGTGCGCAAATTGCTGTGTTTGTGGGGTGCGTGCGCTCCCGCTCCGGGCGTGGTGCCGATATCGAAAAACAGCGGCGAGCCATGGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                                 <  1:450960/244‑1 (MQ=255)
                         cAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                                                                                                                                                                                                              >  1:198950/1‑251 (MQ=255)
                           aCGTTGAGATTTTTGTTGAAAATGCGAAAATTGTTGTGGTGGAGCAGGGGGGGCGCTACCGTTTCCGGGGTGGTGCCGAAATCGAAAAACAGCGCCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                            <  1:107349/251‑1 (MQ=255)
                                             aCAATGCGCAAAATGCTGTGATTGAGCAGAGCGTGCGCTACCGCTTCCGGCGTGGTGACGATATCGATAAACAGCGTCGAGCCATGGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGGGCTCTTTTTCTTCGGTCGGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGTGCCGCACCGCCATGATGG‑CGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGACGCGGCGAAt                                                                                                                                                                                                          >  2:353617/1‑251 (MQ=255)
                                                 ggcgcAAATTGTTGTGATTGAGCAGTGCGGGCGCTCCCGCTTCCGGGGTGGTGCCGATACCGATAAACGGCGTCG‑GCCCTGGGGGAATTTG‑CCCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTGCGGTCGGGGGGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                                                                                                                                      <  1:154636/250‑1 (MQ=255)
                                                    gcAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGCCTGc                                                                                                                                                                                                   >  2:90636/1‑251 (MQ=255)
                                                             tgtgATTGAGCAGAGGGTGCGCTACCGCTTCCGGCGTGGAGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCCGGCTCAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCGGCGACTAGTCTGCGGGCTgc                                                                                                                                                                                            >  2:429656/1‑248 (MQ=255)
                                                                    gAGCAGGGGGTGCGCTACCGCTTCGGGGGGGGGGCGGATATCGATAAACGGCGCCGAGCCATTGGGGT‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                                   <  1:171848/251‑1 (MQ=255)
                                                                                cgcTCCCGCTTCCGGGGGGGTCCCGATATCGATAAACGGCGGGGGGCC‑TTGGGGT‑TTTGGCTCC‑GCCACTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                                                                                                                                       <  1:82281/251‑1 (MQ=255)
                                                                                       gcctCCGCCGTGGTGCCGATATCGATAAACAGCGCCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGATGG‑CGCAGGCTCCGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGCGTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCggg                                                                                                                                                                >  2:329222/4‑250 (MQ=255)
                                                                                       gCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                                >  1:76521/1‑251 (MQ=255)
                                                                                          tccggcttggtcccgcattCGATAACCAGCGTCG‑GGCATTTGGGT‑TTTG‑CTCCCGCCTCTTTGCGGGCGATGCGCTCTTTGTCTTCGTTCTGGGTGGCCTGGCGTTCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                                                              <  1:201942/232‑1 (MQ=255)
                                                                                          ccgggggtggtgCCGATCTCAAAAAACAGCGTCGGGCCATTGGGGT‑TTTG‑CTCC‑GCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                                                                                                                             <  1:328807/245‑1 (MQ=255)
                                                                                             ggCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                                                                                                                          >  1:451627/1‑251 (MQ=255)
                                                                                                         atatCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                                              >  2:300154/1‑251 (MQ=255)
                                                                                                           atCGATAAACAGCGTCGGGCCATTGGGGA‑TTTG‑CCCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                                            <  1:300154/251‑1 (MQ=255)
                                                                                                                           gagccctttggtatttgg‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTGCTCCGTTTTGGTTGGCCGTGCAATCTTGCCCCGGCGTGTTAACCGACCTGGAAGGCAGCTCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTATCCAGTTCGat                                                                                                                               <  1:230229/232‑1 (MQ=255)
                                                                                                                           gAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCACCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGttaa‑t                                                                                                                           >  1:116358/1‑251 (MQ=255)
                                                                                                                            ggccctggggttttttG‑CTCC‑GCCCCGTGGCGGGCGATTGGCTCTTTTTCTTCGGTTTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa‑t                                                                                                                           <  2:198950/238‑1 (MQ=255)
                                                                                                                                     ggtgtttgccccc‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCTTGTTACCCGACCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTtg                                                                                                                    <  1:165679/239‑1 (MQ=255)
                                                                                                                                              cTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTtgtg                                                                                                          >  2:377640/1‑251 (MQ=255)
                                                                                                                                                         gtGCGGGCGATGCGCTTTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTTCCACGTCGTGTTAACCGAACTGGAGGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCa                                                                                                 <  1:197115/249‑1 (MQ=255)
                                                                                                                                                            cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                              >  2:23563/1‑250 (MQ=255)
                                                                                                                                                             gggCGATGCGCTTTTTTTCGGGTGTCTGGGTGGCCTGGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAGGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                            <  2:114786/251‑1 (MQ=255)
                                                                                                                                                              ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                            >  1:95764/1‑250 (MQ=255)
                                                                                                                                                                 gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                                        >  1:290860/1‑251 (MQ=255)
                                                                                                                                                                          tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCCGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGGGTTTGTTTCATTTATAAATCCCTGGAat                                                                               >  1:422067/1‑251 (MQ=255)
                                                                                                                                                                              tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                            <  1:335836/250‑1 (MQ=255)
                                                                                                                                                                              tctgcgggcggggTGGCCTTGCGATCGTGCACCGGCGTGTTAACCTAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCTAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTCAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                            <  1:35589/240‑1 (MQ=255)
                                                                                                                                                                                    gTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGGACTGCCTTAACCCTGCTGGTGAACCAGTTAG‑TAACTACCGTTGTGACGTTGTGTTTGTTTTATTTATAAATTCCTGGAATTAttttcgtt                                                                     >  2:15799/1‑251 (MQ=255)
                                                                                                                                                                                       tGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGAATTAttttcgttttc                                                                  >  2:38414/1‑251 (MQ=255)
                                                                                                                                                                                        gggtgtcctggggagcgttccACGGGGTGTTAACCTAACGGGAAGCAAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCGACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGGTGTGTTTGTTGCATGGATAAATCCCTGGAATGAttttcgttttcg                                                                 <  2:22838/232‑1 (MQ=255)
                                                                                                                                                                                        gggTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                                    <  1:85000/248‑1 (MQ=255)
                                                                                                                                                                                           gggCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                              <  1:358367/250‑1 (MQ=255)
                                                                                                                                                                                           gggCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                              <  2:61781/250‑1 (MQ=255)
                                                                                                                                                                                                  gcggtcgggccccGGCGTGTTAACCGAACTGGAAGGCCGCGCCTCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                                       <  2:456798/239‑1 (MQ=255)
                                                                                                                                                                                                  gCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                                       <  1:362245/251‑1 (MQ=255)
                                                                                                                                                                                                   cGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCg                                                      <  1:324442/251‑1 (MQ=255)
                                                                                                                                                                                                               ggCGTGTTAACCGAACTGGAAGGCAGCCCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCATCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAg                                          <  1:382119/251‑1 (MQ=255)
                                                                                                                                                                                                                  gtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                        >  2:459158/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGATCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGAGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCAATCGCt                                  >  2:280188/1‑251 (MQ=255)
                                                                                                                                                                                                                         ccGAACTGGAGGGCAGCGCGGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGATCGAATCAACAAAAGCCATCGCTaa                                <  1:329222/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGc                       >  1:340192/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccgccgcCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAccc                 <  2:116358/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCg             >  1:362987/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTtgt                                                                                                                   >  2:445979/1‑146 (MQ=255)
                                                                                                                                                                                                                                               aCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTtgt                                                                                                                   <  1:445979/146‑1 (MQ=255)
                                                                                                                                                                                                                                                       gatgG‑CGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTcc  <  1:303366/251‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGA‑TTTG‑CTCC‑GCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAA‑TACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560206‑3560699

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