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

TCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCT  >  NC_000913/4193576‑4194067
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
ccGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTTCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGGTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAgcgc                                                                                                                                                                                                                                                     <  2:217485/250‑1 (MQ=255)
  gccaaccacTGCAATTTCCAGCAACTGAGCCGTTGCCAAACGCAAGTCCGCGGCTTGAGTAATGGCGCTGAAGACGGCGATACTGCCGACACCCGTTGCTATGACCGCAGGCGCGCGTGACAGACTGATGCCGCCAATCGCCACGGGGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCAGTTCCAGCCCCTGCGGTGAAGACGGCCTCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGcgcg                                                                                                                                                                                                                                                    >  2:53092/9‑250 (MQ=255)
     aaCTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGatat                                                                                                                                                                                                                                                >  1:576890/1‑251 (MQ=255)
         ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCCATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGCACACATGTCTCAGCGCGAtata                                                                                                                                                                                                                                               >  1:100145/1‑248 (MQ=255)
         ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAag                                                                                                                                                                                                                                             <  1:472114/250‑1 (MQ=255)
            gCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAgg                                                                                                                                                                                                                                          <  1:291805/250‑1 (MQ=255)
             cAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAgggg                                                                                                                                                                                                                                        <  1:243068/251‑1 (MQ=255)
              aaTTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGc                                                                                                                                                                                                                                       >  1:473500/1‑251 (MQ=255)
               aTTTCCAGCAACTGTGCCGTTGCCAATCGCAAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCg                                                                                                                                                                                                                                      <  2:505905/251‑1 (MQ=255)
                    cAGCAACTGTGCCGTTGCCCAACCCCAGTCTGCGGCTTGAGTCATGGCGCTGACGCCGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCCCGGTGGGATCATCCGCCCGTCGCTCAACATGCCGTGCCCGCTGTTCCAGCCCCTGCGGTGCAGAATGCATCTGTTTGGTTTGCGTCGGCAACACCTGTCTCAgcgctatctaaccgggtctttcct                                                                                                                                                                                                                                 >  2:135261/1‑231 (MQ=255)
                    cAGCAACTGTGCCGTTGCCAAACGCCAGTCGGCGGCTTGAGTAAGGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGtgctg                                                                                                                                                                                                                                 <  2:338682/251‑1 (MQ=255)
                       cAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGccc                                                                                                                                                                                                                              >  2:129082/1‑250 (MQ=255)
                         aCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAgc                                                                                                                                                                                                                            <  1:453009/251‑1 (MQ=255)
                            gtgCCGTTGCCAAACGCCAGTCTGCGGCTTGGGTAATGGCGCTGACGACGGCGCTCCGGCCGACACCCGTTGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTGCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCTCGATATAAGAGGGGCGTGCTGCCAgcgc                                                                                                                                                                                                                          <  2:589360/250‑1 (MQ=255)
                             tgCCGTTGCCAAACGCCAGTCTGCGGGTTGCGTATTGGCACTGACGCCGGCGATACTGCAGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATACTGCGCCAGTCGCTCAACATGCTGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCATCGCGCTATACGAGGGGCGTGCTGCCAGcgcgc                                                                                                                                                                                                                        >  2:64253/1‑250 (MQ=255)
                                    gCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAt                                                                                                                                                                                                                 <  2:16044/251‑1 (MQ=255)
                                     ccaagctccAGTCTGCGGCTTGAGTAATGGCGCTGCCGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAtt                                                                                                                                                                                                                <  1:190041/244‑1 (MQ=255)
                                     ccAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAtt                                                                                                                                                                                                                <  1:460487/251‑1 (MQ=255)
                                     ccAAACGACACTCTGCGGCTTGAGTAATGGCGCTGACGACGCCGATACTGCCGCAACCAGTTGCTATCACCGCAGGCGAGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATACTCCGCCCGCCGCTCAACATGCCGAGCCAGCTGTTCCAGACCCTGCGGTGCAGCAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCCGCGCGGTATAAGAGGGGCGTGCTGCCCGCGCGACGCCGAtt                                                                                                                                                                                                                >  2:378499/1‑251 (MQ=255)
                                      caaccGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                                                                                                                                               <  2:481917/247‑1 (MQ=255)
                                      cAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTACCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                                                                                                                                               >  2:263377/1‑251 (MQ=255)
                                       aaaCGCCAGGCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGCCACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTc                                                                                                                                                                                                              >  1:425398/1‑251 (MQ=255)
                                          cGCCAGTCAGCGGGTTGAGTAAGGGCGCTGACGACGGCTATTCTGCCGACACCCGTTGCTAACACCGCAGGCGCGTGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGAGAATCCGCCAGGCGCTAAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCAACTGTTAGGTTTGCGTCGGGAACACATGTCTCAGTGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCat                                                                                                                                                                                                           <  1:44708/251‑1 (MQ=255)
                                           gCCAGTTTGCGGCTTAAGTAATGGTGCTGACGTTGGCTATACGGTCGCCACCCGGGTGTATCACCGCCGGCGCCGGTTCCAGACTGATGCCGCCAATCTCCACGGTTGGATAATCCTCCAGTCGCTCAACATGCCGCGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGTATCTGTGTGGTTTGCGTCTGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGGCGAtgaccata                                                                                                                                                                                                          <  2:541795/251‑8 (MQ=255)
                                               gTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGCCACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGt                                                                                                                                                                                                     >  1:210875/1‑251 (MQ=255)
                                                tCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAAAGGGGCGAGCTGCCAGCGCGACGTCGATTTCCCTATCGTc                                                                                                                                                                                                    >  1:401953/1‑251 (MQ=255)
                                                 cTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTc                                                                                                                                                                                                    >  2:569788/1‑250 (MQ=255)
                                                  ggCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAt                                                                                                                                                                                                  <  1:401030/250‑1 (MQ=255)
                                                   gCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAtg                                                                                                                                                                                                 <  1:368248/251‑1 (MQ=255)
                                                    cGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAtgt                                                                                                                                                                                                >  1:213780/1‑251 (MQ=255)
                                                        ttGAGTAAGGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTc                                                                                                                                                                                             <  2:440716/250‑1 (MQ=255)
                                                                    gctgaagccgGCGATACTGCCGACACCCATTGCTATCACCGCAGGCGCGCGAGCCAGCCTGCTGCCGCCAATCGCCACGGGGGGACAATACGCCGGTCGCCCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGCCGTCTAGTCTCATATCGTCATATGTCGAAACGCCCCGcc                                                                                                                                                                                >  2:570979/9‑251 (MQ=255)
                                                                            cgGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGGCATGTGTCGAAACGCCCAGCCGCAGGCCt                                                                                                                                                                        <  1:483063/251‑1 (MQ=255)
                                                                             ggCGATACTGCCGTCACCCGTGGCGATCACCGCGGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGGGTCGAAACGCCCAGCCGCATGCCTg                                                                                                                                                                       <  2:583866/251‑1 (MQ=255)
                                                                               cGATCCTGCCGACCCCCGTTGCGATCCCCGCAGGCGCGCGTGCCAGACGGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGcc                                                                                                                                                                     <  2:165818/251‑1 (MQ=255)
                                                                                aaTACTGCCGACACCCGTTGATATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGGCGCTCAACATGCCGTGCCAGATGTTCCAGCCCCTGCGGTGCAGGAGCCATCGGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGTGGCGTGCTGCCAGCGCGACGTCGCTTTCCCTATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGccc                                                                                                                                                                    >  2:481798/2‑250 (MQ=255)
                                                                                    cTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGAGCCAGCTGTTCCAGACCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCcgcgcg                                                                                                                                                                >  1:212143/1‑251 (MQ=255)
                                                                                         gACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAACAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATCGTCATGTGTCGAACCGCCCAGCCGCAGGCATGCCCCGCGGCTgg                                                                                                                                                           >  2:486090/1‑251 (MQ=255)
                                                                                              ccGTTGCGATCACCGCAGGCGGGCGTGCCAGACTGATGCCGACAAGCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTg                                                                                                                                                      <  1:331313/251‑1 (MQ=255)
                                                                                                gTTGCGATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCCCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTgag                                                                                                                                                    <  2:461339/251‑1 (MQ=255)
                                                                                                    gCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                                                                                                                                                 <  1:430389/251‑1 (MQ=255)
                                                                                                    gCTATCACCGCAGGCGCTCTTGCCAGACGGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCCTAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGGCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCTTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                                                                                                                                                 <  2:187380/251‑1 (MQ=255)
                                                                                                    gCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                                                                                                                                                 <  1:243289/251‑1 (MQ=255)
                                                                                                        tCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTg                                                                                                                                             >  2:137907/1‑251 (MQ=255)
                                                                                                          cccGCAGGCGCGCGTGCCAGACTGTTGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGc                                                                                                                                           <  1:559215/250‑1 (MQ=255)
                                                                                                          aCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAACCTCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGc                                                                                                                                           >  1:71286/1‑251 (MQ=255)
                                                                                                             gCAGGCGCGCGTGCCAGACTGATGCAGCCAATCCCCACGGTGGGATAATCCACCAGTCGCTCAACATGCCGTGCCAGCTGTTCCCGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGACCACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGCGATCGGTGGCTTg                                                                                                                                        >  2:574211/1‑251 (MQ=255)
                                                                                                              cAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGc                                                                                                                                       >  1:556188/1‑251 (MQ=255)
                                                                                                               aGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGCTAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGc                                                                                                                                       >  1:451203/1‑250 (MQ=255)
                                                                                                                   gctcgTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATc                                                                                                                                  >  1:461831/4‑251 (MQ=255)
                                                                                                                   gcgcgTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATc                                                                                                                                  <  2:140949/251‑1 (MQ=255)
                                                                                                                   gcgcgTGCCAGACTGATGCCGACAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGCAATGGGTCGAAACGCCCAGCCGCAGGCCTGCCGAGCGGGTGGCATTGAGATCGGTGGCTTGCaactc                                                                                                                                  >  2:193453/1‑248 (MQ=255)
                                                                                                                      ggtcccaggcagtgtccGCCAATCGCCTCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTc                                                                                                                               <  2:175292/236‑1 (MQ=255)
                                                                                                                         gCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTg                                                                                                                           <  1:111914/251‑1 (MQ=255)
                                                                                                                           cAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGcc                                                                                                                         <  1:137907/251‑1 (MQ=255)
                                                                                                                             gACTGATGCCGACAATCGCCGCGGGGGGATAATCCGCCAGGCGCTCAAAATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGAAATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGACCAGCCGCAGGCCTGCAGAGCGGATGGCATTGCGATCGGTGGCTTGCAAATCTTCCTGcccc                                                                                                                       >  2:46950/1‑251 (MQ=255)
                                                                                                                                     cctccAATCGCCACGGTGGGATAAACCGCCAGTCGCTCAACAAGCCGTGCCAGCTGTTCCAGCCCCTGCGGGGCAGAAGGCATCTGTTTGGTTTGCGTAGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCCGCCCCAAATGGAc                                                                                                               <  2:35290/248‑1 (MQ=255)
                                                                                                                                     ccGCCAATCGCCACGGTGGGTAAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGAc                                                                                                               <  2:277536/251‑1 (MQ=255)
                                                                                                                                       gCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCACCATGCCGTGCCAGCTGTTCCAGCACCTGCGGTGCAGAAGGCATCTGTTTGGTTTACGTCGGGACAACATGTCTCACCGCGATCTAAGAGGGGCGGGCTGCCCGCGCGACGTCGATTTCCCTATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGACGCGCGGCTGGCCTTGAGCTCGGTGGCTTGCAACTCTTCCTGCCCCAAATGGACGc                                                                                                             >  2:565544/1‑251 (MQ=255)
                                                                                                                                       cctaaTCGCCACGGTGGGATAATCCGCCAGTCGCTCAACTTGCCGTGCCAGCTGTTCCAGCTCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGc                                                                                                             <  2:433756/248‑1 (MQ=255)
                                                                                                                                            gcGCCACGGTGGGATAACCCGCCAGTCGCGCACCATGCCGTGCCAGCTGTTCCAGCCCCTTCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGTAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATAc                                                                                                        <  2:116492/250‑1 (MQ=255)
                                                                                                                                             cGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACg                                                                                                       >  1:77333/1‑251 (MQ=255)
                                                                                                                                                   ggTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGAt                                                                                                 <  1:218559/251‑1 (MQ=255)
                                                                                                                                                    gagggATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATg                                                                                                >  2:37262/3‑251 (MQ=255)
                                                                                                                                                     ggggaaaacccGCCAGTCGCTAAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACAGGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCAGATCGTCATGTGGCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGc                                                                                               <  2:395554/242‑1 (MQ=255)
                                                                                                                                                      gggATAACCGGCCGGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCGGTTTGGTTTGCGTCGGGAACACATGTCGCAGCGGGATATAAGAGGTGCGTGCTGCCAGCGCGACGTCGATTTCCATAGCTTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCTCTGATGGCATTGAGAGCGGTGGCTTGCAAATCTTCCTGCCCCAAAGGGACGCCATACTCCTGATGCt                                                                                              <  2:38230/251‑1 (MQ=255)
                                                                                                                                                       ggATAATCCGCCATTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGGCGGTAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCTCTACGTCGATTTCCATGTCGGCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGAGCGGTTGCTTGCAAATCTTCCTGCCCCAAATGGACTCCATACGCCTGATGCtt                                                                                             <  2:296657/251‑1 (MQ=255)
                                                                                                                                                       ggATAATCCGCCAGTCACTCAACATGCCGTGCCAGCTGCTCAGGCACCTGCGGGGCAGAAGGCCTCTGTTTGGTTTGCGCCGGGAACACCTGTATCAGCGCGATATAAGCGGGGCGTGCTGCCCGCGCGACGTCGATTTCCATGTCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCAGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCtt                                                                                             >  2:35177/1‑251 (MQ=255)
                                                                                                                                                              ccGCCAGTCGCTCAACAGGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGcc                                                                                      <  2:35301/251‑1 (MQ=255)
                                                                                                                                                                  cAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAgcc                                                                                  <  1:375927/251‑1 (MQ=255)
                                                                                                                                                                      cGCTCAACATGCCGTGCCAGCGGTTCCAGCCCCGGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAGCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCa                                                                              <  1:569788/251‑1 (MQ=255)
                                                                                                                                                                          aaacatcccgccccatcTGTCCCAGCCCCGGAGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGAGGCTTGATCGCCAGCCGCCAGTaa                                                                          <  2:423043/235‑1 (MQ=255)
                                                                                                                                                                                     tccagcgtttCCAGCCCCTGCTGTGCCGAAGGCATCTGGTTGGGTTGCGGCGGGAACACATGTCGCAGCTGGAGATAAGAGGGTCGTGCTGCCAGCGCGACGTCTATGTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGAGCGGTGGCTGGCAAATCTTCCGGCCCCAAATGGACGCCATACGCCGGATGCTTGACCGCCGCCCGCCAGTAATCGTTGATaa                                                                <  2:546091/242‑1 (MQ=255)
                                                                                                                                                                                        aGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACa                                                             >  1:399996/1‑250 (MQ=255)
                                                                                                                                                                                         gCTGTTCCAGCCCCTGCGGTGCAGAAAGCATCTGTTTGGTTTGCGACGGGACCACATGACTCAACCCGATATAAGAGGCGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGTAACGCCCAGACGCCAGCATGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTTATCGCCAGCCGCCAGTAATCGTTGATAAACAAt                                                           >  2:235325/1‑251 (MQ=255)
                                                                                                                                                                                             tgccAGCCCCTGTGGTGCAGAAGTCATCTGGTTGTTTTGCGTCGGGAACACATGTCTCAGCGGGATATAAGAGGGGCGGGCTGGCAGTGCGACGTCGATTTCCATATCGTCATGTGTCGAAAGGCCGAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGGTTGAAATTGTCCGTGCCCCAAATGGACGCCTTCCGGCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATcgc                                                        <  1:527237/248‑1 (MQ=255)
                                                                                                                                                                                                cAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGCCGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGt                                                    >  2:132755/1‑251 (MQ=255)
                                                                                                                                                                                                  gCCCCTGCGGGGCAGAGGGCATCGTTTTGGTTTGCGTGGGGACCACATGTCTCCGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTta                                                  <  2:14926/251‑1 (MQ=255)
                                                                                                                                                                                                   ccccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTtat                                                 <  1:129082/251‑1 (MQ=255)
                                                                                                                                                                                                     ccggctgtgcggatggCATCTGTTTGGTGTCCGTCGGGATCACATTTCTCAGCACGATATAAGAGGGGCGAGCTGCCAGCGCGACGTCGTGGTCCTATCCGTCATGTGGCGAAACGCCAAGCCGCAGGCCTTCCGCGCGGATGTAATTGAGCTCTGTGGCTAGCAAATCTTCCTGCCCTAAATGGACGCCATCCGCCTGATGCTGGATCGCCAGCCGCCAGTAATCGTTGATAAACGATCGCGCGatattg                                               <  2:497539/237‑3 (MQ=255)
                                                                                                                                                                                                       tGCGGTGCAGACGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGCTGCTTGATCGCCAGCCGCCAGTACTCGTTGATAAACCATCGCGCGTTATAgcg                                             >  2:122768/1‑251 (MQ=255)
                                                                                                                                                                                                        gCGGTGCAGAAGGCATCTGTTTGTTTTGCTTCGGGAACACATGACACAGCGCGATAAAAGAGGGGCGTGCTGCCAGCGCGACGTCGATGTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAgcgg                                            <  1:560702/251‑1 (MQ=255)
                                                                                                                                                                                                         cGGTGCAGACGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAACTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAgcggc                                           >  1:83249/1‑251 (MQ=255)
                                                                                                                                                                                                          ggTGCAGCAGGTATCTGTTCGGTTTGCGTCGGGACCACATGTCTCAGCGCGATATCAGAGGGGCGTGCTGCAAGCGCGACGTCGATTTCCATATCGTCATGTGACGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTCATCATTGATAAACAATCGCGCGTTATAacggcg                                          >  2:233857/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGGCAGCTGTTTGGTTTGGGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAgcg                                  <  1:333392/251‑1 (MQ=255)
                                                                                                                                                                                                                            ttggtttgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgc                        >  1:361240/1‑251 (MQ=255)
                                                                                                                                                                                                                            ttggtttgAGTCGGGAACACCTGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGCAAAGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCCAATCTTCCTGCCCCAAATGGACGCCATACGCCTGCTGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTgccgc                        >  1:396931/1‑251 (MQ=255)
                                                                                                                                                                                                                             tggtttgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACCCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgcc                       >  1:67154/1‑251 (MQ=255)
                                                                                                                                                                                                                             tggtttgCGCCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgc                        <  1:132755/250‑1 (MQ=255)
                                                                                                                                                                                                                              gtttggctgcGGGAACACAGGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGGCGACGATTTCCATAGCGTCATGTGTCGAAACGCCCAGCCGCAGGCCGGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGAGCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAGTGCCGCCa                      <  2:100145/242‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACAt                <  1:263377/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCGAGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGCTCGGTGGCTTGCAAATCTTCCTGCCCCAAAAGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAAACGCGAGTTATAGCGGCGGCCCAGCGCCATTGCCGCCCCGACAt                >  2:340299/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gggAACATATTTCTCAGCGCGATATAAGAGGGGCGTGCGGCCTGCGCGACGTCGATTCCCATATCGTCAGGTGTCGAAACGCCCAGCCGCGGGCCTGCCGCGCGGGTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCTATACGCCTGATGCTTGATCGCCGTCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc            <  1:100376/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gggAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc            >  2:209664/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggAACACAGGTCTCAGCGCGATATAAGAGGGGCGTGCGGCCAGCGCGACGTCGATTTCCATAGCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTGGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCGGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCt           <  2:147336/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCtt          >  2:454026/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaCACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGCTCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCACCGCAATTGCCGCCACggactcggattc         >  2:188533/1‑240 (MQ=255)
                                                                                                                                                                                                                                                  gTCCCAGCGGGATAAAAAGGGGGCGTGCTGCCAGCGCGCAGTGGATTCCCATATCGTCATGTGGCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACctct  <  2:292695/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
TCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCT  >  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: