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.G97.8% 148.1 / ‑5.6 46intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (19/26);  minor base T (0/1);  total (19/27)
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.66e-01

TTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTC  >  NC_000913/3560212‑3560698
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGAt                                                                                                                                                                                                                                               >  1:234937/1‑251 (MQ=255)
    ccATCAACGTGTTAGAAACGTTGAGATTGTTGGTGACAATGCGCCAATTGCTGTGATTGAGCAGTGAGTGCGCTAACGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTACGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAgg                                                                                                                                                                                                                                           >  2:118472/1‑251 (MQ=255)
     cATCTACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGGGCGTGCGCTCCCGCTGCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCTTTGGGGAGTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGt                                                                                                                                                                                                                                          <  2:133176/251‑1 (MQ=255)
              gtTAGCAACGTTGAGATTGTTGGTGACAAGGCGCAAATTGCGGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATGAACAGCGGCGAGCCATTGGGGATTTGCTCCGGCACTTTGCGGGCGATGCGCTCTTTTGCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGAGCAGGTGTTGCTcc                                                                                                                                                                                                                                  <  1:34251/250‑1 (MQ=255)
                       ggggatagtgggggtgcctaggcgcAAATTGCTGGGATTAAGCATGGCGTGCGCTACCGCTGCGGGCTGGGGCCGGATTGCGATAAACAGCGTGGTGCCTTGGGGGTGTGTCTCCTCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                          <  2:257128/228‑1 (MQ=255)
                          gagaTTCTCGGTGACAATGCGCAACTTGATGTGATTGAGCAGTGCGCGCGCTACCGCTTCCGGCGTGGAGACGATATCGATACACAGAGTCGCGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGGACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTgg                                                                                                                                                                                                                     >  2:189232/1‑250 (MQ=255)
                                 ttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTAGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                                                                                                                              >  1:306372/1‑251 (MQ=255)
                                                aaaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCAGATATCGATAAACAGCGTCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGTCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTAAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                               >  2:17321/1‑251 (MQ=255)
                                                      cggttaTTGAGGAGTGCGTGCGTTACCGCTTCCGGGGTGGTGCCGATATTGAAAAACGGCGTCGAGCCTTTGGGTATTGTCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                                                                                                                         <  2:279209/246‑1 (MQ=255)
                                                       ggggtgggggcggtccgtgcgCGCCCGCTTCCGGCGGGGGGCCGATATCGATAAAGAGCGTCGAGCAATGGGGGATTGGTGCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                        <  2:232467/236‑1 (MQ=255)
                                                        ggtaTTGGGCAGTGCTGGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                        <  1:261002/247‑1 (MQ=255)
                                                                                  cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                              >  1:226215/1‑250 (MQ=255)
                                                                                       ggCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                                 >  2:57250/1‑194 (MQ=255)
                                                                                       ggCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                                 <  1:57250/194‑1 (MQ=255)
                                                                                        gCGTGGTGCCGATATCGATAAACAGCGTCGCGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGTCATGATGGCGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCTGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                                                                                                                                       <  1:53079/251‑1 (MQ=255)
                                                                                          gtggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGTAAGGCAGCGCCGCAGCGCCATGATTGCTCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGTTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                                     <  1:189232/251‑1 (MQ=255)
                                                                                            ggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                                                                                                                                   >  2:135831/1‑251 (MQ=255)
                                                                                              ggCCGATATCGATAAACAGCGGCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                                 <  1:124971/250‑1 (MQ=255)
                                                                                                ccGATATCGATAATCAGCGTCGAGCCATTGTGGATTTGCGCCGCCACTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGTGCGGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                                                                                                                               <  2:313738/251‑1 (MQ=255)
                                                                                                       ccGATAAACAGCGCTGAGCATTTGGTGATTTGCTCCGCCGCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGTGCTGACGGAGAAATGCTCTGCCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                         <  2:27049/250‑1 (MQ=255)
                                                                                                        caaTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGTTACTGACATAACCCTGCTGTTTa                                                                                                                                        <  1:239764/249‑1 (MQ=255)
                                                                                                               atgtgtcgcaccATTGGGGATTTGTTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTATCCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGATAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                                 <  2:230649/241‑1 (MQ=255)
                                                                                                                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGAGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTATACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCgg                                                                                                                              >  1:276224/1‑250 (MQ=255)
                                                                                                                       ggCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                                                                                                                         <  1:135831/250‑1 (MQ=255)
                                                                                                                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                                            <  1:182662/243‑1 (MQ=255)
                                                                                                                            ttGGGCATTTGCTCCGCCACTTTGCGGGCGATGAGCTATTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGTCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGATCGCGGCGAATAGTCTGCGGGCTGCCGGCGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTTCTGTTTAACCAGTTCGatactctcgg                                                                                                                     >  2:182662/1‑243 (MQ=255)
                                                                                                                               gtgATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCTCCATGATGGCGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGTTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg                                                                                                                 <  2:235426/249‑1 (MQ=255)
                                                                                                                                 gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGTCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                                               <  1:40325/251‑1 (MQ=255)
                                                                                                                                            gCCAATTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAATGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGTCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                     <  1:101989/251‑1 (MQ=255)
                                                                                                                                             ccACTTTGCGGTCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgt                                                                                                    <  1:306791/251‑1 (MQ=255)
                                                                                                                                                cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAACGTTGTGACGTTGTGTTTGTTTc                                                                                                 >  2:300573/1‑251 (MQ=255)
                                                                                                                                                         cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                         >  1:217702/1‑202 (MQ=255)
                                                                                                                                                         cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                         <  2:217702/202‑1 (MQ=255)
                                                                                                                                                           aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGTTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                                      <  1:204896/251‑1 (MQ=255)
                                                                                                                                                           aTGCGCTCTTTTTCTTCGGGCGGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                                      <  2:234937/251‑1 (MQ=255)
                                                                                                                                                               gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTg                                                                                  >  1:107523/1‑251 (MQ=255)
                                                                                                                                                                           cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttc                                                                       <  1:142230/250‑1 (MQ=255)
                                                                                                                                                                             gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtt                                                                    >  1:257350/1‑251 (MQ=255)
                                                                                                                                                                                             gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCTGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                  <  1:255560/141‑1 (MQ=255)
                                                                                                                                                                                             gATCGTCCCACGGCGTGTTACCCGCACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                  >  2:255560/1‑141 (MQ=255)
                                                                                                                                                                                                        ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                                       <  1:60941/189‑1 (MQ=255)
                                                                                                                                                                                                        ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                                       >  2:60941/1‑189 (MQ=255)
                                                                                                                                                                                                          cGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                       <  1:157807/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCAGAGCCGCAGCGCCACGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCACCAAAAGCCATCGATAAACCCACgtct                     >  2:281435/1‑248 (MQ=255)
                                                                                                                                                                                                                             ggCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt                     >  1:308829/1‑250 (MQ=255)
                                                                                                                                                                                                                                               tgatgGCGCAGGATCAGGTTTTGCTCCGACAGCTCATTGCGGTCGCGGCGAATAGCCTGCGGGCTGACGGAGACATGATCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGCGTGCTc  >  2:152080/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
TTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTC  >  NC_000913/3560212‑3560698

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