Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,494 G→A 15.8% pseudogene (112/150 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,4940GA15.8% 142.9 / 23.1 73pseudogene (112/150 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/18);  new base A (8/3);  total (51/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

TTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCGACATCCCAAACAAATCAA  >  NC_000913/3560272‑3560728
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
ttGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCCCATTGAGGCCGCGGGGAATAGTC‑GCggg                                                                                                                                                                                                                              >  2:492085/1‑239 (MQ=255)
             gcgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCGGCTCATTGAGGTCCCGGCGAATAGCCTGCGGGGTGACGGGGaaa                                                                                                                                                                                                                  >  1:749914/1‑239 (MQ=255)
             gcgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                                                        >  1:617305/1‑233 (MQ=255)
                ctcccGCTTCGGGGTTGTGGCGGATTTCGATAAACGGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                                                                <  1:176552/235‑1 (MQ=255)
                  aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                                                             >  1:839789/1‑239 (MQ=255)
                      cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                                                         >  1:223444/1‑239 (MQ=255)
                         ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTATCCGACCTGGAAGGCAGCGCCGAACCGCCATGATGG‑CGCCGGATCAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCAGCGAATACTCTGCGGTCTGACGGGGAAATGCTccaccac                                                                                                                                                                                                       >  2:638399/1‑237 (MQ=255)
                         ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCCCCATGCAGG‑CGCAGGATCCCGCTTTGCTCCGCCACCCCATTTGGGTCGCGGCGCACAGACTGCGGGCGGCCGGGCAAATGcctaacccc                                                                                                                                                                                                       >  2:723918/1‑230 (MQ=255)
                         ccGGCGTGGTGCCGATATCGATAAACAGCGACGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCACCGCCCGATCGCCATGATGG‑CGCAGGTTCAGGCTTTACCCTCCCAGCCCATTCAGGTGCGGGGGAATGGCCTGGGGGCTGCCGGTGAAAAGCTTTAcac                                                                                                                                                                                                        >  2:32742/1‑235 (MQ=255)
                               tggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTTAGGTCGCGGCGAATAGTCTGCGGGCTGTCGGAGAAATGCTCTACCAGGTCTTCg                                                                                                                                                                                                >  1:267044/1‑239 (MQ=255)
                                 gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATTCGGt                                                                                                                                                                                              >  1:702402/1‑239 (MQ=255)
                                 gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTATGCGGGCTGACGGAGAAATGCTCTACCAACTCTTCGGt                                                                                                                                                                                              >  1:372742/1‑239 (MQ=255)
                                            gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCCGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAAGAGTCTGAGGGCGGACGGAAAAATGGGCTACCAACTCTTCGGTCCTT‑ACaaaa                                                                                                                                                                                   >  2:23843/1‑238 (MQ=255)
                                            gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATaa                                                                                                                                                                                   >  1:187778/1‑238 (MQ=255)
                                                     cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGGACTG‑ACATAACCCTGCTGtt                                                                                                                                                                         >  1:875439/1‑239 (MQ=255)
                                                     cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGTAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGtt                                                                                                                                                                         >  1:144470/1‑239 (MQ=255)
                                                     cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCCCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTG‑ACCTAACCCTGCTGtt                                                                                                                                                                         >  1:670087/1‑239 (MQ=255)
                                                          aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGGGAATAGCCTGCGGGCTGAGGGGGAAATGCCCTACCAGCCCTTCGCGACTG‑ACTAAACCCTGCTGTTTAAcc                                                                                                                                                                    >  2:561793/1‑239 (MQ=255)
                                                          aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTCCGGGCTGCCGGAGAAAGGCTCTCCCAGCCCTCCGGTACTG‑AAAAAACCCCGCTGTTTAAca                                                                                                                                                                    >  2:495842/1‑238 (MQ=255)
                                                             tatgggggTTTGGCCCCGGCAGTTGGCGGGCGATGCGTTTTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAg                                                                                                                                                                  <  1:561793/234‑1 (MQ=255)
                                                             ctttGGGTTTTTGCCCCGCCATTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAg                                                                                                                                                                  <  1:573447/236‑1 (MQ=255)
                                                               ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCCACTCTTCGGGACGG‑CCAAAACCCGCCCGTTTAACCCGTTc                                                                                                                                                               >  1:419222/1‑239 (MQ=255)
                                                               ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGGGAAATGCTCTACCAGCTCTTCGGTACTG‑ACCTAACCCTGCTGTTTAAc                                                                                                                                                                     >  1:525607/1‑233 (MQ=255)
                                                               ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGACCACGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAACAGTCTGCGGGCGGGCGGCAAAATGCGCTACCATCTCTTCGGGACGT‑AAATAACCCTGTTGTTTAACCCGt                                                                                                                                                                 >  2:758610/1‑237 (MQ=255)
                                                                     ttttGCCCCCCCATTTGGGGGGCCGTCCCCTTTTTTGCTTGGGTCGGGGGGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAAGTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGataa                                                                                                                                                          <  2:894732/237‑1 (MQ=255)
                                                                      tttcccccGCCCTTTTGGGGGGGATGCGTTCTTTTTCTTCGGTCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGataata                                                                                                                                                        <  1:775097/233‑1 (MQ=255)
                                                                       ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCGGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCACA‑AACCCCGCTGTTTAACCAGttggtaaaaa                                                                                                                                                        >  2:158047/1‑230 (MQ=255)
                                                                         gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCACTTCGGTACTG‑ACATAACCCTGCTGGTTAACCAGTTCGTTAATACCg                                                                                                                                                     >  1:387641/1‑239 (MQ=255)
                                                                              gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt                                                                                                                                                 >  1:761730/1‑238 (MQ=255)
                                                                               ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGACCTG‑ACATAACCCAGCTGTTTAACAAGTTCGAAAAAAACGtttggg                                                                                                                                               >  2:417764/1‑235 (MQ=255)
                                                                                          gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGggtttg                                                                                                                                    >  1:95985/1‑239 (MQ=255)
                                                                                             gtGCGCTCTTTTTCCTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGGTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                                                                  <  1:78239/237‑1 (MQ=255)
                                                                                                     tttttCTTGGGTCGGGGTCCCCTTGCGCTCGGGCCGCGGCTTTTTAACCGAGCTGGAGGGCAGCCCCGCACCGCCATGATGG‑CGCAAGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                                                         <  1:32742/239‑1 (MQ=255)
                                                                                                        ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGTTCACGGAGAAATGCTCTACCAGCTCTCCGGTACTA‑ACAGAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGGCGGTGGGTTTGTTTCATTTATaaa                                                                                                                       >  2:203004/1‑238 (MQ=255)
                                                                                                               gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTTTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCCGTTCGAAATTACCGTTTTGACGTTGTGTTTTTTTCATTTATAAATCCCTAGa                                                                                                               >  2:228747/1‑239 (MQ=255)
                                                                                                                     gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTT‑ACCTAACCCTGCTGTTTAACCCGTTCGATAATACCTTTGTGCCGTTGTTTTTTTTTCTTTTTAAAAACCCTGGActtatt                                                                                                         >  2:595521/1‑239 (MQ=255)
                                                                                                                              cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTGTGTTTGTTTCATTTCTAAATCCCTGGAATTAttttggtttt                                                                                                 >  1:897123/1‑238 (MQ=255)
                                                                                                                                agtgggCCCCGGCGTTTTAACCGAATGGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCCGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcattttcg                                                                                               <  1:492085/233‑1 (MQ=255)
                                                                                                                                   gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCa                                                                                           >  1:712787/1‑239 (MQ=255)
                                                                                                                                          ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGGGAAATGCTCTACCAGCTCTTCGGTACTG‑CCCTAACCCTGCTGTTTAACCAGTTCGGAAATACCGTTGTGACGTTGTGTTTGTTTCTTTTATAAATCCCTTGCATTTTTTTCGTTTTAGCGCCTTGAGCg                                                                                    >  2:605956/1‑239 (MQ=255)
                                                                                                                                             gtgtTAGCCGAGTGGGAGGGCGGCGCCCCGCCGCCGTGATGG‑CGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                                                                 <  2:702402/239‑1 (MQ=255)
                                                                                                                                              tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCGTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCACCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTTGTTCATTTATAAATCCCCTGAATTATTTTCGTTTTCTCGCaatatggcaaca                                                                                >  2:342309/1‑228 (MQ=255)
                                                                                                                                                      gAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTTTAAATCCTTGGAATTATTTTCGTTTTCCCCCTTTGGGCGAATCAACaaaa                                                                         >  1:347825/1‑238 (MQ=255)
                                                                                                                                                        acgggAAGGCAGGGCCGCACCGCCATGATGG‑CGCAGGATTAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:723918/236‑1 (MQ=255)
                                                                                                                                                        acgggAAGGCAGCGCCGCACCGCCATGATGGTG‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:495842/236‑1 (MQ=255)
                                                                                                                                                        accgggaGGCAGGGCCGCGACGCCATGATGG‑CGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCCCAGCTCTTCGGTACTA‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:203004/233‑1 (MQ=255)
                                                                                                                                                        aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTTTGACGTTTGGCTTGTTTCATTTATAAATGCCTGGAATTATTTTCGTTTTCCGTCATTTAGCGCACAAACaaaaac                                                                       >  2:35107/1‑236 (MQ=255)
                                                                                                                                                        aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGAACCG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTTTGTTTGTTTCCTTTATAAATCCCTCGCAATATTTTCGTTTTCTCGCATTTAGCAAAACAaaaaaaaca                                                                      >  2:812647/1‑236 (MQ=255)
                                                                                                                                                        aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTG‑ACATAACCCTGCTGGTTAACCAGTTCGTTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTTGGCCATTGAGGGAATAAACAAAAGcc                                                                      >  2:250260/1‑239 (MQ=255)
                                                                                                                                                              aGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattatt                                                                                                         >  1:189570/1‑198 (MQ=255)
                                                                                                                                                                   cccccgccccgccATGTTGG‑CGGAGTTTCAGGTTTTGCTGCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcc                                                           <  2:875439/236‑1 (MQ=255)
                                                                                                                                                                    cgccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcca                                                          >  1:45437/1‑238 (MQ=255)
                                                                                                                                                                      ccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                                                    <  2:713154/179‑1 (MQ=255)
                                                                                                                                                                      ccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                                                    >  1:713154/1‑179 (MQ=255)
                                                                                                                                                                         cACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                                                       >  1:597298/1‑237 (MQ=255)
                                                                                                                                                                           cccccATGATGG‑CGCAGGTTCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa                                                   <  2:187778/236‑1 (MQ=255)
                                                                                                                                                                                  gatgG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGc                                            <  1:417764/239‑1 (MQ=255)
                                                                                                                                                                                     tgggtcAGGATTAGGTTTTGGCTCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTG‑ACATTACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGa                                          <  1:783059/234‑1 (MQ=255)
                                                                                                                                                                                        cgcAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGaggg                                       >  1:731310/3‑236 (MQ=255)
                                                                                                                                                                                          cAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGTtgt‑gt                                    >  1:298190/1‑239 (MQ=255)
                                                                                                                                                                                          cAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGCCTAACCCGGCGTtgt‑gt                                    >  1:811729/1‑239 (MQ=255)
                                                                                                                                                                                                  ggTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCCCTAAACCCACGGCTAACCCGGCCAGGG‑TTGCTCCGt                             >  2:899036/1‑238 (MQ=255)
                                                                                                                                                                                                       tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGACTCCACAAAACCCATCGCTAAACCCACGGCTAACCC‑GCCACGTGGTGCTCCGTTCGcac                       >  2:723354/1‑237 (MQ=255)
                                                                                                                                                                                                       tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGcc                        >  1:528488/1‑238 (MQ=255)
                                                                                                                                                                                                          tCCGCCTGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCg                    <  1:228747/239‑1 (MQ=255)
                                                                                                                                                                                                          tCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAGCCCTGCTGTTTAGCCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCg                    <  1:638399/239‑1 (MQ=255)
                                                                                                                                                                                                           ccGGCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCGa                   <  1:158047/239‑1 (MQ=255)
                                                                                                                                                                                                             gCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGAATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGACGT‑GTGCTCCGTTCGCCATCgcc                  >  1:785332/1‑236 (MQ=255)
                                                                                                                                                                                                              ccAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAATCCATCGCTAAAACCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCTACa                 >  1:747595/1‑238 (MQ=255)
                                                                                                                                                                                                               cAGCTCATTGAGGTCGCCGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCGACATc               <  1:672063/239‑1 (MQ=255)
                                                                                                                                                                                                                           gTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGt                                                                                                                                                                 >  1:8529/1‑82 (MQ=255)
                                                                                                                                                                                                                           gTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGt                                                                                                                                                                 <  2:8529/82‑1 (MQ=255)
                                                                                                                                                                                                                            tCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAACCCATCGCTAAACCCACGCCTACCCCGGCGATGT‑GTGCTCCCTTCCCCATCGACATCCCAAACAAATCaa  >  2:438258/1‑239 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
TTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTG‑ACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT‑GTGCTCCGTTCGCCATCGACATCCCAAACAAATCAA  >  NC_000913/3560272‑3560728

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