Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,204,468 C→G 5.8% intergenic (+179/‑128) yehK → / → yehL uncharacterized protein/putative hexameric AAA+ MoxR family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,204,4680CG 5.8% 223.4 / 3.7 69intergenic (+179/‑128)yehK/yehLuncharacterized protein/putative hexameric AAA+ MoxR family ATPase
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base G (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCTTGGCGACGAAGCGCAAAATATCAGCCGTA  >  NC_000913/2204169‑2204767
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
gTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTcc                                                                                                                                                                                                                                                                                                             >  1:405500/1‑300 (MQ=255)
        atataAGGGATGGGTTTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATACAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAgg                                                                                                                                                                                                                                                                                                    <  1:19271/301‑1 (MQ=255)
           taAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGttt                                                                                                                                                                                                                                                                                                 >  2:226670/1‑301 (MQ=255)
              gagaTGAGTGTATTGATTAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGtttttt                                                                                                                                                                                                                                                                                              <  1:230519/301‑1 (MQ=255)
                     gtgtATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAAAACCTGCTGGTGAGGCAGGTTCTtttatttat                                                                                                                                                                                                                                                                                       >  1:298272/1‑300 (MQ=255)
                        tATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAAAACCTGCTGGTGAGGCAGGTTCTTTTATTTACtat                                                                                                                                                                                                                                                                                    <  2:298676/301‑3 (MQ=255)
                           tgatgaGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGAt                                                                                                                                                                                                                                                                                 <  2:460588/301‑1 (MQ=255)
                            gatgaGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATc                                                                                                                                                                                                                                                                                >  2:398806/1‑301 (MQ=255)
                               gaGCGGAAATTACGAGGTGAATCACCGGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACGCTAAATAATTCAAATTGGTCCGAGCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCGTAATTGCAACGGGGTCCAGCCGTGGGTATAATTCGGTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGc                                                                                                                                                                                                                                                                             <  1:481012/301‑1 (MQ=255)
                                    gAAATTACGAGGTTAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTcc                                                                                                                                                                                                                                                                        <  2:406063/301‑1 (MQ=255)
                                     aaaTTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCa                                                                                                                                                                                                                                                                       <  1:226427/301‑1 (MQ=255)
                                          aCGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGAAGGTTTTTTTATTTCCTGTGATCTGCTTTCCAAatat                                                                                                                                                                                                                                                                  >  1:237673/1‑301 (MQ=255)
                                                         tGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGACTCACCAGCAGGTTTTTTTATTTACTGTGATATGCTTTCCAGATATTTTTCGCTCAaacaa                                                                                                                                                                                                                                                   >  2:129429/1‑301 (MQ=255)
                                                               tGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAAtt                                                                                                                                                                                                                                             >  2:540217/1‑300 (MQ=255)
                                                                            tttAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGCTCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACAttat                                                                                                                                                                                                                                >  1:601276/1‑299 (MQ=255)
                                                                             ttAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATATCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTaa                                                                                                                                                                                                                               >  2:339763/1‑300 (MQ=255)
                                                                                        attcAGTATTTATTCAACCAAGAAATAGAATATTATACTCTAAATATTGTAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAAGTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTaaaa                                                                                                                                                                                                                    <  1:554952/298‑1 (MQ=255)
                                                                                        atACAGAATATTTTAAACCAAGAAATAGAAGAATATACTCTAAATAATTCAAAGTGGTCCGATCCGGCGAAACGTCCCAATGGCCGGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTGTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTaaaa                                                                                                                                                                                                                    <  1:554945/300‑1 (MQ=255)
                                                                                         tACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAat                                                                                                                                                                                                                   >  1:130303/1‑301 (MQ=255)
                                                                                          aCAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAata                                                                                                                                                                                                                  >  2:477544/1‑301 (MQ=255)
                                                                                             gAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCg                                                                                                                                                                                                               >  2:285391/1‑301 (MQ=255)
                                                                                                  aatttatcctaGAAAGAGAATAATATACGCTAAATAATTTAAATTGGTCCGGTCGTGCGCAACGTCCCAAGGGCCTGGGTTATAAATCACATTAACTTAATTGCAACGGGGTCTAGCCGTGGTTGTAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGGCTGAAAATAGTAAATAAAAGAGCCTGCCTCACCAGCAGGTTTTTTTATATACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAGCATTTATTGCGCGTAAAATATCGTTTAt                                                                                                                                                                                                          <  1:511818/291‑1 (MQ=255)
                                                                                                    ttCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTAttt                                                                                                                                                                                                        <  1:285088/301‑1 (MQ=255)
                                                                                                      cAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTATATTTCa                                                                                                                                                                                                      >  2:590083/1‑301 (MQ=255)
                                                                                                        aCCATGAAATAGAATAATATACTCTAAATAAGTCAAATTGGTCCGATCCGGTGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAAGTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCAtt                                                                                                                                                                                                    <  1:25797/301‑1 (MQ=255)
                                                                                                        aCCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCAtt                                                                                                                                                                                                    <  2:130427/301‑1 (MQ=255)
                                                                                                             gAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATaa                                                                                                                                                                                               >  2:627046/1‑300 (MQ=255)
                                                                                                                  aGAATAATTTACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTc                                                                                                                                                                                          <  1:626169/301‑1 (MQ=255)
                                                                                                                       aataTAATCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGa                                                                                                                                                                                     <  1:539544/301‑1 (MQ=255)
                                                                                                                         atgaCTCTAATTGATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATg                                                                                                                                                                                   <  1:126630/298‑1 (MQ=255)
                                                                                                                                 aaaTAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATACAGTTTATTTCATTAAAACAtttacagaataaaaaaatt                                                                                                                                                                           >  1:97539/1‑285 (MQ=255)
                                                                                                                                 aaaTAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCACCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAatatatt                                                                                                                                                                           >  1:517046/1‑300 (MQ=255)
                                                                                                                                       ttCAAATTGGGCCGATCCTGCGTAACGTCCTAATGTCCTGGATTATAAATGTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTGTTTATTTACTGTGATCTGCTGTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCg                                                                                                                                                                     <  1:390127/301‑1 (MQ=255)
                                                                                                                                         cAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCa                                                                                                                                                                   >  2:326714/1‑301 (MQ=255)
                                                                                                                                               ggTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAACCATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCAACGCAGaaaaa                                                                                                                                                             >  1:350207/1‑301 (MQ=255)
                                                                                                                                                 tCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAgatcggaagagcacacgtctgaactccaga                                                                                                                                                           >  2:274388/1‑272 (MQ=255)
                                                                                                                                                  ccGTTCCGGGGCAACGTCCCAATGGCTTTGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAatca                                                                                                                                                          <  1:589322/301‑1 (MQ=255)
                                                                                                                                                          gcgcAACGTCCCAATGGCCTTGATTATCAATCTCATTATCTGAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAAAAAGATGTTGCAAAGGCTGTGGTGTTAGGAAGGAAAATAGTAAATAAAAGAACCTGCCTAAAAAGCAGGTTGTGTGATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGc                                                                                                                                                  <  1:129305/301‑1 (MQ=255)
                                                                                                                                                                 ttCCCAATGGCCTGGATTATTAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGcc                                                                                                                                           <  2:237940/300‑1 (MQ=255)
                                                                                                                                                                 gTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGcc                                                                                                                                           <  2:517718/301‑1 (MQ=255)
                                                                                                                                                                    ccAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCctgc                                                                                                                                        >  2:482239/1‑301 (MQ=255)
                                                                                                                                                                         tgcttGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCctgctgctg                                                                                                                                   <  2:49363/297‑1 (MQ=255)
                                                                                                                                                                         ggCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCAACTACAACGTCCGCcagctgctg                                                                                                                                   >  1:49295/1‑301 (MQ=255)
                                                                                                                                                                              ggtgtattaatttgtttttCTTAATTGTAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTta                                                                                                                              <  1:231674/283‑1 (MQ=255)
                                                                                                                                                                               gATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTtat                                                                                                                             <  2:388780/301‑1 (MQ=255)
                                                                                                                                                                                    taAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCAGCTGTGTTATACGcc                                                                                                                        >  2:451581/1‑301 (MQ=255)
                                                                                                                                                                                                     aTTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAACATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATACTCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAAtt                                                                                                       >  1:46756/1‑301 (MQ=255)
                                                                                                                                                                                                         aaaCGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTaaaa                                                                                                   <  1:398292/300‑1 (MQ=255)
                                                                                                                                                                                                              gggTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACaaaa                                                                                              <  2:46809/301‑1 (MQ=255)
                                                                                                                                                                                                                   tggagttcagacgtgtgctcttccgatctcacGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATa                                                                                         <  2:11114/272‑1 (MQ=255)
                                                                                                                                                                                                                    aGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTTTTATACGCCGAAGAACTGGCAAAATAAAAACAAAATGaaaa                                                                                        >  1:370765/1‑301 (MQ=255)
                                                                                                                                                                                                                              ttaATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTgc                                                                              <  2:371261/300‑1 (MQ=255)
                                                                                                                                                                                                                               taATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTgcc                                                                             >  1:422845/1‑301 (MQ=255)
                                                                                                                                                                                                                                      ttcgttGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAACTTAAAACAAAATGATAACGCACCTTGCCCGCCCgg                                                                      >  1:99840/1‑301 (MQ=255)
                                                                                                                                                                                                                                        tgttGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGtt                                                                    <  2:97663/300‑1 (MQ=255)
                                                                                                                                                                                                                                            ggAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCa                                                                <  2:350670/301‑1 (MQ=255)
                                                                                                                                                                                                                                             aaaaaggATGGTGAAAAGGTTGTGGTGTTAGGCCTGAAAATATGAAATTAAAGAACTGGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATGTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAg                                                               <  1:326310/295‑1 (MQ=255)
                                                                                                                                                                                                                                                cacGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTACGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAAAAAAATGATaaaatcggaagagcgtcgtgaagggaaaga                                                            >  1:11114/1‑273 (MQ=255)
                                                                                                                                                                                                                                                  cGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAg                                                          >  1:310743/1‑301 (MQ=255)
                                                                                                                                                                                                                                                    aTGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCACGCCCGGTTCGCAGTTaaaat                                                        >  1:343707/1‑298 (MQ=255)
                                                                                                                                                                                                                                                            aaGGCTGTGGTGTTAGGCTTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCg                                                <  1:451017/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttttccgttctacaGTAAATAAAAAAACCTGCTGGTTAGTCAGGTTCTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTa                                   <  1:9474/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCt                               <  2:292012/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    acaGTAAATAAAAAAACCTGCTGGTGAGGCAGGTTCTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTaagatcggaaga                        >  2:9474/3‑290 (MQ=255)
                                                                                                                                                                                                                                                                                        aaaataaaAGAACTTTCCTCACCAGCAGGTTTTTGTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAAGGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGAGGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCTTGGCGACGAAg                    <  1:49468/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         aaataacAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTACTCCGCCAAACATTTATTCCGCCTCAAATCTCGTTTATTTCTTTAATACATTTCAGGCATGAAAATCTGTCACCGCAGAATAATCATCGGCAGCTTCCGCCTGATGACGTGTTATACGCCGATGAACTGGCAAAAGTACAACAAACATATAACGCACCTTTCCCACCCGGTTGGCAGTTACATTTGCATGCGGCCCGTGCTTTTATCCTTGGCGACGaacc                   >  1:569244/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                              aaaGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCATTTAAACTTGGCGACGAAGCACaaa              >  2:608018/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                     cTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCTTGGCGACGAAGCGCAAAATATCAg       <  2:423408/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         cTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCTTGGCGACGAAGCACGAAATATCAGCCGt   >  1:432946/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                          tCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCTTGGCGACGAAGCGCAAAATATCAGCCGTa  >  2:475306/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
GTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTGCTGTGTTATACGCCGATGAACTGGCAAAATTAAAACAAAATGATAACGCACCTTGCCCGCCCGGTTGGCAGTTAAGTTTGCCTGCGGCCCGTGCTTTTATCCTTGGCGACGAAGCGCAAAATATCAGCCGTA  >  NC_000913/2204169‑2204767

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: