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

TGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTAT  >  NC_000913/549800‑550262
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
tGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGAACCGCTGCCTGCTGGGTTa                                                                                                                                                                                                                                    >  1:406373/1‑238 (MQ=255)
  ggcAATTCGGTTGGTGATGCGGTGGATTATGCGGAAATGAACGACCACAGGTGGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACg                                                                                                                                                                                                                                  <  1:864997/237‑1 (MQ=255)
            tGCGTTATGCGGTGGATTATGCGGCAATGAACGCCCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGCGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCCCACGCTGCTGAGTTACGCCAGGTTCgga                                                                                                                                                                                                                       >  2:54259/1‑237 (MQ=255)
            tGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGGGTTACGCCCGtttggg                                                                                                                                                                                                                        >  1:719880/1‑235 (MQ=255)
               gTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACt                                                                                                                                                                                                                    >  1:467708/1‑239 (MQ=255)
                 tATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCCCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGACACCGGGCTTTTCGGGCCACTGGCGTTGCCGGCAAGCGACCCGCTGACTGCTGCGTTACGCCAGTTTCGTCCCtct                                                                                                                                                                                                                  >  1:274191/1‑237 (MQ=255)
                        tGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACAACTGGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGCCATTTCTGTCCCCTGGCATGTGCGGCAAGCGACTCGCTCCCTGCTGCTGTACGCCAGTTTCCTAACTGTTTTCAg                                                                                                                                                                                                            >  2:863843/1‑238 (MQ=255)
                          gATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGAGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATGCGCTGCCTGCTGGGTTACGCCAGTTTCGTCACTGTTTTCAGGcc                                                                                                                                                                                                         >  2:731883/1‑239 (MQ=255)
                               tGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCGGCCTGCTGGGGTACGCCCGTTTCGTCCCTGGTTTCAGGCCccgc                                                                                                                                                                                                     >  1:285960/1‑238 (MQ=255)
                                    cAATGACCCACCACAGGTTTGGGCTCAAGGGATTCCCCTCGGGCGTTTCGCGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAAt                                                                                                                                                                                                <  1:366350/238‑1 (MQ=255)
                                    cAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCCTGGATGCAACGCGCGGAGGCAACCGGGCTTTTCGGGCCACTGGCGTTGGCGCCAGGCGATCCGCGGCCGCCTGAGTTGCGCCCGTTTCTTCCGTGTTTCAAGGCAGCGCTCAATg                                                                                                                                                                                               >  2:307734/1‑239 (MQ=255)
                                       agagcacgtacaGTTTGTGTCGAAAGGGGTTCGCCCGGGCGTTTTCACGGTTAATCCGCCGCGGGGTTATTGTGTGGTGAGTTTTCCGCCGCCTGGGCATCCGCAACCCCTGGCAGGTGCATGGATGCAACGCGCGGATGAAACCGGGGTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCg                                                                                                                                                                                            <  2:763199/230‑1 (MQ=255)
                                                acaGGTTGTGGCTCAAGGGATTCGCCTGGGGGGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCg                                                                                                                                                                                   <  1:806665/239‑1 (MQ=255)
                                                      tgtgGCCCGAGGGTTTCGCCTAGGGCGTTTCCCGGTTAAACAGCCAAGGCGTTATGTGTTGCTGCGTATTACCCCGCCGGCGCATCCTCAACCACGTGCAGCTGCAGGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCGGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTgg                                                                                                                                                                              <  2:341949/238‑1 (MQ=255)
                                                        tgGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCGGCTGAGTGACGCCAGTTTCGTCACTGTTTTCAGGCCCAGCCCACTGGGGTTAAGACCGACggggg                                                                                                                                                                            >  2:685415/1‑236 (MQ=255)
                                                        tgGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAACCTATCCGCTGCCTGCTGCGTTACGCCAGTTTCGTCACTGGTTTCAGGCCGCGGTCAATGGGTGTAAGAACGACTggcgg                                                                                                                                                                           >  2:222960/1‑239 (MQ=255)
                                                          gCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTGTTCAGGCCGCGCTCAATGGCGTTCAGACCGa                                                                                                                                                                                  >  1:763199/1‑230 (MQ=255)
                                                              aaGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGCGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTggcggaccgg                                                                                                                                                                      >  1:339212/1‑233 (MQ=255)
                                                                      cGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCTGGGTaa                                                                                                                                                            >  1:612495/1‑239 (MQ=255)
                                                                            gggCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCATCTGCATGGCTGCAACGCGCGGCGGAAACCGGGCTTTTCGGACCACTGGCGTTGGCGGCAAGCAATCCACTGCCTGCTGAGTTACGCCGTTTTCGTCCCTCTTTTCCGGCCGGCCTCAAAGGCGTTAAGCCCTACTGGCGCCCATTGCTGGGAGAAggga                                                                                                                                                       >  2:653322/1‑236 (MQ=255)
                                                                            gggCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACGTTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCGGGGTAAAGGcc                                                                                                                                                       >  1:556848/1‑238 (MQ=255)
                                                                                  tttCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGAt                                                                                                                                                 >  1:816171/1‑239 (MQ=255)
                                                                                   ttCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCCCAATGGCGTTAAAACCGACTGGCGGCACTGGCGGGGTAAAGGCCCCGGGtt                                                                                                                                                >  2:653754/1‑239 (MQ=255)
                                                                                        ggTTAAACAGCCAAAGCGTTATTGTTTGGTGCGTATTACGCCGCCTGCGCAATCTCAACCACTTGCAGCTGGATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGc                                                                                                                                           <  1:569619/239‑1 (MQ=255)
                                                                                         gTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCTGGGGAAAGGCCCCGGGTTAAcccc                                                                                                                                          >  1:916114/1‑239 (MQ=255)
                                                                                          ttAAACCGCCCAAGTTTTTTTGTTTGTCGGGTATTAGGGCGCCGGGGCACCCTAAACCTCTTGAAGGTGCATGGATGGAACGCGGGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGcc                                                                                                                                          <  2:406373/238‑1 (MQ=255)
                                                                                          atAAACGGCCTCGGGGTTGTTGTTTGCTGGGTATGACGCCGCCTGAGCATGCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAGCAGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGc                                                                                                                                           <  1:489595/236‑1 (MQ=255)
                                                                                                         ttttttttttGCTTCGTTTTCCGCCGCCGGGGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACg                                                                                                                          <  2:816171/232‑1 (MQ=255)
                                                                                                                  tgctgcGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACCTGCTGGGTAAAGGCCCCGGGTTAACACCGGGTCGTGATGAACAGCCGAGCGGa                                                                                                                 >  1:765411/1‑239 (MQ=255)
                                                                                                                   gctgcGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAgg                                                                                                                                                                                                           <  1:606928/149‑1 (MQ=255)
                                                                                                                   gctgcGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAgg                                                                                                                                                                                                           >  2:606928/1‑149 (MQ=255)
                                                                                                                       cGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCcgcg                                                                                                                                                                                                      >  1:91424/1‑150 (MQ=255)
                                                                                                                       cGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCcgcg                                                                                                                                                                                                      <  2:91424/150‑1 (MQ=255)
                                                                                                                          aTTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCTGGGGAAAGGCCCCCGATTAAAGCCCAGACATGGTGCAACCCTGGGGGGaaacacgcc                                                                                                         >  2:488704/1‑232 (MQ=255)
                                                                                                                            tcgcgcgccTGCGCATCATCAACCGTGTGAAGCTGCAGGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACGGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTgg                                                                                                       <  1:834506/234‑1 (MQ=255)
                                                                                                                                 cgccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGcc                                                                                                  <  1:653754/239‑1 (MQ=255)
                                                                                                                                 ccccTGCGCCTCCGCCACCCTTTGCAGCTGCGTGGATGCAACGCGGGGAGGAAACCGGGCTTTTCGGGCCCCTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGcc                                                                                                  <  2:285960/237‑1 (MQ=255)
                                                                                                                                   ccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCCCTGGCTGGGGAAAGGCCCCGGCTTAAACCCGGGCCATGATGAAACCCTGAGCGGAAAGCTGCGGGGGGcccg                                                                                                >  2:380126/1‑237 (MQ=255)
                                                                                                                                   ccTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGTCTGGGGAAAGGCCCCGGATTAAAGCCGAGGCCTGGTGACAAGCTGGGCGGAATGCGgcgggcgggcgg                                                                                                >  2:371294/1‑235 (MQ=255)
                                                                                                                                      gcgcATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAAGCCCCGGATTAAAGCCGAGTCATGATGACACCCCGAACGGAATGCTGCTGGGGGCCTggga                                                                                             >  1:897431/1‑237 (MQ=255)
                                                                                                                                            ccTCAACCTATTGCAGCTGCAGGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATgg                                                                                       <  1:697277/239‑1 (MQ=255)
                                                                                                                                            ccTCAACCACTTGGAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATgg                                                                                       <  1:488704/239‑1 (MQ=255)
                                                                                                                                                 aCCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGCTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGGATTATGGCGCtt                                                                                  >  1:410714/1‑239 (MQ=255)
                                                                                                                                                       tGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCACCTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATg                                                                            <  1:784431/239‑1 (MQ=255)
                                                                                                                                                            cTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGCCACGCTGAGCGGAATGCTGCTGGCGGCCTGGGATTATTGCGCTTTAGATgcgcg                                                                        >  1:835721/1‑238 (MQ=255)
                                                                                                                                                                aTGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTTACGCCGAGTCATGCTGACACCCTCACGGGAATGTGGCCGGCGGCCCGGGTTTATGGCGCTTTAGTGGCGCGCCCCg                                                                   >  2:54996/1‑239 (MQ=255)
                                                                                                                                                                         aCGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACACCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGCATTATTGGGCTTTAGATGCGCGCTCCGGTGTTCCGt                                                          >  1:602025/1‑239 (MQ=255)
                                                                                                                                                                            cgcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTa                                                                                                                                                             <  2:486901/137‑1 (MQ=255)
                                                                                                                                                                            cgcgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTa                                                                                                                                                             >  1:486901/1‑137 (MQ=255)
                                                                                                                                                                              cgGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGCTTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGACTGGTATTATGGCGCTTTAGATGGGCGCCCCGGTCGGCCGgttttt                                                     >  2:410472/1‑239 (MQ=255)
                                                                                                                                                                                gaggaAACCGGGCTTTTCGGCCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTTCGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTg                                                    <  1:805877/238‑1 (MQ=255)
                                                                                                                                                                                 aggaAACCGGGTTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGc                                                   <  1:380126/238‑1 (MQ=255)
                                                                                                                                                                                         ggggTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGGGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGa                                          <  1:685415/239‑1 (MQ=255)
                                                                                                                                                                                         gggTTTTTCGGGCCCCTGGGGTTGGCGGCAAGCGTTCCGCTGCGTGCTGAGTTACGCCAGTTTGGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGCCCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGa                                          <  2:916114/239‑1 (MQ=255)
                                                                                                                                                                                           gCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCGCCGGTCGGCCGTTTTTTGCCTGTTCCGACa                                        >  1:263733/1‑239 (MQ=255)
                                                                                                                                                                                           gCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGGCATGATGACACGCTGAGCGGGATGCTGCTGGCGGCCTGGTATTATGGGGCTTTAGATGCCGGGCCCGGTCGCCCGTTTTTTGCCTGGTCCGACa                                        >  2:19802/1‑239 (MQ=255)
                                                                                                                                                                                            gttttGGGGCCCCTGGGGGTGGGGGCAAGGGATCCGGTGCCTGTTGAGTTACGGCAGTTTCGTCACTGGTTTCAGGCCGCGCTTAATGGCGTTAAGACCGAATGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACa                                        <  2:467708/237‑1 (MQ=255)
                                                                                                                                                                                            cTTTTCGGGCCACTGGCGTTGGCGGCAAGCGAGCCGCTGCCTGCTGAGGTACGCCAGTTTCGTCACTGTGTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGAGTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGTCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACa                                        <  1:692186/238‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCCGCTGGCGTTGGGGGGAAGCCATCCGCTGCCTGCTGAGTTTCGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGGCCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGc                                  <  1:547152/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGGCCCCTGGCTTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGc                                  <  1:715734/239‑1 (MQ=255)
                                                                                                                                                                                                   ggcccatGGCTTTGGGCGCAAGCGATCCGCTGCCTGCTTGTTTCCGCCAGTTTCGTCACGGTTTTCCGGCCGCGCTCAATGGGTTTAAGCACGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  2:556848/233‑1 (MQ=255)
                                                                                                                                                                                                    gcagcTGGCGTTGGCGGCAAGGGATCCGCTGCCTGCTGGGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:584693/234‑1 (MQ=255)
                                                                                                                                                                                                    gCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGGTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:876113/238‑1 (MQ=255)
                                                                                                                                                                                                    gCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAATTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:730093/238‑1 (MQ=255)
                                                                                                                                                                                                     ccACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGGGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGGTAAAGGCCCCGGATTAACGCCGCGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCaa                                <  1:87820/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCCATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTg              <  1:210500/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTg              <  1:19802/238‑1 (MQ=255)
                                                                                                                                                                                                                           ttCCGCTGGCTGCTGAGTTGCGCCAGTTTCGTCAGTGTTTTCAGGCCGCGGTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGt         <  1:797230/237‑1 (MQ=255)
                                                                                                                                                                                                                           ttCCGCTGCCGTCTGGGTTTCGCCCGTTTGGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc        <  2:765411/238‑1 (MQ=255)
                                                                                                                                                                                                                           aTCCGCTGCTTGCTGGGTTGGGCCAGTTTCGGCACTGTTTTCGGGGCGCGCTCAATGGCGTTAAGACCGACTGGCGGGACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGCCACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc        <  2:624210/239‑1 (MQ=255)
                                                                                                                                                                                                                            tCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTc        <  1:757587/238‑1 (MQ=255)
                                                                                                                                                                                                                             ccGCTGCGTGCTGGGGTTACCGAGTTGCGTCCCTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAg      <  2:719880/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tggctgctgAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGGTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTAt  <  1:371294/236‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
TGCGCAATTCGATGCGTTATGCGGTGGATTATGCGGCAATGAACGACCACAGGTTGTGGCTCAAGGGATTCGCCTCGGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCCTGGTATTATGGCGCTTTAGATGCGCGCTCCGGTCGTCCGTTTTTTGCCTGTTCCGACAATCTGCAACTCGTTACCACAGCGGTGAGCGTCAGTTAT  >  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: