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.3% 191.2 / NA 58intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (25/32);  minor base T (1/0);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAT  >  NC_000913/3560212‑3560690
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAt                                                                                                                                                                                                                                       >  1:58170/1‑251 (MQ=255)
 ttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgaaggatc                                                                                                                                                                                                                                      >  2:162002/1‑243 (MQ=255)
     cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                                                                                                                                                                                                  >  1:344759/1‑251 (MQ=255)
               ttAGCAACGATGAGATTGTAGGTGAAAAAGAGCAAATTGCTTTGGTTAAGCAGTAAGTGTGCTAACGCTTCCGGCGTGGTGCGGGTCTCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTTGGCATTGCGCTCGTGCAACGGCGTGTTAACCGAACTGGAGGGCAGCGCGGCACCGCCATGATGGCGCAGGATAAGGATTTGCTctgc                                                                                                                                                                                                                        >  2:266971/1‑248 (MQ=255)
                       gTTGAGATTGTTCGTGACAATGCGCAGATTGCTGTGATTGAGCAGAGCGTGCGCTCACGCTTCCGGCGTGGTGCCGAAATCGATAAAAAGCGTCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTCGGTGGCCTTGCGATCGTGCCACGGCGTGTTCACCGACCTGGAAGGCAGCGCCGCACAGCCATGATGGCGCAGGCTCAGGGTTTGCTCCGcccacacct                                                                                                                                                                                                                >  2:212864/1‑244 (MQ=255)
                                   ggTGAAAAGGCGAAAATTGCTGGGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                     <  1:103045/250‑1 (MQ=255)
                                        cAATGCGTAATTTGCTGTGTTGGAGCAGTGCGTGCGCTCCCGCTCCCGGGTGGGTGCCGATAGCGAAAAACAGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                                                                                                               <  1:158313/251‑1 (MQ=255)
                                        cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTCCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAACCGGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                                                                                                               <  1:315454/251‑1 (MQ=255)
                                            gcgcAAATTGTTGTGATTGACCGGTGCGTGCGCTACCGCTCCGGGCGTGGTGCCGATATCGATAAACAGCGTCGAGTCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGCCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCt                                                                                                                                                                                           <  1:312280/251‑1 (MQ=255)
                                                 aaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTACGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                      >  2:363922/1‑251 (MQ=255)
                                                   ttGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                    >  1:186153/1‑251 (MQ=255)
                                                     gcggtgATTGAGAGGTGCGGGCGTTCCCGCTTCCGGCGTGGGGCCGATACCGATAAGCAGCGTCGAGCCATTGGGGATTGGCCCCGCCACTTTGCGGGCGATTCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGGGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa                                                                                                                                                                                  <  1:20425/248‑1 (MQ=255)
                                                         tgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAACCAGCGTCGAGCCAATAGGGATTTGCACTGCCACTTTCCGGGCGCCGCCCTCCTTTTATTCGCTCTGGCTGGCCTTGCGACCGTGCAACGCCGTGTTAACCCAACTGGAAGGCAACGCCGCACCCCACCGATGGGGCAGGTTAAGAATTTGCTCCGACAGCTCATTTCGACCGCGGCGCATAGTCTGCgggatcgcggt                                                                                                                                                                              >  1:196247/1‑243 (MQ=255)
                                                         tgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                               >  1:298755/1‑250 (MQ=255)
                                                           aTTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGCTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTTAGGTCGCGGCGACTAGTCTGCGGGCTGACgg                                                                                                                                                                               >  2:116466/1‑248 (MQ=255)
                                                             tGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCTCAGGATCAGGTTTTGCTCCGCCAGCCCATTGGGGGCGCGGCGGATAGTCTGCTGGCCGACGGAGaa                                                                                                                                                                           >  2:15055/1‑250 (MQ=255)
                                                                   gtgcgtgcGCTACCGCTCCGGGCGTGGTGCCGATTTCGATAAACAGCGTCGCGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctct                                                                                                                                                                    <  1:180498/251‑1 (MQ=255)
                                                                         gcgctcccctttCCGGCGGGGTCCAGATATCGATAACCAGTGCCGACCCTTGGGGGTTTTGCCCCGCCTCTTTGCGGGCTATGCGCTCTTTTTCTTCGGTCTGGTTGGCCTGGCAATCGTGCCACGGCGTGTAAACCGAACTGGAATGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                              <  2:359311/241‑1 (MQ=255)
                                                                                    tCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                                                                                                                   >  1:314445/1‑251 (MQ=255)
                                                                                            ggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCGTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                                                                                                                           <  1:55632/251‑1 (MQ=255)
                                                                                             gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                          >  1:140121/1‑251 (MQ=255)
                                                                                                 ccGATATCGATAAACCGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGTGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACActg                                                                                                                                       >  2:8660/1‑251 (MQ=255)
                                                                                                      atCGATAAACGGCGCCGAGCCGTTGGGTATTTCCTCCGCCACTTTGCGGGCGATGCGCTTTTTTTCTTCTTTCTGGGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                                  <  1:305947/251‑1 (MQ=255)
                                                                                                       tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTATGCGGGCTGACGGAGAAATGATCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                 >  2:70164/1‑251 (MQ=255)
                                                                                                            aaccaccgcggggcgatggggtatgtGCTCCGCCACTTGGGGGGCGGTGCGCTCTTTTTCTGCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGACCTGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                                                                                                                            <  2:31000/226‑1 (MQ=255)
                                                                                                                     ggAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                                                                                                                   <  1:54780/250‑1 (MQ=255)
                                                                                                                        gCCATTGGGGATTTGCCCCGCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                <  1:89612/251‑1 (MQ=255)
                                                                                                                          cATTGGGGATTTGCTCCGCCACTTTGCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                              >  1:353404/1‑251 (MQ=255)
                                                                                                                             tGGGGATTTGCTCCGCCTCTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                                                                                                            <  1:77978/250‑1 (MQ=255)
                                                                                                                                             cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                           <  1:96752/251‑1 (MQ=255)
                                                                                                                                                tgtGCGGGCGACGCGCCTTGTTTTTCCGGTCTGGGGGGCCTTGCGTTCGGGCCTCGGCGTGTTAACCGAACTGGAAGGCAGGCCCGCACCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGGACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTAGtgtttgttac                                                                                         <  2:417581/248‑3 (MQ=255)
                                                                                                                                                 gtGCGGGCGATGCTCCCTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTGACATAACCCGGCTGTTTAACCAGTTCGATAATACCTTTGTGACTTTTTGTTTGTTTc                                                                                         <  2:59172/248‑1 (MQ=255)
                                                                                                                                                   gggggggatacgcTCTTTTTCTTCGGGCTGTGTGGCCTTGCGAGCGTGCCCCGGCGTGTAAACCGAACGGGAGGGCAGCCCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCTGGCTGACGGAGAAATTCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGTCGTTGTGTTTGTTTCAttt                                                                                     <  2:416730/241‑1 (MQ=255)
                                                                                                                                                   gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                     <  1:239396/251‑1 (MQ=255)
                                                                                                                                                    cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                     >  1:129986/1‑250 (MQ=255)
                                                                                                                                                           tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                                                                                                                   >  2:322514/1‑181 (MQ=255)
                                                                                                                                                            gcgcTCTTTTTCTTCGGGCGGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGAGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTGGTTTCATTTATAAATcc                                                                             <  1:77214/250‑1 (MQ=255)
                                                                                                                                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                           >  1:330377/1‑251 (MQ=255)
                                                                                                                                                                   ttttCTTCGGTTGGGGTGGCCTGGCGAGCGCGCCCCGGGGTGTTAACCGAACGGGAAGGCTGCGCCGCACCTCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTTCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAat                                                                      <  2:286038/250‑1 (MQ=255)
                                                                                                                                                                     tttCTCCGGTCTGGGGGGCCTGGCGATCGGGCCGCGGCGTGTTAACCGAAGTGGAAGGCAGCGCCTCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGAGAAATGCTCTACCAGCTCTTCGGTACTGACATTACCCTTCTGTTTAACCAGTGCGATAATGCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                     <  1:262014/250‑1 (MQ=255)
                                                                                                                                                                      ttcttcGGTATGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACAAGCTCTTCGGTACTGACATAAACCTGCTGTTTAACCAGTTCGATAATAACGTTGTGACGTTGTGTTTGTTTCATTTATAAATACCTGGAattat                                                                   >  2:207696/1‑251 (MQ=255)
                                                                                                                                                                         ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGACAAATGCTTTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGGCGTTGTGTTTGTTTCATTTATAATTCCCTGGAATTAtttt                                                                >  1:60418/1‑251 (MQ=255)
                                                                                                                                                                           cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                              <  1:143311/251‑1 (MQ=255)
                                                                                                                                                                           cGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttc                                                               <  1:287953/250‑1 (MQ=255)
                                                                                                                                                                                ggggggccgtggcgagcgggCCAGGGCGTGTAAACCGAGCTGGAGGGCACCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGGTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                          <  2:232907/231‑1 (MQ=255)
                                                                                                                                                                                ggggTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttc                                                         <  1:162002/250‑1 (MQ=255)
                                                                                                                                                                                       ccTTGCGATCGTGCCACGGCGTGTTAACCTAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                  <  1:87616/251‑1 (MQ=255)
                                                                                                                                                                                       ccTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                  >  2:262015/1‑251 (MQ=255)
                                                                                                                                                                                              aTCGGGCCACGGGGGGTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGaa                                           <  2:187978/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                   <  2:129986/250‑1 (MQ=255)
                                                                                                                                                                                                       cggtgtgtTAACCGACCTGGAAGGCAGCCCCGCACCGCCATGATGGCGCATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                                  <  1:358818/247‑1 (MQ=255)
                                                                                                                                                                                                             gtTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAGTTATAAAGCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                             <  2:58170/250‑1 (MQ=255)
                                                                                                                                                                                                                    gAACTGGAAGGCAGCGACGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTAGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAc                     >  1:163994/1‑251 (MQ=255)
                                                                                                                                                                                                                      aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc                   >  1:175796/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggAAGGCACCGCCCCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCa                  <  1:261999/249‑1 (MQ=255)
                                                                                                                                                                                                                             ggCAGCGCCCCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa            <  2:298755/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGCATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg      >  1:214655/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCg    >  2:350304/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAt  >  2:276337/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
TTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAT  >  NC_000913/3560212‑3560690

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