Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,619,080 G→T L284L (CTG→CTT)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,619,0800GT98.9% 277.1 / NA 88L284L (CTG→CTT) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (39/48);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑GGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTG  >  NC_000913/4618834‑4619321
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
cGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTAtct                                                                                                                                                                                                                                                >  1:530473/1‑251 (MQ=255)
cGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTAtct                                                                                                                                                                                                                                                <  2:530473/251‑1 (MQ=255)
 gTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTAtctc                                                                                                                                                                                                                                               <  1:459372/251‑1 (MQ=255)
  ggAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGTTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTTAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTTGGTGATGTCACGATGGCGCTGGGCGGGGAGATGTTTATCTcc                                                                                                                                                                                                                                              <  2:591890/250‑1 (MQ=255)
   gAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACTATGGCGCTGTGCGTGGAGATGCTTATCTCCg                                                                                                                                                                                                                                             <  1:536594/251‑1 (MQ=255)
     aaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGc                                                                                                                                                                                                                                           >  1:477577/1‑251 (MQ=255)
     aaGTGGGTAGCGGCGCGTCTATGCCGACCTCCGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCCCCACCGCGCTGCTCACCGACATAAATCAGGGACTGGCCTCCCGTGCAGGTAACGCGGTTGCAGTTCGTGAAGCGGGGCAGGTCCTGACGGGGGAGTATAGTAACCCGCGGCTGTTTGCGGTCCCGATGGCGATGTGCGGGGAGATGCTGAGCTCCGGc                                                                                                                                                                                                                                           >  2:2064/1‑251 (MQ=255)
       gTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCaa                                                                                                                                                                                                                                         <  1:509016/251‑1 (MQ=255)
          ggTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGCAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGACTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGGGGAGATGCTTATCTCCGGCAAACt                                                                                                                                                                                                                                      >  2:504908/1‑251 (MQ=255)
          ggTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACt                                                                                                                                                                                                                                      >  1:133835/1‑251 (MQ=255)
                ggcgCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGa                                                                                                                                                                                                                                 >  1:190454/1‑250 (MQ=255)
                 gcggcttgaTGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGTCAGGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTGGAATTGCGTGAAGCGGTGCTGTTCCTGACGGTTGAGTATCGTAACCCGCTTCTGTTTGATGTCACGATGGCGCTGTGCGTGGTGATGCTTAGCTCCGTCAAACTGGCGaaa                                                                                                                                                                                                                               <  2:372847/243‑1 (MQ=255)
                   gcgTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGaa                                                                                                                                                                                                                                >  2:108980/1‑248 (MQ=255)
                      tttATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATg                                                                                                                                                                                                                           >  2:339168/1‑250 (MQ=255)
                       ttATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGa                                                                                                                                                                                                                          <  1:419566/250‑1 (MQ=255)
                       ttATGCCGACCTACGAACTCGCTGAAGCCCGTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCGCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGa                                                                                                                                                                                                                          <  2:95857/250‑1 (MQ=255)
                           gcctaCCTACGAACGCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGGCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGcc                                                                                                                                                                                                                      <  2:172763/246‑1 (MQ=255)
                              gACCTACGCACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGCGGAGATGCTTATCTCCGGCACACTGGCGAAAGATGACGCCGa                                                                                                                                                                                                                    >  2:77422/1‑249 (MQ=255)
                                   tcGAACTCTCTGAAGCCCTTGCCGAAGCGATTGGTGGCGGGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGGTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGGGCGTGGAGAGGCTTATCTCCGGGAAACTGGCGAAAGATGACGCCGAAgcgc                                                                                                                                                                                                               <  2:48214/248‑1 (MQ=255)
                                     gAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgcgc                                                                                                                                                                                                           >  2:223451/1‑251 (MQ=255)
                                      aaCTCTCTGAAGCCCGTGCCGAAGCGAGTGGTGGCGTGGCTAACGGCGCGGGCGTGCGCACCACCGCGCGGCTGACCGACATGAATCAGGGACGGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTTTTGAAGCGGGGCAGTTGCTGACGGGTGAGTATCGTAACCCGCGTCTGGTTGATGTCAGGATGGCGCTGTGCGTGGATATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGcgcgcgcg                                                                                                                                                                                                          <  1:621855/251‑1 (MQ=255)
                                         tctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCCGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCCCGATGGCGCTGTGCGTGGAGATTCTTATCTCCGGAAAACTGGCGAAATATGACGCCGAAGCGCGCGCGaa                                                                                                                                                                                                        >  2:488618/1‑250 (MQ=255)
                                             tGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGc                                                                                                                                                                                                   >  1:515860/1‑251 (MQ=255)
                                              gAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCa                                                                                                                                                                                                  <  1:156873/251‑1 (MQ=255)
                                               gtgCCCTTGCCGAGTCGATCGTTGGCGTGGCTGTCGGAGTGGGCGTGCGCACTACCCCCCTGCTACCCGACATGAAGCATGTACTGGCCTCCAGTGCAGGTAACGCGTTTGAAGGTCGTGAAGCGGTGCATTTCCTGACGGGTGAGGATCGTTACCCGCGTCTGTTTGATGTCACGATGGCGCGGTGCGGTGAGATGCTTATCTCCTGCAAACTGGCTAAAGATGActattaagtgcgcgcgacatttct                                                                                                                                                                                                  <  2:367100/248‑24 (MQ=255)
                                                  cccTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCaagca                                                                                                                                                                                              <  1:223451/251‑5 (MQ=255)
                                                      tGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGAAGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGc                                                                                                                                                                                         >  1:229651/1‑251 (MQ=255)
                                                         cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCAAGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGACGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGCTGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTgg                                                                                                                                                                                      >  1:411657/1‑251 (MQ=255)
                                                          gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGACGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTgg                                                                                                                                                                                      >  1:83665/1‑250 (MQ=255)
                                                          gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTgg                                                                                                                                                                                      >  1:46069/1‑250 (MQ=255)
                                                              agATTGTTGGCGTGGCTAACGGCGCTGGCGTGAGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGACAGATGACGCAGAAGCGCGCGCGCAATTGCA‑AGCATGGCTGGACAAc                                                                                                                                                                                 >  2:332586/2‑251 (MQ=255)
                                                               gATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACg                                                                                                                                                                                <  1:347686/251‑1 (MQ=255)
                                                                  tgttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGATGCTCCCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGAGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGt                                                                                                                                                                              >  2:553892/1‑250 (MQ=255)
                                                                     tgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTaaa                                                                                                                                                                           <  1:507785/250‑1 (MQ=255)
                                                                       gCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCg                                                                                                                                                                        <  2:167158/251‑1 (MQ=255)
                                                                        cGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCgg                                                                                                                                                                       <  1:553892/251‑1 (MQ=255)
                                                                               aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAg                                                                                                                                                                 <  1:548592/250‑1 (MQ=255)
                                                                               aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGt                                                                                                                                                                <  1:356066/251‑1 (MQ=255)
                                                                                aCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGCAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGaagaaggc                                                                                                                                                               >  2:222720/1‑249 (MQ=255)
                                                                                 cGGCGCTGGCGTGCGCACCAACGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGACGCGGAGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGATGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAACGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCCGTGCTGGACAACGGGAAAGCGGCAGACGTCt                                                                                                                                                              >  2:542879/1‑251 (MQ=255)
                                                                                     gcTGGCGTGCGAAACACCGAGCTGCTCACAGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGATGAAGTTCGTGAAGCGGTGCCGTTCCTGGCGGGTGCGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGGGCGTGGAGATGCTTATCTACGGCAAACTGGCGCAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGAAGAAGTATTTgg                                                                                                                                                          >  2:558002/1‑251 (MQ=255)
                                                                                      cTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGCTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTTGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTCGGt                                                                                                                                                         >  1:108740/1‑251 (MQ=255)
                                                                                                       gcgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGcaca                                                                                                                                        <  2:241846/251‑1 (MQ=255)
                                                                                                                aCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGATGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGAAGCGGGACAAAAAAGccc                                                                                                                               >  1:59242/1‑251 (MQ=255)
                                                                                                                       tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGGTCGTGAAGCGGTGCAGTTCATGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGATGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGCCGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTAGGGGAGCGGCACaaaacagccccacggac                                                                                                                        >  2:195975/1‑238 (MQ=255)
                                                                                                                              ggTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGGTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTg                                                                                                                 <  2:229651/251‑1 (MQ=255)
                                                                                                                                 ccTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCGGACGGGTGAGTATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTgaga                                                                                                              <  2:352220/250‑1 (MQ=255)
                                                                                                                                 aCTGGCCTCCAGTGCAGGTAACGCCGTTGAAGTTCGTGACGCGGTCCAGTTCATGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCGCGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGCAAGCTGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGAATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTg                                                                                                                 >  2:618561/1‑248 (MQ=255)
                                                                                                                                   tGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGCTGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAAc                                                                                                            >  2:42293/1‑251 (MQ=255)
                                                                                                                                                        cgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCTTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACa                                                                                       <  1:431748/251‑1 (MQ=255)
                                                                                                                                                         ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGAc                                                                                        >  1:450275/1‑249 (MQ=255)
                                                                                                                                                         ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGAc                                                                                        >  1:541385/1‑249 (MQ=255)
                                                                                                                                                             gAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGAt                                                                                  <  1:108980/251‑1 (MQ=255)
                                                                                                                                                                  tCGTGAAGCGGTGCAGTTCTTGACGTGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGGGTGGAGATGCGTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAAGCGATTTCGTTGGGAACTACGCGAAGTATCTGCCCACAGCGATGCTg                                                                              <  2:429027/250‑1 (MQ=255)
                                                                                                                                                                   cGTGAAGCGGTGCAGTTCCTGACGGTTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGTTTTCGTTGAGAACTACGCGAAGTAGCTGCCGACAGCGATGCTGAc                                                                            <  1:503959/251‑1 (MQ=255)
                                                                                                                                                                   cGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGCAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                                                                            >  1:257609/1‑251 (MQ=255)
                                                                                                                                                                   agTGAAGCGGTGCAGTTCCTGAAGGGTGAGTCTCGTAACACGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGAACATAAAGGCCCAACCGACTTCGTTGAGAACTCAGCGAAGTATAGACCGACAGCGATGCTGa                                                                             >  2:416431/2‑250 (MQ=255)
                                                                                                                                                                                aGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAt                                                               >  1:351969/1‑251 (MQ=255)
                                                                                                                                                                                    ccTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAATATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGTCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTg                                                           <  2:492749/251‑1 (MQ=255)
                                                                                                                                                                                        gcGGGTGAGTGTCGTAACCCGCGGCTGTGTGATGTCACGATGGCGCTGTGCGGGGAGATGCGTAGCTCCGGTAAAGTGGCGAAAGATGACGGCGAAGCGCGCGCGAAATGGCAAGGAGG‑GTTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTTGTATGGTAGCTGCACAAAAAGGCCCAAGCGTCTTCGTTGAGAACTGCGCGAAGTATCGGCCGACAGTGAGGCTGTCGAAAGCACTCTTTGCTGATAc                                                       <  2:568568/250‑1 (MQ=255)
                                                                                                                                                                                         cGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAcc                                                      <  2:83665/251‑1 (MQ=255)
                                                                                                                                                                                          gggTGAGTAGCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTTGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATGTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCg                                                     <  1:266720/251‑1 (MQ=255)
                                                                                                                                                                                           ggTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGa                                                    <  1:432227/251‑1 (MQ=255)
                                                                                                                                                                                              gattatCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAg                                                  <  2:254860/247‑1 (MQ=255)
                                                                                                                                                                                                     gTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTc                                          >  1:482498/1‑251 (MQ=255)
                                                                                                                                                                                                          ccGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTg                                      >  2:616974/1‑250 (MQ=255)
                                                                                                                                                                                                           cgcgTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGTCGAAGCTCGCGCGAAATTGCA‑TGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCATCCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATTCTGATACCGAAGGTTTTGTTAGTGaa                                    <  2:123000/251‑1 (MQ=255)
                                                                                                                                                                                                            gcgTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGaaa                                   <  2:153408/251‑1 (MQ=255)
                                                                                                                                                                                                            gcgTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGaaa                                   <  1:190327/251‑1 (MQ=255)
                                                                                                                                                                                                                cggTTTGATGTCTCGATGGCGCGGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGTTGATACCGAAGGTTTTGTCAGTGAAATgg                                <  1:209565/248‑1 (MQ=255)
                                                                                                                                                                                                                  gTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                              >  1:600478/1‑250 (MQ=255)
                                                                                                                                                                                                                  gTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAGGCTGATACCGAAGGTTTTGTCAGGGAAATgg                                <  1:162722/248‑1 (MQ=255)
                                                                                                                                                                                                                     tGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGTCAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAccc                          >  1:532582/1‑251 (MQ=255)
                                                                                                                                                                                                                      gATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAccc                          <  1:540913/250‑1 (MQ=255)
                                                                                                                                                                                                                       aTGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTCAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgc                        >  2:39987/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgcg                       >  2:503947/1‑251 (MQ=255)
                                                                                                                                                                                                                           cACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgcgcg                     <  1:542879/250‑1 (MQ=255)
                                                                                                                                                                                                                            aCGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCt                   <  1:488618/251‑1 (MQ=255)
                                                                                                                                                                                                                              gagtgcgcTGTGCGTGGAGATGCTTATCTCCGGCAACCTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGTTAGCGGCACAAAAAGGCCCAACCGTCTTCGTTGAGAACGACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCTGTCTATGCTGATACCTAAGGTTTTGTCAGTGAAATGGATACCCCCGCGCTgg                 <  2:368932/247‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTgggg               <  1:77422/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cgcTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGGCGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGTTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATg            <  1:359188/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tgtgCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                              >  2:295829/1‑182 (MQ=255)
                                                                                                                                                                                                                                      tgtgCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                              <  1:295829/182‑1 (MQ=255)
                                                                                                                                                                                                                                        ggCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGt       <  2:482498/250‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGGAGATGCTTATCTACGGCAAACTGGCGAGAGATGACGGCGAAGCGCGCGAGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGACTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGACGACAGCGATGCTGCCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAACCGCGCGCTGGGGATGGCAGTGGt    >  2:400796/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ggAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGTGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAATGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg  <  1:630110/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ggAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg  >  2:9074/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ggAGATGCTTAGCTCCGGCAAACTGGCGAAAGAGGACGCCGAAGCGTGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACGACGCGAAGGATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTGGTCATTGAAATGGATACCCGCTCGCTGGGGATGGCAGTGGTTg  <  2:434436/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
CGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑GGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTG  >  NC_000913/4618834‑4619321

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