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

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560206‑3560686
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
tcttctTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcc                                                                                                                                                                                                                                        >  2:75905/1‑251 (MQ=255)
   ttttttCCCATCAGCGTTTTAGCAACGTTAAAATTGTTGGTGAAAATGCGAAATTTGTTGTGATGGAGCAGGGCTTGCGCTCCCGCTTCCGGCGTGGTGCCGATACCGAAAAACAGCGTGGGGCCATTGGGGTTTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgca                                                                                                                                                                                                                                       <  1:246323/248‑4 (MQ=255)
   tcttTTCCCATCAGCGTGTTGGCATCGTTGAGATTGTGGGGGACAGTGCGCAAATTGCTGGGATTGAGCAGTGCGGGCGTTCCCGCTCCCGGGGGGGGGCCGATATGGAAAAACGGCGTCGAGCCTTTGGGGTTTTGCCCGGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCGGCACCGCCATGATGGcg                                                                                                                                                                                                                                         <  1:148365/248‑2 (MQ=255)
       ttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                  >  1:420120/1‑251 (MQ=255)
        acccaaatttgtgttgcaaaCGTTGAGATTGTTGGTGACAATGCATAAATTGTTTAGATTGAGCAGGGCGTGCGCTACCCTTTCGGGGGTGTGCCATTTATCTAAAAACAGCGCCGGCCCTTGGGTGGTTTGTCCCGCCACTTTGCGGGCGATGCGCCCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAATGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                  <  2:92545/232‑1 (MQ=255)
         cccATCAACGTGTGAGCAATGTTAAGATTTTTGTTCAAAATGCGCAATTTGCTGGAATTGAGCGGTGCGTGCCCACCCGCTCCCGGCGTGGTGCTGTTACCAATAACCGGCTCCGACCCATGGGGTTTTTGCTCCGCCATTTTGCGGGCGATGCGCTCTTTTTCTTCGTTCTGGTTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAg                                                                                                                                                                                                                                <  1:89841/250‑1 (MQ=255)
              cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGAGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGAAAAACAGCGACGAGCCATTGGGGATTTGCTCCGCCACTTTGAGGGCGATGCGCTATTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGAGCCGCACCGCCATGATGGCGCAGGATAAGGgttt                                                                                                                                                                                                                           >  1:57183/1‑251 (MQ=255)
                   tgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                      >  1:235528/1‑251 (MQ=255)
                     taacaaacgttggttttgttgtggACAATGCGCAATTTGCTTTGATGGAGCGTTGTGTGCGTCCCGGCTTCCGGCGTGTTCCAGATATCAATAAACGGCGTCAGGCCTTTGGGAATTTCCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                                                                                                                     <  1:352594/227‑1 (MQ=255)
                        gCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                   >  1:65180/1‑249 (MQ=255)
                          aaggTTGAGATTGTGGGTGAAAATGCGAAATTTGGGGTGATTGAGCAGTGGGGGGCCTCCCGCTTCGGGCGGGGTGCCGATTTCGAAAAACAGCGCGGAGCCTTTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                               <  1:20260/248‑1 (MQ=255)
                          aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGAGCAGGATCCGGGTTTGCTCCGCCAgac                                                                                                                                                                                                               >  2:482165/1‑249 (MQ=255)
                          aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                               >  1:212909/1‑251 (MQ=255)
                            ggTTGAGATTTTTGGTGAAAATGCGAAATTTGTTGTGTTTGAGCAGTGCGGGCCCTCCCCTTTCCGGGGTGGTGCCGATATCGATAAACAGCGTGGACCCTTTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                                                                                                             <  1:120302/250‑1 (MQ=255)
                                    ttgttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                     >  1:51393/1‑251 (MQ=255)
                                            gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                                                                                                              >  1:159503/1‑250 (MQ=255)
                                            gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                                                                                                              >  1:417455/1‑250 (MQ=255)
                                                    gaaatttgtggtgATTGGGCAGGGCGTGCGGTACCGCTTCGGGCGGGGTGCCGATACCGATAAAAAGCGTCGAGCAATTGGGGATTTGCTCCGCCATTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                                                                                                                     <  1:40656/241‑1 (MQ=255)
                                                      aaaTTGCTGTGATTGAGCAGTGCGTGCGCTATCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCCTTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                   >  2:128025/1‑251 (MQ=255)
                                                             tgtgATGGAGCAGGGGGTGCGCTCCCTTTTCCGGGGGGGGGCCGATTTCGATAACCAGCGTCGGCCAATGGGGGTTTGCCTCCGCCACTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCTCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                               <  1:376145/248‑1 (MQ=255)
                                                                  attgggcgttggcggcgtcaccctttCCGGGGTGGTGCCAATACCGAAAAACGCCGTCGAGCCATTGGGGTTTTTCTCCGCCACTTTGCGGGCGATGCCCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                         <  1:226803/226‑1 (MQ=255)
                                                                    gcgcatggcggcggctcccGCTTCCGGCGTGGGGTCGAATGCAATAATCAGCGTCTAGCCTTTGGGGATTTGCCCCGCCCTTTTGGTGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTTTTAACCGAACTGGAAGGCAGCGCTGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                                                                                                                      <  1:457722/233‑1 (MQ=255)
                                                                       cagttctggcgcTCCCGCTTCCGGCGGGGGCCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                  <  1:330353/243‑1 (MQ=255)
                                                                            cgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCAGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                             >  1:443935/1‑251 (MQ=255)
                                                                                   tACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                      >  1:222099/1‑251 (MQ=255)
                                                                                    acccgctgctgggtgggtcccGATATTGATAACCAGCGTCGAGCCATGGGGGATTTGCTCCCCCACTTTGCGGGCGAGGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                      <  1:60412/232‑1 (MQ=255)
                                                                                           tCCGGCGGGGTGCCGATATCGATAAACAGCGTCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                                                                                                 <  1:75905/249‑1 (MQ=255)
                                                                                              ggCGTGGTGCCGATATAGATAAACAGCGTCGAGCCATTGGGAATTTGCTCCACCACTTTGCGGGCGACGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAGTGGAAGGCAGCGCCGCACCGCCCTGATGGCGCAGGAACAGGTTTTGCTCCGGCAGCTCATTTAGGTCGAGGCGAATAGACTGCGGGCTGACGTCGAAATGCTCTACCAGCTCTTCGGTAATGAc                                                                                                                                            >  2:286730/1‑251 (MQ=255)
                                                                                                              aacGATAAACAGCGTAGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCCGGTTTTGCTCCGCCACCGCATTGCGGTCACTGCGAATAGTCCGCGGGCTGACGGACAACGGATCTACCAGCACTTCGGCACTGCCATACACCTGCTGtt                                                                                                                              >  2:32506/3‑251 (MQ=255)
                                                                                                                      aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                      >  1:57834/1‑251 (MQ=255)
                                                                                                                              gggagcaggtgtgtttttcttCCGCCACTTGGGGGGCGAGGCGTTTTTTTTCTTCGGTGTGGGTGCCCTGGCGACGGTGCCACGGGGTGTTAACCGACCTGGAGGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGGTGAACCAGTTCGata                                                                                                                <  2:317986/231‑1 (MQ=255)
                                                                                                                                   aTTGGGGATTTGCTCCGCGACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                         >  2:403101/1‑251 (MQ=255)
                                                                                                                                      ggggGTTTGCCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                       <  1:110635/251‑1 (MQ=255)
                                                                                                                                             ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                 <  1:380951/250‑1 (MQ=255)
                                                                                                                                                   cccccACTTTGCGGGCGTGGCGCTTTTTTTCTTGGGTCGGGGGGGCCTTGCGACGGTGCCACGGCGTGTTACCCGAACGGGAGGCCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGTGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                           <  2:417455/248‑1 (MQ=255)
                                                                                                                                                         cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                     >  1:29697/1‑251 (MQ=255)
                                                                                                                                                          tttGCGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                     <  1:54907/250‑1 (MQ=255)
                                                                                                                                                            tGCGGGCGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCCGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAt                                                                                   <  1:422752/250‑1 (MQ=255)
                                                                                                                                                               gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                               >  1:217023/1‑251 (MQ=255)
                                                                                                                                                               gggCGATGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAATGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                               >  1:334162/1‑251 (MQ=255)
                                                                                                                                                                       cgcTCTTTTTATTCGGTCTGGTGGGCGTTGGGAGTGTGCCCCGGGTGGTTAACCGAACTGGGGGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                       <  2:193186/251‑1 (MQ=255)
                                                                                                                                                                                     ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGATTTAttttcgt                                                         >  2:159592/1‑251 (MQ=255)
                                                                                                                                                                                              ggCCTTGCGCTGGGGCCAGGGCGTGTCAGTCGCACTGGAAGGCTTCGCCGCCCGGCCATGATGGCGCAGGTTCAGTTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATTATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATGTTCGTTTTAGCGCa                                                <  2:153739/251‑1 (MQ=255)
                                                                                                                                                                                                ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                              >  1:294732/1‑251 (MQ=255)
                                                                                                                                                                                                         cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                     <  1:144120/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                                >  1:76628/1‑184 (MQ=255)
                                                                                                                                                                                                                       gttaccggaccTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                         <  1:2193/240‑1 (MQ=255)
                                                                                                                                                                                                                        ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                        <  1:368178/250‑1 (MQ=255)
                                                                                                                                                                                                                            ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                               >  1:465200/1‑191 (MQ=255)
                                                                                                                                                                                                                                 cTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGATTTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc               >  2:125539/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggcgcgcccgcACCGCCATGATGGCGCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCGCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACACAAGCCATCGCTAAACCCACGGCTa        <  2:39590/244‑1 (MQ=255)
                                                                                                                                                                                                                                        gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                              >  2:345157/1‑180 (MQ=255)
                                                                                                                                                                                                                                        gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                              <  1:345157/180‑1 (MQ=255)
                                                                                                                                                                                                                                        gCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGACGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATCACCCTGCTGTTTAACCAGCTCGATAATACCGGTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTACGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa        >  2:280061/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aGCGCCGCAACGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTCTTTTCTTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc     >  2:105829/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gccgcACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg  >  1:394085/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGa                                                                    >  1:44779/1‑181 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  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: