Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,949 A→G 7.6% intergenic (‑207/+451) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,9490AG 7.6% 186.3 / 3.9 69intergenic (‑207/+451)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base G (4/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTC  >  NC_000913/4325662‑4326229
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                          
ttCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTCGGCACC‑At                                                                                                                                                                                                                                                                                           >  2:757150/1‑289 (MQ=255)
          aaCGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGCCGGCACC‑ATTTCCGGATgg                                                                                                                                                                                                                                                                                 >  1:756079/1‑288 (MQ=255)
            cgtacgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATgcgg                                                                                                                                                                                                                                                                               <  2:609294/289‑1 (MQ=255)
                cgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTa                                                                                                                                                                                                                                                                           <  2:756079/289‑1 (MQ=255)
                    gATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAAc                                                                                                                                                                                                                                                                        >  2:277/1‑288 (MQ=255)
                          aGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                 <  1:673292/289‑1 (MQ=255)
                           ggCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGGCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTGGTAGGGCTGATAAGACGCGCAAGCGGCGCAACAGGCAGTTGGCACC‑ATTGCCGGATGCGGCGGAAAGGCCTTAt                                                                                                                                                                                                                                                                 >  2:735352/1‑288 (MQ=255)
                            gCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTacc                                                                                                                                                                                                                                                                >  1:122613/1‑286 (MQ=255)
                               aTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                                                                                                                             >  2:600780/1‑288 (MQ=255)
                               aTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                                                                                                                             >  2:140400/1‑288 (MQ=255)
                                   cgTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                         <  1:778184/288‑1 (MQ=255)
                                    gTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCAGCCTACg                                                                                                                                                                                                                                                       <  2:11586/289‑1 (MQ=255)
                                               tACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                                                                                                                                                              <  1:140400/287‑1 (MQ=255)
                                                      caaagcggtgTAGGGACGTAACGGCAGTGTGGTTAGCTAAGTGAAGAGAGGCCTGTTAAGTGTAGTTACAAATTGATACAGATGAGGGGTGGCCTTTTG‑AGACGCGCGAGCGGCGCATCAGGCAGTCGGCACTGTGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                                                                                                                                                                      <  1:780138/279‑1 (MQ=14)
                                                        aaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                                                                                                                                                                   <  1:345371/289‑1 (MQ=255)
                                                        aaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                                                                                                                                                                   >  2:345371/1‑289 (MQ=255)
                                                        aaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTGCGAAAGTGAGACAGATGAGAGGCGGTCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGCATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                                                                                                                                                                   >  2:656706/1‑289 (MQ=255)
                                                              ggTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                                                                                                                                                             >  2:440958/1‑289 (MQ=255)
                                                                                   tgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCATC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGAAAGCGTCGCa                                                                                                                                                                                                         >  2:499542/1‑288 (MQ=255)
                                                                                       gACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGGCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                                                                                                                                                                    <  2:693609/289‑1 (MQ=255)
                                                                                       gACCTCAGTGAAGAGAAGCCAGTTAAGTGTAGTTGGAAATGTATACAGATGAGAGGCGGCCTGATA‑AGAAGCGCGAGCGTCGCATGAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                                                                                                                     <  1:14164/288‑1 (MQ=255)
                                                                                        aCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGCATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAgg                                                                                                                                                                                                    >  2:633638/1‑288 (MQ=255)
                                                                                         ccTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                                                                                                                                                                   <  2:710700/288‑1 (MQ=255)
                                                                                             aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                                                                                                                                                                                              >  2:214557/1‑289 (MQ=255)
                                                                                                    agaAGCCAGTTAAGTGTAGTGCGAAAATTATACAGATGGGGGGCGGCCTGATAGAGA‑GCGCGGGCGGCGCATCAGGCAGTGGGCACTGTTGCCGGATGCGTGGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑g                                                                                                                                                                                       <  1:727473/288‑1 (MQ=17)
                                                                                                       aGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGGGCGGGAATCTGTAGGCCTGATAAGACGCGCAAGCGTCGCACCCGGCAATCCGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTCCGGAGGGGCCGGGAACTTGTAGGCCTGCTAAGACGCGCACGCCTCCCATCGAGCAGTCGGCCCG‑Gttc                                                                                                                                                                                    >  1:18035/1‑287 (MQ=255)
                                                                                                       aGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                 <  1:96326/212‑1 (MQ=255)
                                                                                                       aGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                                                                                                                                 >  2:96326/1‑212 (MQ=255)
                                                                                                       aGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTg                                                                                                                                                                                    >  2:223848/1‑288 (MQ=255)
                                                                                                              tAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTGCGGGTGCGGGGATTTGTACGCCTGATAATCAGCGCAAGCGTCGCATCAGGCATCCTGCACC‑CGTGCCGGATgc                                                                                                                                                                            >  2:424507/1‑289 (MQ=14)
                                                                                                                       ttCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAAAGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGAGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGAAGCGGCGTAAAc                                                                                                                                                                   >  2:526029/1‑289 (MQ=255)
                                                                                                                                  taCAGCGGAGGGGGGGCCTGCGA‑AGACGGGCGGGGGTCGCATCAGGGAGTCGGCACTTTGGCCGGAGGCGGCGTAAACGCCTTTTCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                        <  1:757150/289‑1 (MQ=11)
                                                                                                                                      gATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCt                                                                                                                                                     >  2:228148/1‑288 (MQ=255)
                                                                                                                                      gATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAATTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGGCTTATCAGGCCt                                                                                                                                                     >  2:716983/1‑288 (MQ=255)
                                                                                                                                           gagGCGGCCGGAGA‑AGACGCGGGAGCGTGGCATCAGGCAGTCGGAACTGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                <  1:553688/288‑1 (MQ=255)
                                                                                                                                           agggcggcCTGATA‑CGAGCGGGGAGCGGCGCATCAGGAAGTCGGCAGTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                <  1:276234/286‑1 (MQ=11)
                                                                                                                                            agGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCCGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAATACGCGACAGCGTCGCCTCAGGCAGTCGGCACC‑TTTGCCGGCTGCGGCGTAAACGCCTTATCCGGCCTACGGTTg                                                                                                                                              >  2:641020/1‑289 (MQ=255)
                                                                                                                                            agGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCACC‑ATTGCCGGATTCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGTCGGGAATTTGTAGGCCTGCTAAGACGCGCAAGCGTCCCCTCAGGCAGTCGGCACC‑CTTGCCGGATGCGGCTGAAACGCCTTATCCGGCCTACGGCTg                                                                                                                                              >  2:781663/1‑289 (MQ=11)
                                                                                                                                                 gcCTGATT‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGACTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGATGGTGCg                                                                                                                                         >  2:10640/1‑289 (MQ=255)
                                                                                                                                                 gcCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                                                                           >  1:587564/1‑192 (MQ=34)
                                                                                                                                                 gcCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                                                                           <  2:587564/192‑1 (MQ=34)
                                                                                                                                                 gcCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGc                                                                                                                                                                 <  1:126263/265‑1 (MQ=37)
                                                                                                                                                 gcCTGATA‑AGACGCGCGAGCGTCGCACCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGTAGGGCGCCGGAATTTGTAGGCCTGACAAGACGCGCCAGCGCCGCACCAGGCAGGCGCACCC‑GTTGGCCCAGGCGGCGTAAACTCCCTACCCGGccgcagcagggcgcg                                                                                                                                         >  1:635267/1‑276 (MQ=11)
                                                                                                                                                    tGATA‑AGACGGGCTGGAGTGTCATCAGGAAGTCGGAACTGTTTCCGGATGGGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGa                                                                                                                                      <  1:214557/289‑1 (MQ=11)
                                                                                                                                                           agcgcgcgAGGGTTGCATGAGGGGATGGGCGCAGTTGCCGGGTGGGGTGAGAACGCCTTATGCGGCCTACGGAGGGGGCGGGAATTTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTGATCCGGCGTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                                                                <  1:633638/287‑1 (MQ=11)
                                                                                                                                                            aCGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                                       >  2:126263/1‑265 (MQ=37)
                                                                                                                                                               ggcgggtgtggcagacaggagttggaagtGTTGCCGGATGAGGCGGAAACGCCTTGTCCGGGCTAAGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                                                              <  1:174442/259‑1 (MQ=11)
                                                                                                                                                                                            gTTGCCGGAGGGGGAGTAAGCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                <  1:600780/288‑1 (MQ=255)
                                                                                                                                                                                                      gcggcgTAAAGTCCTTTTCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                      <  1:223848/288‑1 (MQ=37)
                                                                                                                                                                                                      gcggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                      <  2:390986/288‑1 (MQ=37)
                                                                                                                                                                                                                  cgttATCGGGCCGAGGGAGGGGGCGGGAATGTGTAGTCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGAACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGAt                                                                         <  1:277/286‑1 (MQ=255)
                                                                                                                                                                                                                    ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATACGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATTCGGCGGAAACGCCTTGTCCGGCCTCCGGACGGTCCGGGTCTTTGTAGGGCTGATAAGACCCGCCAGCGTCGCCTCCGGCACTTCGCA‑CTGTTGCCggtgtag                                                                      >  1:604627/1‑284 (MQ=11)
                                                                                                                                                                                                                    ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATgcg                                                                      <  1:440958/289‑1 (MQ=38)
                                                                                                                                                                                                                     tATCCGGCCTACGGTGGGGGCAGGAATGTGTAGGCCGGATAAAACGCGCAAGCGGCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCCGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATgcg                                                                      <  1:656706/288‑1 (MQ=255)
                                                                                                                                                                                                                          ggCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATGCGGCGTaa                                                                <  2:122613/289‑1 (MQ=38)
                                                                                                                                                                                                                                cGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                                                                                                                                                             >  1:314076/1‑144 (MQ=32)
                                                                                                                                                                                                                                aaccgatccggcaccaatGTAGGCCTTATACGACGCGCCAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATGCGGCGTAATCGCCTTATCCGGCCTAC‑‑‑GGTTCGGACTATTTGTAGGCACGATAAGACGCGACAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTTCG‑‑‑GACTATTTGTAGGCACGATAAGACGCGACAGCGTCGCATCAGGCATCATGCACCAATTGCCGGATGCGGCGTAAACGCCTTAt                                                       <  2:736947/272‑1 (MQ=17)
                                                                                                                                                                                                                                    ggggGCGGGAAATTGTGGGCCTGCTAAGACGCGCAAGCGTCGCATCAGGGAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATGCGGCGTAAACGCCTTAt                                                       <  1:415586/288‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                                                                                                     >  2:565512/1‑99 (MQ=33)
                                                                                                                                                                                                                                     ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                                                                                                     <  1:565512/99‑1 (MQ=33)
                                                                                                                                                                                                                                      gTGCGGGAATTTGTAGTCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATCTTGCGGCCTGATAAGACGCGCAAGCGTCGCATCACCCAGTCGGCCCC‑TTTGCCGGATGCGCCGTACACTCCTTATCCTG‑CTACGGGTGGCTCGGGACTTTGTAGTCCTGCTACGACCCTCCCGCCTTCCGTCAGGCATTTCGCA‑CTGTTGCTCGACGCGGCCTCAAccccctccccc                                                    >  1:35325/1‑281 (MQ=2)
                                                                                                                                                                                                                                          gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGCCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGACGGCGCGGGAATTTGTAGGCCTGCTAAGACGCCCCAGCGTCGCACCAGGCCGTCGGCA‑CCGTTGCCAGCCGCCGGGTAAACGCCTTGTcccccc                                                 >  1:231810/1‑284 (MQ=11)
                                                                                                                                                                                                                                             acttggaaGGGCTGATAAGACGCGCAAGGGGCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCGGAGGTTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                             <  1:716983/282‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                                                                                                                                                                      >  2:710368/1‑102 (MQ=32)
                                                                                                                                                                                                                                                 tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                                                                                                                                                                      <  1:710368/102‑1 (MQ=32)
                                                                                                                                                                                                                                                           aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGCATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACCCCCCAGCGTCGCATCCGGCCGTCGGCA‑CTGGTGCCGGTCTCGGCCTAAACGCCTTTTCCGGCCCACGGACGGTGGGGGAAt                               >  1:571644/1‑289 (MQ=11)
                                                                                                                                                                                                                                                                 cgcgCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGGGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAGAAGACGCGCAAGCGTCGGAGCAGGCAGTTGGCA‑CTGTTGCCGGAGGCGGCGTCATAGCCTGTTCCGCCTTACGGAGGGGGCGGGGATTTGTAg                         >  1:245399/1‑289 (MQ=11)
                                                                                                                                                                                                                                                                             gCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                                                                           >  2:738674/1‑69 (MQ=35)
                                                                                                                                                                                                                                                                                        tCGGCACT‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCACCAGGCACTTGGCA‑CTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCGGGTGGGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc  >  1:319897/1‑289 (MQ=16)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                          
TTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATA‑AGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑ATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACC‑GTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA‑CTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTC  >  NC_000913/4325662‑4326229

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑