Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,494 G→A 6.9% 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,4940GA 6.9% 183.6 / 4.5 72pseudogene (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 (47/19);  new base A (3/2);  total (50/22)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAA  >  NC_000913/3560266‑3560728
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
cTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGGGGGGAAAAg                                                                                                                                                                                                                                     >  2:378709/1‑238 (MQ=255)
  gtgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGGGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                                                                                                                                                     >  1:840054/1‑236 (MQ=255)
          gCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGGGGTCGCGGCGAAAAGTCTGCGGGCTg                                                                                                                                                                                                                          >  1:363446/1‑239 (MQ=255)
              tgcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑GCCC‑CGGGGTCAGGTTTTGCTCCGCCAGCCCATTGGGGGCCCGGCGAATAGTGTGGGGGCTGACgg                                                                                                                                                                                                                      >  2:588318/1‑239 (MQ=255)
               gcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGCGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCCTGATGG‑GG‑CAGGGCCAGGTTTTGCCCCGCCAGGCCATTGGGGTCGCGGCGAAACGTCTGGGGGcgtggggg                                                                                                                                                                                                                     >  2:119791/1‑232 (MQ=255)
                cgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGCAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCACCCCATTGAGGTCCCGGCGAATAGTCTCCGGGCTGACGgag                                                                                                                                                                                                                    >  1:405478/1‑239 (MQ=255)
                cgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                                                                    >  1:1032872/1‑239 (MQ=255)
                cgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCACCG‑CCG‑CCCCGCCATGAT‑G‑CCGCGGGATAAGGTTTTGCTCCGCCACCGCATTGAGGCCGCGGCGAATAGCCTGGGGGCTGGCGgag                                                                                                                                                                                                                    >  2:712952/1‑239 (MQ=255)
                   gcgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                                                       >  1:432817/1‑233 (MQ=255)
                        aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCGGCCCATTAAGGTCGCGGCGAAAAGTCTGCGGGCTGCCGGAGCAATGc                                                                                                                                                                                                              >  1:593407/1‑237 (MQ=255)
                        aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGGGGGCTGACGGGGAAATGc                                                                                                                                                                                                              >  1:1052287/1‑237 (MQ=255)
                        aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGGGAAAAGTCTGCGGGCGGCCGGGGAAAGGcc                                                                                                                                                                                                             >  2:1058587/1‑237 (MQ=255)
                            cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCCCCAGCCCATTGAGGTCCCGGGGAAAAGTTTGGGGGGTGAGGGGGAAAGTCTCCAcc                                                                                                                                                                                                        >  2:765823/1‑239 (MQ=255)
                               ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTACTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGCGGAATGCTATACCa                                                                                                                                                                                                       >  1:160715/1‑237 (MQ=255)
                                    ggggggcgtatatCAAAAAACGGGGTCGAGCCTTTGGGGTTTTG‑CCCCCCCCCTTGGCGGGCAATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTG‑CCACGGGGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                                                 <  1:639033/229‑1 (MQ=255)
                                       gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAATGCAGCG‑CCG‑CCCCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAAAATTCTCCGGGCTGCCGGGGAAATACTCTACAACCACTTCggg                                                                                                                                                                                             >  2:457801/1‑238 (MQ=255)
                                         gCCAATATGGAAAAACGGCGTCGGGCCTTTGGGGATTTC‑CTCCGCCACTTGGCGGGCGAGGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                                                                                                                                            <  1:765823/238‑1 (MQ=255)
                                             atatCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGGGAACTGCTCTACCAGCCCTTCGGTCCTgcc                                                                                                                                                                                       >  1:802716/1‑237 (MQ=255)
                                                 cGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCCCTACCAGCTCTTCGGTACTGACATa                                                                                                                                                                                    >  1:1007436/1‑238 (MQ=255)
                                                  gATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTTGGGTCGCGGCGAATAGTCTGCGGGCTGACGGGAGAATGCTCTACCAGCTCTTCGGAACTAACATaa                                                                                                                                                                                   >  1:691442/1‑238 (MQ=255)
                                                  gATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTGCCATaa                                                                                                                                                                                   >  1:706372/1‑238 (MQ=255)
                                                      aaCCGCGTCGTGCTTTTGGGGGTTTGTCCCCGCCTCTTGGGGGGGGTGGCGCTCTTTTTCTGCGGGGTGGGGGGCCTTGCGATCGTG‑CCACGGCGTGTAAACCGAACTGGAAGGCAGCG‑CCG‑CCCCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                                               <  2:672874/239‑1 (MQ=255)
                                                           cGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCCCGGCGAATAGTCTGCGGGCTGACGGGGAAATGCGCTACCACCTCTTCGGTCCTGACCTAACCCTGCTGtt                                                                                                                                                                         >  1:1065782/1‑239 (MQ=255)
                                                             tCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                                                                                                                                            >  1:181856/1‑218 (MQ=255)
                                                             tCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTCCGGGCTGCGGGAGAAATGCCTTACCCGCTCTTCGGCACTGACATAACCCTGCTGTTTa                                                                                                                                                                       >  2:251541/1‑239 (MQ=255)
                                                                aGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTCCTGTTTAAcc                                                                                                                                                                    >  1:699526/1‑239 (MQ=255)
                                                                   ctttGGGGTTTGT‑CCCCGCCCGTTGGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGGGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                                                  <  1:601828/236‑1 (MQ=255)
                                                                            tttG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGGACTGCCATAACCCTGCTGGTTAAACAGTTCGaaaaa                                                                                                                                                         >  2:1003729/1‑237 (MQ=255)
                                                                            tttG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGAGAAAGCCTCTACCGGCCCTTCGGAACTGACATAACCCTGCAGTTTAACCAGTTCGttaat                                                                                                                                                         >  1:1067179/1‑238 (MQ=255)
                                                                             ttG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTTACATAAACCCGCGGTTTAACCAGGTCGaaaata                                                                                                                                                        >  2:7204/1‑238 (MQ=255)
                                                                             ttG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTTACAAAACCCTGCTGTTTAACCTGtctcttataaa                                                                                                                                                       >  2:927307/1‑229 (MQ=255)
                                                                                      ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCCGTAATACCGTtctg                                                                                                                                                >  1:693774/1‑238 (MQ=255)
                                                                                         gtttGGGGGGGATGCGCTCTTTTTCTTGGGTCTGGGTGGCCTTGCGATCGGG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                                                            <  1:215822/238‑1 (MQ=255)
                                                                                         cttggggggCGATGCGCTTTTTTTCTGCGGTCTGGTGGGCCTGGCGATCGTC‑ACACGGCGTGTTAACCGAAGTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                                                            <  1:470762/233‑1 (MQ=255)
                                                                                            tGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCAAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGGATAGTCTGCGGGCTGACGGAGGAATGCTCTACCCCCTCTTCGGTACTGCCCTAACCCTGCTGTTTAACCAGTTCGATAAAACCGTTGTGACGtttg                                                                                                                                         >  2:877474/1‑237 (MQ=255)
                                                                                                   gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCCGCTGTTTAACCAGGTCGTAAAACCCGTTGGGACGTTGgggttgtt                                                                                                                                  >  2:788624/1‑239 (MQ=255)
                                                                                                       ccctttttttCTGGGGTCGGGGGGGCTTTGCG‑GCGTGCCCAGGGCG‑GTGAACCCGACTGGGAGGGCGCGCCCG‑CACCGCCATGAT‑G‑GCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCCCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAt                                                                                                                              <  2:405478/234‑1 (MQ=255)
                                                                                                             ttttcgcgggTCGGGGTGGCCTGGCGATCGTG‑CCGCGGCTTGTTAACCGACCGGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                                                        <  1:119791/231‑1 (MQ=255)
                                                                                                                  tgcGGTTTGGGTGGCTTTGGGCTCGTG‑CCAGGGCGTGTTAACCGAGCTGGAAGGCGGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCGGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                                                                   <  2:322936/237‑1 (MQ=255)
                                                                                                                    cGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGGT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAAAACCGTTGTGACGTTGTGTTTTTTTCATTTTTAAATCCCTg                                                                                                                 >  2:413580/1‑239 (MQ=255)
                                                                                                                      gTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTTTACCAGCCCTCCGGAACTGACATAACCCTGCTGTTTAACCAGGTCGATAAAACCGTTGGGACGTGTGGTTTTTTTCATTTTAAAATCCCCGGa                                                                                                               >  2:353030/1‑239 (MQ=255)
                                                                                                                        cggggTGGCCTTGGAATGGTG‑CCAGGGGGTGTTACCCGAACGGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGaa                                                                                                              <  1:788624/236‑1 (MQ=255)
                                                                                                                                ccTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTTTAAATCCCTGGAATTAttttcg                                                                                                     >  1:701772/1‑239 (MQ=255)
                                                                                                                                     cGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttc                                                                                                >  1:396645/1‑239 (MQ=255)
                                                                                                                                         cGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGAAAAAACCGTTGTGACGTTGGGGTTGTTTAATTTTTAAATCCCCGGAATAATTTTCGTTTTCgcgc                                                                                            >  2:875367/1‑239 (MQ=255)
                                                                                                                                          gTG‑CCACGGGGTGTTACCCGAACGGGAAGGCAGGG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCa                                                                                           <  1:378709/239‑1 (MQ=255)
                                                                                                                                          gTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGGGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCa                                                                                           >  1:545583/1‑239 (MQ=255)
                                                                                                                                                      ggtTAACCGAACGGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                                                                 <  1:7204/237‑1 (MQ=255)
                                                                                                                                                       gtTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTTGGACGTTGTGTTTGTTTCATTTTTAAATCCCTGGAATTATTTTCGTTTTCGGCCATTTAGCGAATca                                                                               >  2:682661/1‑239 (MQ=255)
                                                                                                                                                              gAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCCTTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                                                                         >  1:838455/1‑238 (MQ=255)
                                                                                                                                                                atgggAAGGCAGGC‑CCG‑CACCGCAAGTAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGGATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  1:877474/236‑1 (MQ=255)
                                                                                                                                                                  ttgggggCAGCC‑CCGCCCCCGCTCTGAT‑G‑GCGAAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGGTGGGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                                                      <  2:160715/233‑1 (MQ=255)
                                                                                                                                                                          aGCG‑CCG‑CACCGCCATGAT‑G‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaaa                                                             >  1:84295/1‑238 (MQ=255)
                                                                                                                                                                           ccg‑ccc‑cCCCGCCATGGT‑G‑GCGCAGGAACAGGTTTGGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑cc                                                          <  2:1052287/238‑1 (MQ=255)
                                                                                                                                                                                              gattgcgcaggttcaggttTTGCTCCCCCAGCCCATTGAGGTCCCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGc                                           <  2:363446/232‑1 (MQ=255)
                                                                                                                                                                                               atggcgcaggttcaggttTTGGCCCGCCAGCCCATTGAGGGCGCGGGGGATAGTCTGGGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGc                                           <  2:1032872/232‑1 (MQ=255)
                                                                                                                                                                                                  gcgcAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                                                                   <  1:856694/209‑1 (MQ=255)
                                                                                                                                                                                                  gcgcAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                                                                   >  2:856694/4‑212 (MQ=255)
                                                                                                                                                                                                     cAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCCCCGGCTAACCCGGCGAtgtg                                     >  2:460648/1‑239 (MQ=255)
                                                                                                                                                                                                       ggATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGAtgtgtg                                   <  1:251541/238‑1 (MQ=255)
                                                                                                                                                                                                                  tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                                                 <  2:276088/70‑1 (MQ=255)
                                                                                                                                                                                                                  tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                                                 >  1:276088/1‑70 (MQ=255)
                                                                                                                                                                                                                  tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCTTTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTCTGCTCCGTTCGcc                        >  1:328972/1‑238 (MQ=255)
                                                                                                                                                                                                                  tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGcc                        >  1:1051126/1‑238 (MQ=255)
                                                                                                                                                                                                                  tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGcc                        >  1:27158/1‑238 (MQ=255)
                                                                                                                                                                                                                  tGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGcc                        >  1:190658/1‑238 (MQ=255)
                                                                                                                                                                                                                    cTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAATCCATCGCTAAA‑CCCCAGGCTAACCCGGCGATGTGTGCTCCGTTCGCCCTc                     >  1:951271/1‑239 (MQ=255)
                                                                                                                                                                                                                    cTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑ACCACGCCTAACCCGGCGATGTGTGCTCCGTTCCCCATc                     >  1:21907/1‑239 (MQ=255)
                                                                                                                                                                                                                        gCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCg                                          <  2:174460/214‑1 (MQ=255)
                                                                                                                                                                                                                        gCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCg                                          >  1:174460/1‑214 (MQ=255)
                                                                                                                                                                                                                               taTTGAGGTCGCGGCGAAATGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTTTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCcaaa          <  1:1003729/238‑1 (MQ=255)
                                                                                                                                                                                                                                      gTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCAATCGACATCCCAAACAAATCa   >  1:668126/1‑239 (MQ=255)
                                                                                                                                                                                                                                        cgcgGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGCCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCaa  >  1:304270/1‑238 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCG‑CCG‑CACCGCCATGAT‑G‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAA‑CCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAA  >  NC_000913/3560266‑3560728

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