Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,493 G→A 16.2% pseudogene (111/149 nt) glpR ← pseudogene, DNA‑binding transcriptional repressor,regulator, Energy metabolism, carbon: Anaerobic respiration, repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4930GA16.2% 149.4 / 24.0 77pseudogene (111/149 nt)glpRpseudogene, DNA‑binding transcriptional repressor,regulator, Energy metabolism, carbon: Anaerobic respiration, repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base G (43/21);  new base A (8/4);  total (51/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

ATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAA  >  NC_000913/3560269‑3560727
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
ccggggcgcggggtgcgCTACCGCTTCCGGGGGGGGGGCGATTACGATAACCGGGCGGGGGGCCTGTGGGGTTTGGCCCCCCCACTGTGCGGGCGAGGCGCTTTTTGTCTTGGGTTTGGGGGGCCTTGGGACGGGGCCCGGGCGGGTTAACCAAACTGGAAGGCAGCGCCGCGCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                                                                                                                                                                                                                  <  2:617305/226‑1 (MQ=255)
 ttGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCCCATTGAGGCCGCGGGGAATAGTCGCggg                                                                                                                                                                                                                              >  2:492085/1‑239 (MQ=255)
              gcgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCGGCTCATTGAGGTCCCGGCGAATAGCCTGCGGGGTGACGGGGaaa                                                                                                                                                                                                                  >  1:749914/1‑239 (MQ=255)
              gcgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                                                        >  1:617305/1‑233 (MQ=255)
                 ctcccGCTTCGGGGTTGTGGCGGATTTCGATAAACGGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                                                                <  1:176552/235‑1 (MQ=255)
                   aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                                                             >  1:839789/1‑239 (MQ=255)
                       cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                                                         >  1:223444/1‑239 (MQ=255)
                          ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTATCCGACCTGGAAGGCAGCGCCGAACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCAGCGAATACTCTGCGGTCTGACGGGGAAATGCTccaccac                                                                                                                                                                                                       >  2:638399/1‑237 (MQ=255)
                          ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCCCCATGCAGGCGCAGGATCCCGCTTTGCTCCGCCACCCCATTTGGGTCGCGGCGCACAGACTGCGGGCGGCCGGGCAAATGcctaacccc                                                                                                                                                                                                       >  2:723918/1‑230 (MQ=255)
                          ccGGCGTGGTGCCGATATCGATAAACAGCGACGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCACCGCCCGATCGCCATGATGGCGCAGGTTCAGGCTTTACCCTCCCAGCCCATTCAGGTGCGGGGGAATGGCCTGGGGGCTGCCGGTGAAAAGCTTTAcac                                                                                                                                                                                                        >  2:32742/1‑235 (MQ=255)
                                tggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTTAGGTCGCGGCGAATAGTCTGCGGGCTGTCGGAGAAATGCTCTACCAGGTCTTCg                                                                                                                                                                                                >  1:267044/1‑239 (MQ=255)
                                  gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATTCGGt                                                                                                                                                                                              >  1:702402/1‑239 (MQ=255)
                                  gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTATGCGGGCTGACGGAGAAATGCTCTACCAACTCTTCGGt                                                                                                                                                                                              >  1:372742/1‑239 (MQ=255)
                                             gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAAGAGTCTGAGGGCGGACGGAAAAATGGGCTACCAACTCTTCGGTCCTTACaaaa                                                                                                                                                                                   >  2:23843/1‑238 (MQ=255)
                                             gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                                                                                                                                                                   >  1:187778/1‑238 (MQ=255)
                                                       cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGGACTGACATAACCCTGCTGtt                                                                                                                                                                         >  1:875439/1‑239 (MQ=255)
                                                       cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGTAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                                                                         >  1:144470/1‑239 (MQ=255)
                                                       cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCCCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACCTAACCCTGCTGtt                                                                                                                                                                         >  1:670087/1‑239 (MQ=255)
                                                            aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGGGAATAGCCTGCGGGCTGAGGGGGAAATGCCCTACCAGCCCTTCGCGACTGACTAAACCCTGCTGTTTAAcc                                                                                                                                                                    >  2:561793/1‑239 (MQ=255)
                                                            aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTCCGGGCTGCCGGAGAAAGGCTCTCCCAGCCCTCCGGTACTGAAAAAACCCCGCTGTTTAAca                                                                                                                                                                    >  2:495842/1‑238 (MQ=255)
                                                               tatgggggTTTGGCCCCGGCAGTTGGCGGGCGATGCGTTTTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                                                  <  1:561793/234‑1 (MQ=255)
                                                               ctttGGGTTTTTGCCCCGCCATTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                                                  <  1:573447/236‑1 (MQ=255)
                                                                 ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCCACTCTTCGGGACGGCCAAAACCCGCCCGTTTAACCCGTTc                                                                                                                                                               >  1:419222/1‑239 (MQ=255)
                                                                 ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGGGAAATGCTCTACCAGCTCTTCGGTACTGACCTAACCCTGCTGTTTAAc                                                                                                                                                                     >  1:525607/1‑233 (MQ=255)
                                                                 ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGACCACGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAACAGTCTGCGGGCGGGCGGCAAAATGCGCTACCATCTCTTCGGGACGTAAATAACCCTGTTGTTTAACCCGt                                                                                                                                                                 >  2:758610/1‑237 (MQ=255)
                                                                       ttttGCCCCCCCATTTGGGGGGCCGTCCCCTTTTTTGCTTGGGTCGGGGGGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAAGTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                                                                                                                                                          <  2:894732/237‑1 (MQ=255)
                                                                        tttcccccGCCCTTTTGGGGGGGATGCGTTCTTTTTCTTCGGTCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                                                                                                                                                        <  1:775097/233‑1 (MQ=255)
                                                                         ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCGGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCACAAACCCCGCTGTTTAACCAGttggtaaaaa                                                                                                                                                        >  2:158047/1‑230 (MQ=255)
                                                                           gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCACTTCGGTACTGACATAACCCTGCTGGTTAACCAGTTCGTTAATACCg                                                                                                                                                     >  1:387641/1‑239 (MQ=255)
                                                                                gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt                                                                                                                                                 >  1:761730/1‑238 (MQ=255)
                                                                                 ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGACCTGACATAACCCAGCTGTTTAACAAGTTCGAAAAAAACGtttggg                                                                                                                                               >  2:417764/1‑235 (MQ=255)
                                                                                            gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGggtttg                                                                                                                                    >  1:95985/1‑239 (MQ=255)
                                                                                               gtGCGCTCTTTTTCCTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGGTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                                                                  <  1:78239/237‑1 (MQ=255)
                                                                                                       tttttCTTGGGTCGGGGTCCCCTTGCGCTCGGGCCGCGGCTTTTTAACCGAGCTGGAGGGCAGCCCCGCACCGCCATGATGGCGCAAGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                                                         <  1:32742/239‑1 (MQ=255)
                                                                                                          ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGTTCACGGAGAAATGCTCTACCAGCTCTCCGGTACTAACAGAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGGCGGTGGGTTTGTTTCATTTATaaa                                                                                                                       >  2:203004/1‑238 (MQ=255)
                                                                                                                 gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTTTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGTTCGAAATTACCGTTTTGACGTTGTGTTTTTTTCATTTATAAATCCCTAGa                                                                                                               >  2:228747/1‑239 (MQ=255)
                                                                                                                       gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTTACCTAACCCTGCTGTTTAACCCGTTCGATAATACCTTTGTGCCGTTGTTTTTTTTTCTTTTTAAAAACCCTGGActtatt                                                                                                         >  2:595521/1‑239 (MQ=255)
                                                                                                                                cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTGTGTTTGTTTCATTTCTAAATCCCTGGAATTAttttggtttt                                                                                                 >  1:897123/1‑238 (MQ=255)
                                                                                                                                  agtgggCCCCGGCGTTTTAACCGAATGGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCCGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcattttcg                                                                                               <  1:492085/233‑1 (MQ=255)
                                                                                                                                     gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCa                                                                                           >  1:712787/1‑239 (MQ=255)
                                                                                                                                            ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGGGAAATGCTCTACCAGCTCTTCGGTACTGCCCTAACCCTGCTGTTTAACCAGTTCGGAAATACCGTTGTGACGTTGTGTTTGTTTCTTTTATAAATCCCTTGCATTTTTTTCGTTTTAGCGCCTTGAGCg                                                                                    >  2:605956/1‑239 (MQ=255)
                                                                                                                                               gtgtTAGCCGAGTGGGAGGGCGGCGCCCCGCCGCCGTGATGGCGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                                                                 <  2:702402/239‑1 (MQ=255)
                                                                                                                                                tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCGTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCACCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTTGTTCATTTATAAATCCCCTGAATTATTTTCGTTTTCTCGCaatatggcaaca                                                                                >  2:342309/1‑228 (MQ=255)
                                                                                                                                                        gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTTTAAATCCTTGGAATTATTTTCGTTTTCCCCCTTTGGGCGAATCAACaaaa                                                                         >  1:347825/1‑238 (MQ=255)
                                                                                                                                                          acgggAAGGCAGGGCCGCACCGCCATGATGGCGCAGGATTAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:723918/236‑1 (MQ=255)
                                                                                                                                                          acgggAAGGCAGCGCCGCACCGCCATGATGGTGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:495842/236‑1 (MQ=255)
                                                                                                                                                          accgggaGGCAGGGCCGCGACGCCATGATGGCGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCCCAGCTCTTCGGTACTAACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:203004/233‑1 (MQ=255)
                                                                                                                                                          aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTTTGACGTTTGGCTTGTTTCATTTATAAATGCCTGGAATTATTTTCGTTTTCCGTCATTTAGCGCACAAACaaaaac                                                                       >  2:35107/1‑236 (MQ=255)
                                                                                                                                                          aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGAACCGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTTTGTTTGTTTCCTTTATAAATCCCTCGCAATATTTTCGTTTTCTCGCATTTAGCAAAACAaaaaaaaca                                                                      >  2:812647/1‑236 (MQ=255)
                                                                                                                                                          aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTGACATAACCCTGCTGGTTAACCAGTTCGTTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTTGGCCATTGAGGGAATAAACAAAAGcc                                                                      >  2:250260/1‑239 (MQ=255)
                                                                                                                                                           tgggagggagccgcgcccccGCCATGATGGGGGAGGATGAGGGTTTGCTCCGCCCGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCa                                                                     <  2:95985/221‑1 (MQ=255)
                                                                                                                                                                aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattatt                                                                                                         >  1:189570/1‑198 (MQ=255)
                                                                                                                                                                 ggccccccccccccgcatttgtgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                                                                <  2:873336/217‑1 (MQ=255)
                                                                                                                                                                     cccccgccccgccATGTTGGCGGAGTTTCAGGTTTTGCTGCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcc                                                           <  2:875439/236‑1 (MQ=255)
                                                                                                                                                                      cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcca                                                          >  1:45437/1‑238 (MQ=255)
                                                                                                                                                                        ccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                                                    <  2:713154/179‑1 (MQ=255)
                                                                                                                                                                        ccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                                                    >  1:713154/1‑179 (MQ=255)
                                                                                                                                                                           cACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                                                       >  1:597298/1‑237 (MQ=255)
                                                                                                                                                                             cccccATGATGGCGCAGGTTCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa                                                   <  2:187778/236‑1 (MQ=255)
                                                                                                                                                                                    gatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGc                                            <  1:417764/239‑1 (MQ=255)
                                                                                                                                                                                     ggggcgcggggttgggTTTTGCCCCGCCAGTTCATTGGGGTCGCGGGGAATAGTCTGCGGGCTGAGCGAAAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCg                                           <  2:298190/225‑1 (MQ=255)
                                                                                                                                                                                      tgggtcAGGATTAGGTTTTGGCTCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATTACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGa                                          <  1:783059/234‑1 (MQ=255)
                                                                                                                                                                                         cgcAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGaggg                                       >  1:731310/1‑236 (MQ=255)
                                                                                                                                                                                           cAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGTtgtgt                                    >  1:298190/1‑239 (MQ=255)
                                                                                                                                                                                           cAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGCCTAACCCGGCGTtgtgt                                    >  1:811729/1‑239 (MQ=255)
                                                                                                                                                                                                   ggTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCCCTAAACCCACGGCTAACCCGGCCAGGGTTGCTCCGt                             >  2:899036/1‑238 (MQ=255)
                                                                                                                                                                                                        tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGACTCCACAAAACCCATCGCTAAACCCACGGCTAACCCGCCACGTGGTGCTCCGTTCGcac                       >  2:723354/1‑237 (MQ=255)
                                                                                                                                                                                                        tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGcc                        >  1:528488/1‑238 (MQ=255)
                                                                                                                                                                                                           tCCGCCTGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCg                    <  1:228747/239‑1 (MQ=255)
                                                                                                                                                                                                           tCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAGCCCTGCTGTTTAGCCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCg                    <  1:638399/239‑1 (MQ=255)
                                                                                                                                                                                                            ccGGCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGa                   <  1:158047/239‑1 (MQ=255)
                                                                                                                                                                                                              gCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGAATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGACGTGTGCTCCGTTCGCCATCgcc                  >  1:785332/1‑236 (MQ=255)
                                                                                                                                                                                                               ccAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAATCCATCGCTAAAACCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCTACa                 >  1:747595/1‑238 (MQ=255)
                                                                                                                                                                                                                cAGCTCATTGAGGTCGCCGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATc               <  1:672063/239‑1 (MQ=255)
                                                                                                                                                                                                                            gTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                                                                 >  1:8529/1‑82 (MQ=255)
                                                                                                                                                                                                                            gTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                                                                 <  2:8529/82‑1 (MQ=255)
                                                                                                                                                                                                                             tCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAACCCATCGCTAAACCCACGCCTACCCCGGCGATGTGTGCTCCCTTCCCCATCGACATCCCAAACAAATCaa  >  2:438258/1‑239 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
ATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAA  >  NC_000913/3560269‑3560727

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