Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,820 C→T G129R (GGA→AGA)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,8200CT97.9% 269.1 / NA 98G129R (GGA→AGA) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (42/51);  minor base G (0/2);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCT  >  NC_000913/4193576‑4194067
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
ccGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTTCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGGTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAgcgc                                                                                                                                                                                                                                                     <  2:217485/250‑1 (MQ=255)
  gccaaccacTGCAATTTCCAGCAACTGAGCCGTTGCCAAACGCAAGTCCGCGGCTTGAGTAATGGCGCTGAAGACGGCGATACTGCCGACACCCGT‑TGCTATGACCGCAGGCGCGCGTGACAGACTGATGCCGCCAATCGCCACGGGGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCAGTTCCAGCCCCTGCGGTGAAGACGGCCTCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGcgcg                                                                                                                                                                                                                                                    >  2:53092/9‑250 (MQ=255)
     aaCTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGatat                                                                                                                                                                                                                                                >  1:576890/1‑251 (MQ=255)
         ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCCATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGCACACATGTCTCAGCGCGAtata                                                                                                                                                                                                                                               >  1:100145/1‑248 (MQ=255)
         ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAag                                                                                                                                                                                                                                             <  1:472114/250‑1 (MQ=255)
            gCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAgg                                                                                                                                                                                                                                          <  1:291805/250‑1 (MQ=255)
             cAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAgggg                                                                                                                                                                                                                                        <  1:243068/251‑1 (MQ=255)
              aaTTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGc                                                                                                                                                                                                                                       >  1:473500/1‑251 (MQ=255)
               aTTTCCAGCAACTGTGCCGTTGCCAATCGCAAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCg                                                                                                                                                                                                                                      <  2:505905/251‑1 (MQ=255)
                    cAGCAACTGTGCCGTTGCCCAACCCCAGTCTGCGGCTTGAGTCATGGCGCTGACGCCGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCCCGGTGGGATCATCCGCCCGTCGCTCAACATGCCGTGCCCGCTGTTCCAGCCCCTGCGGTGCAGAATGCATCTGTTTGGTTTGCGTCGGCAACACCTGTCTCAgcgctatctaaccgggtctttcct                                                                                                                                                                                                                                 >  2:135261/1‑231 (MQ=255)
                    cAGCAACTGTGCCGTTGCCAAACGCCAGTCGGCGGCTTGAGTAAGGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGtgctg                                                                                                                                                                                                                                 <  2:338682/251‑1 (MQ=255)
                       cAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGccc                                                                                                                                                                                                                              >  2:129082/1‑250 (MQ=255)
                         aCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAgc                                                                                                                                                                                                                            <  1:453009/251‑1 (MQ=255)
                            gtgCCGTTGCCAAACGCCAGTCTGCGGCTTGGGTAATGGCGCTGACGACGGCGCTCCGGCCGACACCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTGCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCTCGATATAAGAGGGGCGTGCTGCCAgcgc                                                                                                                                                                                                                          <  2:589360/250‑1 (MQ=255)
                             tgCCGTTGCCAAACGCCAGTCTGCGGGTTGCGTATTGGCACTGACGCCGGCGATACTGCAGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATACTGCGCCAGTCGCTCAACATGCTGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCATCGCGCTATACGAGGGGCGTGCTGCCAGcgcgc                                                                                                                                                                                                                        >  2:64253/1‑250 (MQ=255)
                                    gCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAt                                                                                                                                                                                                                 <  2:16044/251‑1 (MQ=255)
                                     ccaagctccAGTCTGCGGCTTGAGTAATGGCGCTGCCGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAtt                                                                                                                                                                                                                <  1:190041/244‑1 (MQ=255)
                                     ccAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAtt                                                                                                                                                                                                                <  1:460487/251‑1 (MQ=255)
                                     ccAAACGACACTCTGCGGCTTGAGTAATGGCGCTGACGACGCCGATACTGCCGCAACCAGT‑TGCTATCACCGCAGGCGAGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATACTCCGCCCGCCGCTCAACATGCCGAGCCAGCTGTTCCAGACCCTGCGGTGCAGCAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCCGCGCGGTATAAGAGGGGCGTGCTGCCCGCGCGACGCCGAtt                                                                                                                                                                                                                >  2:378499/1‑251 (MQ=255)
                                      caaccGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                                                                                                                                               <  2:481917/247‑1 (MQ=255)
                                      cAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTACCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                                                                                                                                               >  2:263377/1‑251 (MQ=255)
                                       aaaCGCCAGGCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGCCACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTc                                                                                                                                                                                                              >  1:425398/1‑251 (MQ=255)
                                          cGCCAGTCAGCGGGTTGAGTAAGGGCGCTGACGACGGCTATTCTGCCGACACCCGT‑TGCTAACACCGCAGGCGCGTGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGAGAATCCGCCAGGCGCTAAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCAACTGTTAGGTTTGCGTCGGGAACACATGTCTCAGTGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCat                                                                                                                                                                                                           <  1:44708/251‑1 (MQ=255)
                                           gCCAGTTTGCGGCTTAAGTAATGGTGCTGACGTTGGCTATACGGTCGCCACCCGGGTG‑TATCACCGCCGGCGCCGGTTCCAGACTGATGCCGCCAATCTCCACGGTTGGATAATCCTCCAGTCGCTCAACATGCCGCGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGTATCTGTGTGGTTTGCGTCTGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGGCGAtgaccata                                                                                                                                                                                                          <  2:541795/251‑8 (MQ=255)
                                               gTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGCCACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑Gt                                                                                                                                                                                                     >  1:210875/1‑251 (MQ=255)
                                                tCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAAAGGGGCGAGCTGCCAGCGCGACGTCGATTTCCCTATC‑GTc                                                                                                                                                                                                    >  1:401953/1‑251 (MQ=255)
                                                 cTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTc                                                                                                                                                                                                    >  2:569788/1‑250 (MQ=255)
                                                  ggCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCAt                                                                                                                                                                                                  <  1:401030/250‑1 (MQ=255)
                                                   gCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCAtg                                                                                                                                                                                                 <  1:368248/251‑1 (MQ=255)
                                                    cGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCAtgt                                                                                                                                                                                                >  1:213780/1‑251 (MQ=255)
                                                        ttGAGTAAGGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTc                                                                                                                                                                                             <  2:440716/250‑1 (MQ=255)
                                                                    gctgaagccgGCGATACTGCCGACACCCAT‑TGCTATCACCGCAGGCGCGCGAGCCAGCCTGCTGCCGCCAATCGCCACGGGGGGACAATACGCCGGTCGCCCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGCCGTCTAGTCTCATATC‑GTCATATGTCGAAACGCCCCGcc                                                                                                                                                                                >  2:570979/9‑251 (MQ=255)
                                                                            cgGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GGCATGTGTCGAAACGCCCAGCCGCAGGCCt                                                                                                                                                                        <  1:483063/251‑1 (MQ=255)
                                                                             ggCGATACTGCCGTCACCCGT‑GGCGATCACCGCGGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGGGTCGAAACGCCCAGCCGCATGCCTg                                                                                                                                                                       <  2:583866/251‑1 (MQ=255)
                                                                               cGATCCTGCCGACCCCCGT‑TGCGATCCCCGCAGGCGCGCGTGCCAGACGGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGcc                                                                                                                                                                     <  2:165818/251‑1 (MQ=255)
                                                                                aaTACTGCCGACACCCGT‑TGATATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGGCGCTCAACATGCCGTGCCAGATGTTCCAGCCCCTGCGGTGCAGGAGCCATCGGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGTGGCGTGCTGCCAGCGCGACGTCGCTTTCCCTATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGccc                                                                                                                                                                    >  2:481798/2‑250 (MQ=255)
                                                                                    cTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGAGCCAGCTGTTCCAGACCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCcgcgcg                                                                                                                                                                >  1:212143/1‑251 (MQ=255)
                                                                                         gACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAACAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATC‑GTCATGTGTCGAACCGCCCAGCCGCAGGCATGCCCCGCGGCTgg                                                                                                                                                           >  2:486090/1‑251 (MQ=255)
                                                                                              ccGT‑TGCGATCACCGCAGGCGGGCGTGCCAGACTGATGCCGACAAGCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTg                                                                                                                                                      <  1:331313/251‑1 (MQ=255)
                                                                                                gT‑TGCGATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCCCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTgag                                                                                                                                                    <  2:461339/251‑1 (MQ=255)
                                                                                                    gCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                                                                                                                                                 <  1:430389/251‑1 (MQ=255)
                                                                                                    gCTATCACCGCAGGCGCTCTTGCCAGACGGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCCTAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGGCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCTTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                                                                                                                                                 <  2:187380/251‑1 (MQ=255)
                                                                                                    gCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                                                                                                                                                 <  1:243289/251‑1 (MQ=255)
                                                                                                        tCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTg                                                                                                                                             >  2:137907/1‑251 (MQ=255)
                                                                                                          cccGCAGGCGCGCGTGCCAGACTGTTGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGc                                                                                                                                           <  1:559215/250‑1 (MQ=255)
                                                                                                          aCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAACCTCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGc                                                                                                                                           >  1:71286/1‑251 (MQ=255)
                                                                                                             gCAGGCGCGCGTGCCAGACTGATGCAGCCAATCCCCACGGTGGGATAATCCACCAGTCGCTCAACATGCCGTGCCAGCTGTTCCCGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGACCACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGCGATCGGTGGCTTg                                                                                                                                        >  2:574211/1‑251 (MQ=255)
                                                                                                              cAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGc                                                                                                                                       >  1:556188/1‑251 (MQ=255)
                                                                                                               aGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGCTAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGc                                                                                                                                       >  1:451203/1‑250 (MQ=255)
                                                                                                                   gctcgTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATc                                                                                                                                  >  1:461831/4‑251 (MQ=255)
                                                                                                                   gcgcgTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATc                                                                                                                                  <  2:140949/251‑1 (MQ=255)
                                                                                                                   gcgcgTGCCAGACTGATGCCGACAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GCAATGGGTCGAAACGCCCAGCCGCAGGCCTGCCGAGCGGGTGGCATTGAGATCGGTGGCTTGCaactc                                                                                                                                  >  2:193453/1‑248 (MQ=255)
                                                                                                                      ggtcccaggcagtgtccGCCAATCGCCTCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTc                                                                                                                               <  2:175292/236‑1 (MQ=255)
                                                                                                                         gCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑Tg                                                                                                                           <  1:111914/251‑1 (MQ=255)
                                                                                                                           cAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGcc                                                                                                                         <  1:137907/251‑1 (MQ=255)
                                                                                                                             gACTGATGCCGACAATCGCCGCGGGGGGATAATCCGCCAGGCGCTCAAAATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGAAATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGACCAGCCGCAGGCCTGCAGAGCGGATGGCATTGCGATCGGTGGCTTGCAAATCTTCC‑TGcccc                                                                                                                       >  2:46950/1‑251 (MQ=255)
                                                                                                                                     cctccAATCGCCACGGTGGGATAAACCGCCAGTCGCTCAACAAGCCGTGCCAGCTGTTCCAGCCCCTGCGGGGCAGAAGGCATCTGTTTGGTTTGCGTAGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑CGCCCCAAATGGAc                                                                                                               <  2:35290/248‑1 (MQ=255)
                                                                                                                                     ccGCCAATCGCCACGGTGGGTAAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGAc                                                                                                               <  2:277536/251‑1 (MQ=255)
                                                                                                                                       gCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCACCATGCCGTGCCAGCTGTTCCAGCACCTGCGGTGCAGAAGGCATCTGTTTGGTTTACGTCGGGACAACATGTCTCACCGCGATCTAAGAGGGGCGGGCTGCCCGCGCGACGTCGATTTCCCTATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGACGCGCGGCTGGCCTTGAGCTCGGTGGCTTGCAACTCTTCC‑TGCCCCAAATGGACGc                                                                                                             >  2:565544/1‑251 (MQ=255)
                                                                                                                                       cctaaTCGCCACGGTGGGATAATCCGCCAGTCGCTCAACTTGCCGTGCCAGCTGTTCCAGCTCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGc                                                                                                             <  2:433756/248‑1 (MQ=255)
                                                                                                                                            gcGCCACGGTGGGATAACCCGCCAGTCGCGCACCATGCCGTGCCAGCTGTTCCAGCCCCTTCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGTAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATAc                                                                                                        <  2:116492/250‑1 (MQ=255)
                                                                                                                                             cGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACg                                                                                                       >  1:77333/1‑251 (MQ=255)
                                                                                                                                                   ggTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGAt                                                                                                 <  1:218559/251‑1 (MQ=255)
                                                                                                                                                    gagggATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATg                                                                                                >  2:37262/3‑251 (MQ=255)
                                                                                                                                                     ggggaaaacccGCCAGTCGCTAAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACAGGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCAGATC‑GTCATGTGGCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGc                                                                                               <  2:395554/242‑1 (MQ=255)
                                                                                                                                                      gggATAACCGGCCGGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCGGTTTGGTTTGCGTCGGGAACACATGTCGCAGCGGGATATAAGAGGTGCGTGCTGCCAGCGCGACGTCGATTTCCATAGC‑TTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCTCTGATGGCATTGAGAGCGGTGGCTTGCAAATCTTCC‑TGCCCCAAAGGGACGCCATACTCCTGATGCt                                                                                              <  2:38230/251‑1 (MQ=255)
                                                                                                                                                       ggATAATCCGCCATTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGGCGGTAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCTCTACGTCGATTTCCATGTC‑GGCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGAGCGGTTGCTTGCAAATCTTCC‑TGCCCCAAATGGACTCCATACGCCTGATGCtt                                                                                             <  2:296657/251‑1 (MQ=255)
                                                                                                                                                       ggATAATCCGCCAGTCACTCAACATGCCGTGCCAGCTGCTCAGGCACCTGCGGGGCAGAAGGCCTCTGTTTGGTTTGCGCCGGGAACACCTGTATCAGCGCGATATAAGCGGGGCGTGCTGCCCGCGCGACGTCGATTTCCATGTC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCAGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCtt                                                                                             >  2:35177/1‑251 (MQ=255)
                                                                                                                                                              ccGCCAGTCGCTCAACAGGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGcc                                                                                      <  2:35301/251‑1 (MQ=255)
                                                                                                                                                                  cAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAgcc                                                                                  <  1:375927/251‑1 (MQ=255)
                                                                                                                                                                      cGCTCAACATGCCGTGCCAGCGGTTCCAGCCCCGGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAGCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCa                                                                              <  1:569788/251‑1 (MQ=255)
                                                                                                                                                                          aaacatcccgccccatcTGTCCCAGCCCCGGAGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGAGGCTTGATCGCCAGCCGCCAGTaa                                                                          <  2:423043/235‑1 (MQ=255)
                                                                                                                                                                                     tccagcgtttCCAGCCCCTGCTGTGCCGAAGGCATCTGGTTGGGTTGCGGCGGGAACACATGTCGCAGCTGGAGATAAGAGGGTCGTGCTGCCAGCGCGACGTCTATGTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGAGCGGTGGCTGGCAAATCTTCC‑GGCCCCAAATGGACGCCATACGCCGGATGCTTGACCGCCGCCCGCCAGTAATCGTTGATaa                                                                <  2:546091/242‑1 (MQ=255)
                                                                                                                                                                                        aGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACa                                                             >  1:399996/1‑250 (MQ=255)
                                                                                                                                                                                         gCTGTTCCAGCCCCTGCGGTGCAGAAAGCATCTGTTTGGTTTGCGACGGGACCACATGACTCAACCCGATATAAGAGGCGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGTAACGCCCAGACGCCAGCATGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTTATCGCCAGCCGCCAGTAATCGTTGATAAACAAt                                                           >  2:235325/1‑251 (MQ=255)
                                                                                                                                                                                             tgccAGCCCCTGTGGTGCAGAAGTCATCTGGTTGTTTTGCGTCGGGAACACATGTCTCAGCGGGATATAAGAGGGGCGGGCTGGCAGTGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAAGGCCGAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGGTTG‑AAATTGTCCGTGCCCCAAATGGACGCCTTCCGGCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATcgc                                                        <  1:527237/248‑1 (MQ=255)
                                                                                                                                                                                                cAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGCCGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGt                                                    >  2:132755/1‑251 (MQ=255)
                                                                                                                                                                                                  gCCCCTGCGGGGCAGAGGGCATCGTTTTGGTTTGCGTGGGGACCACATGTCTCCGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTta                                                  <  2:14926/251‑1 (MQ=255)
                                                                                                                                                                                                   ccccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTtat                                                 <  1:129082/251‑1 (MQ=255)
                                                                                                                                                                                                     ccggctgtgcggatggCATCTGTTTGGTGTCCGTCGGGATCACATTTCTCAGCACGATATAAGAGGGGCGAGCTGCCAGCGCGACGTCGTGGTCC‑TATCCGTCATGTGGCGAAACGCCAAGCCGCAGGCCTTCCGCGCGGATGTAATTGAGCTCTGTGGCTAGCAAATCTTCC‑TGCCCTAAATGGACGCCATCCGCCTGATGCTGGATCGCCAGCCGCCAGTAATCGTTGATAAACGATCGCGCGatattg                                               <  2:497539/237‑3 (MQ=255)
                                                                                                                                                                                                       tGCGGTGCAGACGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGCTGCTTGATCGCCAGCCGCCAGTACTCGTTGATAAACCATCGCGCGTTATAgcg                                             >  2:122768/1‑251 (MQ=255)
                                                                                                                                                                                                        gCGGTGCAGAAGGCATCTGTTTGTTTTGCTTCGGGAACACATGACACAGCGCGATAAAAGAGGGGCGTGCTGCCAGCGCGACGTCGATGTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAgcgg                                            <  1:560702/251‑1 (MQ=255)
                                                                                                                                                                                                         cGGTGCAGACGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAACTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAgcggc                                           >  1:83249/1‑251 (MQ=255)
                                                                                                                                                                                                          ggTGCAGCAGGTATCTGTTCGGTTTGCGTCGGGACCACATGTCTCAGCGCGATATCAGAGGGGCGTGCTGCAAGCGCGACGTCGATTTCCATATC‑GTCATGTGACGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTCATCATTGATAAACAATCGCGCGTTATAacggcg                                          >  2:233857/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGGCAGCTGTTTGGTTTGGGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAgcg                                  <  1:333392/251‑1 (MQ=255)
                                                                                                                                                                                                                            ttggtttgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgc                        >  1:361240/1‑251 (MQ=255)
                                                                                                                                                                                                                            ttggtttgAGTCGGGAACACCTGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGCAAAGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCCAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGCTGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTgccgc                        >  1:396931/1‑251 (MQ=255)
                                                                                                                                                                                                                             tggtttgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACCCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgcc                       >  1:67154/1‑251 (MQ=255)
                                                                                                                                                                                                                             tggtttgCGCCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgc                        <  1:132755/250‑1 (MQ=255)
                                                                                                                                                                                                                              gtttggctgcGGGAACACAGGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGGCGACGATTTCCATAGC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCGGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGAGCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAGTGCCGCCa                      <  2:100145/242‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACAt                <  1:263377/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCGAGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGCTCGGTGGCTTGCAAATCTTCC‑TGCCCCAAAAGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAAACGCGAGTTATAGCGGCGGCCCAGCGCCATTGCCGCCCCGACAt                >  2:340299/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gggAACATATTTCTCAGCGCGATATAAGAGGGGCGTGCGGCCTGCGCGACGTCGATTCCCATATC‑GTCAGGTGTCGAAACGCCCAGCCGCGGGCCTGCCGCGCGGGTGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCTATACGCCTGATGCTTGATCGCCGTCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc            <  1:100376/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gggAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc            >  2:209664/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggAACACAGGTCTCAGCGCGATATAAGAGGGGCGTGCGGCCAGCGCGACGTCGATTTCCATAGC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTGGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCGGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCt           <  2:147336/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCtt          >  2:454026/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaCACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGCTCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCACCGCAATTGCCGCCACggactcggattc         >  2:188533/1‑240 (MQ=255)
                                                                                                                                                                                                                                                  gTCCCAGCGGGATAAAAAGGGGGCGTGCTGCCAGCGCGCAGTGGATTCCCATATC‑GTCATGTGGCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACctct  <  2:292695/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
TCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATC‑GTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCC‑TGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCT  >  NC_000913/4193576‑4194067

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