Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 550,033 A→G 14.3% E136G (GAG→GGG)  ylbF → putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913550,0330AG14.3% 163.7 / 2.4 73E136G (GAG→GGG) ylbFputative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base G (4/7);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTAT  >  NC_000913/549800‑550262
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
tGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGAACCGCTGCCTGCT‑GGGTTa                                                                                                                                                                                                                                    >  1:406373/1‑238 (MQ=255)
  ggcAATTCGGTTGGTGATGCGGTGGATTATGCGGAAATGAACGACCACAGGTGGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACg                                                                                                                                                                                                                                  <  1:864997/237‑1 (MQ=255)
            tGCGTTATGCGGTGGATTATGCGGCAATGAACGCCCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGCGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCCCACGCTGCT‑GAGTTACGCCAGGTTCgga                                                                                                                                                                                                                       >  2:54259/1‑237 (MQ=255)
            tGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GGGTTACGCCCGtttggg                                                                                                                                                                                                                        >  1:719880/1‑235 (MQ=255)
               gTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACt                                                                                                                                                                                                                    >  1:467708/1‑239 (MQ=255)
                 tATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCCCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGACACCGGGCTTTTCGGGCCACTGGCGTTGCCGGCAAGCGACCCGCTGACTGCT‑GCGTTACGCCAGTTTCGTCCCtct                                                                                                                                                                                                                  >  1:274191/1‑237 (MQ=255)
                        tGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACAACTGGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGCCATTTCTGTCCCCTGGCATGTGCGGCAAGCGACTCGCTCCCTGCT‑GCTGTACGCCAGTTTCCTAACTGTTTTCAg                                                                                                                                                                                                            >  2:863843/1‑238 (MQ=255)
                          gATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGAGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATGCGCTGCCTGCT‑GGGTTACGCCAGTTTCGTCACTGTTTTCAGGcc                                                                                                                                                                                                         >  2:731883/1‑239 (MQ=255)
                               tGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCGGCCTGCT‑GGGGTACGCCCGTTTCGTCCCTGGTTTCAGGCCccgc                                                                                                                                                                                                     >  1:285960/1‑238 (MQ=255)
                                    cAATGACCCACCACAGGTTTGGGCTCAAGGGATTCCCCTCGGGCG‑TTTCGCGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAAt                                                                                                                                                                                                <  1:366350/238‑1 (MQ=255)
                                    cAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCCTGGATGCAACGCGCGGAGGCAACCGGGCTTTTCGGGCCACTGGCGTTGGCGCCAGGCGATCCGCGGCCGCCT‑GAGTTGCGCCCGTTTCTTCCGTGTTTCAAGGCAGCGCTCAATg                                                                                                                                                                                               >  2:307734/1‑239 (MQ=255)
                                       agagcacgtacaGTTTGTGTCGAAAGGGGTTCGCC‑CGGGCGTTTTCACGGTTAATCCGCCGCGGGGTTATTGTGTGGTGAGTTTTCCGCCGCCTGGGCATCCGCAACCCCTGGCAGGTGCATGGATGCAACGCGCGGATGAAACCGGGGTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCg                                                                                                                                                                                            <  2:763199/230‑1 (MQ=255)
                                                acaGGTTGTGGCTCAAGGGATTCGCCTGGGGGG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCg                                                                                                                                                                                   <  1:806665/239‑1 (MQ=255)
                                                      tgtgGCCCGAGGGTTTCGCCTAGGGCG‑TTTCCCGGTTAAACAGCCAAGGCGTTATGTGTTGCTGCGTATTACCCCGCCGGCGCATCCTCAACCACGTGCAGCTGCAGGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCGGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTgg                                                                                                                                                                              <  2:341949/238‑1 (MQ=255)
                                                        tgGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCGGCT‑GAGTGACGCCAGTTTCGTCACTGTTTTCAGGCCCAGCCCACTGGGGTTAAGACCGACggggg                                                                                                                                                                            >  2:685415/1‑236 (MQ=255)
                                                        tgGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAACCTATCCGCTGCCTGCT‑GCGTTACGCCAGTTTCGTCACTGGTTTCAGGCCGCGGTCAATGGGTGTAAGAACGACTggcgg                                                                                                                                                                           >  2:222960/1‑239 (MQ=255)
                                                          gCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTGTTCAGGCCGCGCTCAATGGCGTTCAGACCGa                                                                                                                                                                                  >  1:763199/1‑230 (MQ=255)
                                                              aaGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGCGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTggcggaccgg                                                                                                                                                                      >  1:339212/1‑233 (MQ=255)
                                                                      cGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCT‑GGGTaa                                                                                                                                                            >  1:612495/1‑239 (MQ=255)
                                                                            gggCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCATCTGCATGGCTGCAACGCGCGGCGGAAACCGGGCTTTTCGGACCACTGGCGTTGGCGGCAAGCAATCCACTGCCTGCT‑GAGTTACGCCGTTTTCGTCCCTCTTTTCCGGCCGGCCTCAAAGGCGTTAAGCCCTACTGGCGCCCATTGCT‑GGGAGAAggga                                                                                                                                                       >  2:653322/1‑236 (MQ=255)
                                                                            gggCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACGTTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCG‑GGGTAAAGGcc                                                                                                                                                       >  1:556848/1‑238 (MQ=255)
                                                                                  tttCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGAt                                                                                                                                                 >  1:816171/1‑239 (MQ=255)
                                                                                   ttCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCCCAATGGCGTTAAAACCGACTGGCGGCACTGGCG‑GGGTAAAGGCCCCGGGtt                                                                                                                                                >  2:653754/1‑239 (MQ=255)
                                                                                        ggTTAAACAGCCAAAGCGTTATTGTTTGGTGCGTATTACGCCGCCTGCGCAATCTCAACCACTTGCAGCTGGATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGc                                                                                                                                           <  1:569619/239‑1 (MQ=255)
                                                                                         gTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCT‑GGGGAAAGGCCCCGGGTTAAcccc                                                                                                                                          >  1:916114/1‑239 (MQ=255)
                                                                                          ttAAACCGCCCAAGTTTTTTTGTTTGTCGGGTATTAGGGCGCCGGGGCACCCTAAACCTCTTGAAGGTGCATGGATGGAACGCGGGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGcc                                                                                                                                          <  2:406373/238‑1 (MQ=255)
                                                                                          atAAACGGCCTCGGGGTTGTTGTTTGCTGGGTATGACGCCGCCTGAGCATGCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAGCAGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGc                                                                                                                                           <  1:489595/236‑1 (MQ=255)
                                                                                                         ttttttttttGCTTCGTTTTCCGCCGCCGGGGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACg                                                                                                                          <  2:816171/232‑1 (MQ=255)
                                                                                                                  tgctgcGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACCTGCT‑GGGTAAAGGCCCCGGGTTAACACCGGGTCGTGATGAACAGCCGAGCGGa                                                                                                                 >  1:765411/1‑239 (MQ=255)
                                                                                                                   gctgcGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAgg                                                                                                                                                                                                           <  1:606928/149‑1 (MQ=255)
                                                                                                                   gctgcGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAgg                                                                                                                                                                                                           >  2:606928/1‑149 (MQ=255)
                                                                                                                       cGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCcgcg                                                                                                                                                                                                      >  1:91424/1‑150 (MQ=255)
                                                                                                                       cGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCcgcg                                                                                                                                                                                                      <  2:91424/150‑1 (MQ=255)
                                                                                                                          aTTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCT‑GGGGAAAGGCCCCCGATTAAAGCCCAGACATGGTGCAACCCTGGGGGGaaacacgcc                                                                                                         >  2:488704/1‑232 (MQ=255)
                                                                                                                            tcgcgcgccTGCGCATCATCAACCGTGTGAAGCTGCAGGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACGGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTgg                                                                                                       <  1:834506/234‑1 (MQ=255)
                                                                                                                                 cgccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGcc                                                                                                  <  1:653754/239‑1 (MQ=255)
                                                                                                                                 ccccTGCGCCTCCGCCACCCTTTGCAGCTGCGTGGATGCAACGCGGGGAGGAAACCGGGCTTTTCGGGCCCCTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGcc                                                                                                  <  2:285960/237‑1 (MQ=255)
                                                                                                                                   ccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCT‑GGGGAAAGGCCCCGGCTTAAACCCGGGCCATGATGAAACCCTGAGCGGAAAGCTGCGGGGGGcccg                                                                                                >  2:380126/1‑237 (MQ=255)
                                                                                                                                   ccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGTCT‑GGGGAAAGGCCCCGGATTAAAGCCGAGGCCTGGTGACAAGCTGGGCGGAATGCGgcgggcgggcgg                                                                                                >  2:371294/1‑235 (MQ=255)
                                                                                                                                      gcgcATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAAGCCCCGGATTAAAGCCGAGTCATGATGACACCCCGAACGGAATGCTGCTGGGGGCCTggga                                                                                             >  1:897431/1‑237 (MQ=255)
                                                                                                                                            ccTCAACCTATTGCAGCTGCAGGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATgg                                                                                       <  1:697277/239‑1 (MQ=255)
                                                                                                                                            ccTCAACCACTTGGAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATgg                                                                                       <  1:488704/239‑1 (MQ=255)
                                                                                                                                                 aCCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGCTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGGATTATGGCGCtt                                                                                  >  1:410714/1‑239 (MQ=255)
                                                                                                                                                       tGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCACCTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATg                                                                            <  1:784431/239‑1 (MQ=255)
                                                                                                                                                            cTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGCCACGCTGAGCGGAATGCTGCTGGCGGCCTGGGATTATTGCGCTTTAGATgcgcg                                                                        >  1:835721/1‑238 (MQ=255)
                                                                                                                                                                aTGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTTACGCCGAGTCATGCTGACACCCTCACGGGAATGTGGCCGGCGGCCCGGGTTTATGGCGCTTTAGTGGCGCGCCCCg                                                                   >  2:54996/1‑239 (MQ=255)
                                                                                                                                                                         aCGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACACCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGCATTATTGGGCTTTAGATGCGCGCTCCGGTGTTCCGt                                                          >  1:602025/1‑239 (MQ=255)
                                                                                                                                                                            cgcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTa                                                                                                                                                             <  2:486901/137‑1 (MQ=255)
                                                                                                                                                                            cgcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTa                                                                                                                                                             >  1:486901/1‑137 (MQ=255)
                                                                                                                                                                              cgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGCTTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGACTGGTATTATGGCGCTTTAGATGGGCGCCCCGGTCGGCCGgttttt                                                     >  2:410472/1‑239 (MQ=255)
                                                                                                                                                                                gaggaAACCGGGCTTTTCGGCCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTTCGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTg                                                    <  1:805877/238‑1 (MQ=255)
                                                                                                                                                                                 aggaAACCGGGTTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGc                                                   <  1:380126/238‑1 (MQ=255)
                                                                                                                                                                                         ggggTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GGGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGa                                          <  1:685415/239‑1 (MQ=255)
                                                                                                                                                                                         gggTTTTTCGGGCCCCTGGGGTTGGCGGCAAGCGTTCCGCTGCGTGCT‑GAGTTACGCCAGTTTGGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGCCCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGa                                          <  2:916114/239‑1 (MQ=255)
                                                                                                                                                                                           gCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCGCCGGTCGGCCGTTTTTTGCCTGTTCCGACa                                        >  1:263733/1‑239 (MQ=255)
                                                                                                                                                                                           gCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGGCATGATGACACGCTGAGCGGGATGCTGCTGGCGGCCTGGTATTATGGGGCTTTAGATGCCGGGCCCGGTCGCCCGTTTTTTGCCTGGTCCGACa                                        >  2:19802/1‑239 (MQ=255)
                                                                                                                                                                                            gttttGGGGCCCCTGGGGGTGGGGGCAAGGGATCCGGTGCCTGTT‑GAGTTACGGCAGTTTCGTCACTGGTTTCAGGCCGCGCTTAATGGCGTTAAGACCGAATGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACa                                        <  2:467708/237‑1 (MQ=255)
                                                                                                                                                                                            cTTTTCGGGCCACTGGCGTTGGCGGCAAGCGAGCCGCTGCCTGCT‑GAGGTACGCCAGTTTCGTCACTGTGTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGAGTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGTCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACa                                        <  1:692186/238‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCCGCTGGCGTTGGGGGGAAGCCATCCGCTGCCTGCT‑GAGTTTCGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGGCCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGc                                  <  1:547152/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCCCCTGGCTTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGc                                  <  1:715734/239‑1 (MQ=255)
                                                                                                                                                                                                   ggcccatGGCTTTGGGCGCAAGCGATCCGCTGCCTGCT‑TGTTTCCGCCAGTTTCGTCACGGTTTTCCGGCCGCGCTCAATGGGTTTAAGCACGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  2:556848/233‑1 (MQ=255)
                                                                                                                                                                                                    gcagcTGGCGTTGGCGGCAAGGGATCCGCTGCCTGCT‑GGGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:584693/234‑1 (MQ=255)
                                                                                                                                                                                                    gCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGGTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:876113/238‑1 (MQ=255)
                                                                                                                                                                                                    gCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAATTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:730093/238‑1 (MQ=255)
                                                                                                                                                                                                     ccACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GGGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGGTAAAGGCCCCGGATTAACGCCGCGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:87820/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCCATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTg              <  1:210500/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTg              <  1:19802/238‑1 (MQ=255)
                                                                                                                                                                                                                           ttCCGCTGGCTGCT‑GAGTTGCGCCAGTTTCGTCAGTGTTTTCAGGCCGCGGTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGt         <  1:797230/237‑1 (MQ=255)
                                                                                                                                                                                                                           ttCCGCTGCCGTCT‑GGGTTTCGCCCGTTTGGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc        <  2:765411/238‑1 (MQ=255)
                                                                                                                                                                                                                           aTCCGCTGCTTGCT‑GGGTTGGGCCAGTTTCGGCACTGTTTTCGGGGCGCGCTCAATGGCGTTAAGACCGACTGGCGGGACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGCCACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc        <  2:624210/239‑1 (MQ=255)
                                                                                                                                                                                                                            tCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc        <  1:757587/238‑1 (MQ=255)
                                                                                                                                                                                                                             ccGCTGCGTGCTGGGGTTAC‑CGAGTTGCGTCCCTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAg      <  2:719880/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tggctgct‑gAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGGTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTAt  <  1:371294/236‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
TGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCG‑TTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCT‑GAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCT‑GGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTAT  >  NC_000913/549800‑550262

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: ‑