Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,094,457 C→T intergenic (‑185/+266) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,4570CT94.5% 148.7 / 1.3 56intergenic (‑185/+266)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (30/22);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

CAGGCGTAGCGTCAGGTTGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAA  >  NC_000913/4094162‑4094742
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                             
cAGGCGTAGCGTCAGGTTGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGAATGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACg                                                                                                                                                                                                                                                                                           <  1:156273/300‑1 (MQ=255)
               gTTGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCATGCGTCGGTGGTGGCTTTNATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATANGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATc                                                                                                                                                                                                                                                                            <  2:656241/300‑1 (MQ=255)
                    cgccgTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                       <  2:20489/300‑1 (MQ=255)
                      acgTTGCGCTCATCCCAGCCTTTCAGCCAGGNGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGNACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCCGGCCTa                                                                                                                                                                                                                                                                     >  2:491072/2‑300 (MQ=255)
                          tGCGCTCATCCTAGCCTTTCATCCATGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                                                                                                                                                                                                                                 <  2:310316/300‑1 (MQ=255)
                          tGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGTTGGTGGCTTTGATCATTCCCTGGACAAACCAGGATTGATTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTAGGTGTTGTCGGATGCACGACCCGTAGGCTGGATAATGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTATGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                                                                                                                                                                                                                                 <  1:355111/300‑1 (MQ=25)
                             gcTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGt                                                                                                                                                                                                                                                              >  1:212258/1‑300 (MQ=255)
                                  tCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                                                                                                                                         >  2:402724/1‑300 (MQ=255)
                                  tCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATCCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                                                                                                                                         >  1:638413/1‑300 (MQ=255)
                                            tCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGc                                                                                                                                                                                                                                               <  1:111453/300‑1 (MQ=255)
                                                   ggCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCCGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCAGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAAGCCTACAACCGTCGCAGACCCGCAGGCAGGATaa                                                                                                                                                                                                                                        >  2:283450/1‑300 (MQ=21)
                                                     cGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCCGATGTGACGCTGTCGCGTCTTATCAGGCCTCCAACTGTTGCCGACCTGTCGGCCGGATacgg                                                                                                                                                                                                                                      >  2:536161/1‑297 (MQ=255)
                                                           tggtggCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCTGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtg                                                                                                                                                                                                                                <  1:557545/300‑1 (MQ=255)
                                                           tggtggCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtg                                                                                                                                                                                                                                <  1:380776/300‑1 (MQ=255)
                                                           tagctacTTTGATCGCTCCATGGACAAACAAGGACTTGTTCATGTTTTGCATTTTCTGTATTCCGATAAATTATGTGTGTTCGGATGCACGACCCGTAGTCCGGGTAAGTCACTCGCGCCGCATCCGGCAATCAATGCCTGATACGACACGGTCGCTCCTTATCAGGCCGACAACTATTGCCGACCTTTAGGTCCAATAAGGCGCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGTTGTCGCGTCTTATCATGCCTACAACTGTTGCCGACCTTTAGGCCGGATAAGGCACTtg                                                                                                                                                                                                                                <  2:95204/294‑1 (MQ=9)
                                                           gggcggctgttaTCATTCCCTGGACAAGACGGGAATGAGTAATGTTGTGCATGTGCTGTGTTCCTGTAAATTCGGTAGTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCTTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACACTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtg                                                                                                                                                                                                                                <  2:95273/289‑1 (MQ=18)
                                                                                     aaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTACATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAACCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTgc                                                                                                                                                                                                      >  1:164917/1‑299 (MQ=255)
                                                                                     aaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTCCAACTGTTGCCGACCTGTAGGCCGGATACGGCACTTGTGCCGCATCCGGCAATCAATGCCTGa                                                                                                                                                                                                      >  1:517670/1‑300 (MQ=255)
                                                                                     aaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGa                                                                                                                                                                                                      >  1:278371/1‑300 (MQ=255)
                                                                                     aaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCAGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATAACTGCCTGa                                                                                                                                                                                                      >  1:404438/1‑300 (MQ=255)
                                                                                     aaCCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAACTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAACCAATGCCTGa                                                                                                                                                                                                      >  2:284683/1‑300 (MQ=255)
                                                                                     aaCCAGGACTGAGTAATGTCTTGCATATTCTGTGTTCCTGTAAATCCGGTGTTGTCGGATGCACGACCCGTAGGTCGGCTAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCATAAAACTATTGCCGACCTGTAGGCCTGCCAAGGCACTTGTGCCGCATCCGGCAATTAATGCATGATGTGACGCTGTCGCGTCCTATCAGGCCTACAACCGTTGCCGACCTGTAAGCTGGATAAGGCACCTGCGCCACATCAGGcaagcaatacccgc                                                                                                                                                                                                      >  2:218535/1‑294 (MQ=18)
                                                                                        cAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCNTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                                                                                                                                                                                                   <  1:467152/300‑1 (MQ=255)
                                                                                                  tAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACNCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcg                                                                                                                                                                                         <  1:506607/300‑1 (MQ=255)
                                                                                                      gTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTATAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTTGTGCCGCATCCGGcaattcttacctgctgcgacgccacctctct                                                                                                                                                                                     >  1:300229/1‑273 (MQ=17)
                                                                                                                ttccgtgtTCCTGTAAATTCCGTGTTGTCAGATCCACGACCCGTATGCCGGATAAGGCGCTCGCTCCACATCCGGCAAGCAATGCCTGATGCGACGCTGTCGCTTTTCATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGACACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                                                                                                                           <  2:270705/296‑1 (MQ=17)
                                                                                                                     tgtTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCTTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACaa                                                                                                                                                                      <  1:358547/300‑1 (MQ=255)
                                                                                                                        tCCTGTAAATTGTTTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTTGCGCCGCATCCGGCAATCAATGCCTGATGCGACACTGTCGCGTCTGATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCGGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATTCGACGCTGTCGCGTCTTATCAGGCCTACAACTa                                                                                                                                                                   <  2:360601/300‑1 (MQ=14)
                                                                                                                         ccTGTAAATTCGGTGTTGTCGGATGCACGACCCTTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTAt                                                                                                                                                                  <  1:704782/300‑1 (MQ=255)
                                                                                                                                     gtgtTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAACTGTa                                                                                                                                                      >  1:181454/1‑300 (MQ=255)
                                                                                                                                     gtgtTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAAAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCATGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGAATCTGGCAATTAATGCCTGACGCGACGCTGTCGCGTCTTATCAGGCCTAAAACTATTGCCGACCTGTa                                                                                                                                                      >  2:146242/1‑300 (MQ=17)
                                                                                                                                      cgttgtCGGATGCATGACCTGTAGGCCGGATGATGCGTGCGTGCAGCATCCGGTAATAGAGGTCGGATGCGACGTTGTCGCGTCGTATCAGGCCTACAACTATTGCCGATCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAAGGCCTGATGGGACGCTGTCGCGTCGCATCAGGCCTACAACTGTTGTCGACCTGTAGGCCGGATAAGGCACGTGTGCCGCATCTGGCAAGCAGTGCCTGATGCGACGCTGTCGCGTCTGATCAGGCTTACAACTATTGCCGACCTGTAg                                                                                                                                                     <  2:109206/299‑1 (MQ=255)
                                                                                                                                                  cACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGTACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGc                                                                                                                                         <  1:613210/300‑1 (MQ=255)
                                                                                                                                                       cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCGCTGCCCGATGCGACGCTGTCGCGTCTTCTCAGGCCTACAACTATTACCGAACTGTCGGCCTGATAAGGCACTTGTGCAGCATCTGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAAGACTACAACCGTTGCCAACCTGTAGGCAGGATACGGCACTCATGCCGCATCCGGCACTCAATGCCTGACGCTACGCTGTCGCATCTTATCAGGCCTCCAACCATTGCCGACATGTAgacaagataaggaacctg                                                                                                                                    >  1:609602/1‑283 (MQ=255)
                                                                                                                                                       cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCCGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTCCAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCCGGCCTACAACTATTGCAGACCTGTAGGCCGGCTAAGGCACTtg                                                                                                                                    >  1:609432/1‑300 (MQ=255)
                                                                                                                                                                 ggATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATcc                                                                                                                          >  1:707708/1‑300 (MQ=255)
                                                                                                                                                                 ggATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGAAGCTGTCGCGTCTTATCAGGCCTACGACTATTGCCGACCTGTAGGCCTGCTAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTATGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTAGCGTCTTATCGGGCCTACAACTATTGCCGACCTGTAGGCCGGATACGGTACTTGTGCCGCATcc                                                                                                                          >  2:356844/1‑300 (MQ=14)
                                                                                                                                                                                             atcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAACTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc                                                                                              >  1:256338/1‑300 (MQ=255)
                                                                                                                                                                                                   gCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAAGCTGTCGCGTCTTCTCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                        >  1:600864/1‑300 (MQ=255)
                                                                                                                                                                                                       gATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                    <  1:670817/300‑1 (MQ=255)
                                                                                                                                                                                                                  ttcgcgtctcacagtgcctacgccggttGCCGACCTGTATGCTGGATAGGGCCCGTGTGCCGCATCCGGCAATCAAGGTCTGATGTGACGCTGTCGCTTCTTATCAGGCCTACAACGGTTCCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                                         <  2:592114/274‑1 (MQ=11)
                                                                                                                                                                                                                  gTCGCTTCTTATCAGGCCTACACCTATTGCTGGCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGCTTTTACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                                         <  2:599756/300‑1 (MQ=17)
                                                                                                                                                                                                                  gTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTTTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                                         <  1:533526/300‑1 (MQ=255)
                                                                                                                                                                                                                                ggCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGt                                                           >  1:713477/1‑300 (MQ=255)
                                                                                                                                                                                                                                ggCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCAGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTCTCCGGCATACGACTGTTGCCGACCCGt                                                           >  2:78718/1‑300 (MQ=255)
                                                                                                                                                                                                                                 gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCCCTGCCTGATGCGACGCTGTCGCGTCTTATCGGGCATACAACTATTGCCGAACTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAACCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAATGTTgtaaacacgta                                                          >  2:58257/1‑290 (MQ=17)
                                                                                                                                                                                                                                 gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTa                                                          >  1:185919/1‑300 (MQ=255)
                                                                                                                                                                                                                                  ccTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAg                                                         <  1:313162/300‑1 (MQ=255)
                                                                                                                                                                                                                                                    gtGTAGGCATGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAACTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCCGATGCGAAGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCATGATGCGACGCTATCGCGTCTTATCAGGACTACAACTGTCGCCGACACGTAGGACGGCTAGGGCGCTcgc                                       >  2:450848/2‑300 (MQ=14)
                                                                                                                                                                                                                                                    cTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCACTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGCCTTATCAGGCCTACAACTGTTGCCGAACCGTAGGCCGGATACGGCGCTcgc                                       >  1:653979/1‑300 (MQ=255)
                                                                                                                                                                                                                                                    cTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGNCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTcgc                                       >  1:488228/1‑300 (MQ=255)
                                                                                                                                                                                                                                                            cTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATTTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATc                               <  1:26164/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                tAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                           <  1:304928/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cTTGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGttt                    >  2:104469/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                        ttGTGCCGCATCCGGCAATCAATGCCTGATGTGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTa                   >  1:251827/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         gtcaatGCCTGATGTGACGCTGTCGCGTCTTATCGGGCCTAAAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGTTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCAGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCaaaa  >  1:347827/2‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                             
CAGGCGTAGCGTCAGGTTGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTCCCTGGACAAACCAGGACTGAGTAATGTTTTGCATGTTCTGTGTTCCTGTAAATTCGGTGTTGTCGGATGCACGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAA  >  NC_000913/4094162‑4094742

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