Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,140,198 T→G 15.1% H675P (CAC→CCC)  frwA ← putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,140,1980TG15.1% 244.5 / 3.3 110H675P (CAC→CCC) frwAputative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base T (36/47);  new base G (23/3);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 3.91e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

ATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTGCACGGCGTAATCGAGCG  >  NC_000913/4139965‑4140435
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
aTTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCCCTTTCCAGGTGTATGCACTCGGCGGTGACCAGCGGTGCGGCTTGTTGGGTCATGCGGAATTTCGCCAGCAGGTGGTCt                                                                                                                                                                                                                                           >  1:1071011/1‑239 (MQ=255)
 tttGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACAAGCGGTGCGGCTTGTTGGGTCCTGCGGAATTGGCCCAGCAGGGGTTCt                                                                                                                                                                                                                                           >  1:546934/1‑238 (MQ=255)
   tGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGGGCGGCTTGTTGGGTCATGCGGAATTGCGCCCGCAGGGGtttcact                                                                                                                                                                                                                                        >  1:701033/1‑234 (MQ=255)
     cTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGGGCACTTCTTCTTTGCTGCGCCAGGCGATTGCCAGGGAGCTGCCCTCGGCGGCGACCCGCGGTCCGTATTGTTTGTGCGTGCGGAATTTGCGCAGCAGGGGGTCTagctt                                                                                                                                                                                                                                      >  2:299405/1‑235 (MQ=255)
       gtgtgGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGGGCGGCTTGTTGGGTCATTCGGAATTGCGCCCGCAGGGGTTTTACTTcc                                                                                                                                                                                                                                     >  1:542106/1‑238 (MQ=255)
        tgtgGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGACTTGCGCCAGCAGGGGTTCTCCtttccg                                                                                                                                                                                                                                   >  1:281569/1‑235 (MQ=255)
        tgtgGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCATCGCCAATCGCTTTCCAGGGTGATGCACTCGGCGGCTGACCGCCGTGCGGCTGTTTGGGGCATTGGGAATTGCGCCCGCAGGGGTTTTTATTCCAg                                                                                                                                                                                                                                   >  2:625067/1‑239 (MQ=255)
        tgtgGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTCCTGCGCCAGTCGCTTTCCAGGGTGATGCACTCGGCGGTAACCAGCGGTGCGTCTTGTTGGGTCAGGCGGAATTGCGCCCACCGGTGTTCTACATCCa                                                                                                                                                                                                                                    >  1:279776/1‑238 (MQ=255)
               tGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGGGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTa                                                                                                                                                                                                                                          >  1:732048/1‑225 (MQ=255)
                  cAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCACCGGGGCGCCTTGTTGGGTCCTGCGGAATTTCGCCACCAGGTGTTTTCCTTCCAGCGAAGAACg                                                                                                                                                                                                                          >  1:142772/1‑238 (MQ=255)
                   aaatCTAAAGCCCCGCCCGGGAGGCAAAACGGCTTCCGCCGCCCACAGTGCGTCTTCCAGTTTCCGGGGTTAGGGGCAGCGGTCCGCCAGGTGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGGGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTTCTTCCAGCGAAGTACGGc                                                                                                                                                                                                                        <  2:408013/239‑1 (MQ=255)
                   aaaaCTAAAGCCCAGCCCGGTAGAGAAACCGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGc                                                                                                                                                                                                                        <  1:870808/239‑1 (MQ=255)
                    aaaCTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGCCTTGTTGGGGCATTCGGAATTTCGCCATCAGGGATTCTCCTTCCACCGAAGAACggga                                                                                                                                                                                                                       >  1:863075/1‑237 (MQ=255)
                      aCTAAAGCCCAGACGGGTAGAGAAACCGGCCTCGCGCGCCCACAAGTCGGTTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACgg                                                                                                                                                                                                                         <  1:994396/235‑1 (MQ=255)
                          aaGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTTCGCCAGCCGGGGTTCTCCTTCCAGCGCAGTACGGCAGGCCAt                                                                                                                                                                                                                 >  1:19824/1‑239 (MQ=255)
                           aGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGAGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTTTACTTCCACCGAAGTACGGGAGGccct                                                                                                                                                                                                                 >  1:732802/1‑236 (MQ=255)
                              ccAGACGGGTAAAAAAAACGCCCGCGCGCGCCCACAAGTCGGTTTCCAGTTTCCGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGcc                                                                                                                                                                                                             <  1:1243111/239‑1 (MQ=255)
                                 gACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTTACCAGCGGTGCGTCCTTTTGGGTCATGCGGAATTGCGCCAGCAGGAGTTCTACCTCCAGCGAAGTAAGGCAGGCCATTTCCCgg                                                                                                                                                                                                          >  2:852057/1‑239 (MQ=255)
                                 gACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCCTGCGGGATTGCGCCAGCCGGGGTTCTACTTCCAGCGAAGAACGGCAGGCCCTTGCCCgg                                                                                                                                                                                                          >  1:732689/1‑239 (MQ=255)
                                  aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTTTTGGGTCATGCGGAATTGCGCCAGCAGGGGTGCTACTTCCAGCGAAGCACGGCAGGCCATTGCCTGGt                                                                                                                                                                                                         >  1:1125965/1‑239 (MQ=255)
                                  aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGGTCTACTTCCAGCGAAGTACGGCCGGCCATTGCCTggg                                                                                                                                                                                                         >  1:1144492/1‑238 (MQ=255)
                                  aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGt                                                                                                                                                                                                         >  1:301431/1‑239 (MQ=255)
                                       tAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTTGGTCATGCGGCATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGc                                                                                                                                                                                                    >  1:799668/1‑239 (MQ=255)
                                       tAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATTCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGc                                                                                                                                                                                                    >  1:1029250/1‑239 (MQ=255)
                                       tAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTCGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTAGGGCAGGCCCTTGCCTGGTTGAga                                                                                                                                                                                                    >  1:635101/1‑238 (MQ=255)
                                            aaaaCGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGGTCTACTTCCAGCGAAGTAAGGCAGGCCCTTGCCTGGTTGAGCAACtt                                                                                                                                                                                               >  1:1292839/1‑239 (MQ=255)
                                               cccgCCCCGCGCCCCCCAAAGGCGGCTTCCAGTTTGCGGGGGTTGCGGACGCGGCCCGCCCGCCGCAGGTTATCGGTCATGCCTTTGAGCCCTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGGGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCg                                                                                                                                                                                            <  2:767920/236‑1 (MQ=255)
                                               aCGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTCCTTCCAGCGAAGTACGGCAGGCCATTTCCTGGTTGAGCAACTTGc                                                                                                                                                                                             >  1:895992/1‑238 (MQ=255)
                                                 ggCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGCATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGc                                                                                                                                                                                          >  1:117585/1‑239 (MQ=255)
                                                      gcggcccccAAAAGCCGGCTTCCAGTTTCGGCGGATCGCGGCAGGGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTTTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCa                                                                                                                                                                                     <  1:1253023/233‑1 (MQ=255)
                                                      cgcgcgCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCa                                                                                                                                                                                     <  1:161254/239‑1 (MQ=255)
                                                      cgcgcgCCCAAAAGTCGGCTTCCAGTTTCCGCGGATAGCGGCAGCGGCCCGCCAGCAGGAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCa                                                                                                                                                                                     <  1:1223270/239‑1 (MQ=255)
                                                              ccaaggTCGGGTTCCATTTTTGTGGGTTAGCGGGGGCGCCCCGCCAGGAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATc                                                                                                                                                                             <  2:1172675/234‑1 (MQ=255)
                                                                           ccAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTTCGGCACTCACGGCTATCAAGTTGCGCCAt                                                                                                                                                                >  1:47365/1‑238 (MQ=255)
                                                                             aGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCCTCCg                                                                                                                                                             >  1:969804/1‑239 (MQ=255)
                                                                                       gATGGCGGCAGCGGCCCGCCAGGAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTcgc                                                                                                                                                    <  1:1051846/238‑1 (MQ=255)
                                                                                         tAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTTATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCCGCGAAATACGGCAGGCCATTACCCGGGTGAGCAACTTTGGGCACCCCCGGCTATCAACGTTGCCCATCCGAGGttttcgcgg                                                                                                                                                 >  2:553814/1‑233 (MQ=255)
                                                                                          ggcgggggcggcCCGCCGGGGGCGGGTTTGCGGTTATGGCGTTGAGCCCTTCTTCTTGTGTGCGGCCGGCGCTTTCCCGTGGGATGCCCGCGGCGGTGGCCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCgccgcc                                                                                                                                                <  2:630117/232‑1 (MQ=255)
                                                                                          aggggaagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCgccgcc                                                                                                                                                <  1:291381/233‑1 (MQ=255)
                                                                                          agcggcagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTTCGCCATCCGAGCTTTCgccgcc                                                                                                                                                >  1:456658/1‑239 (MQ=255)
                                                                                                  cggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGGGCTTTTCCCCCCGGGATTgg                                                                                                                                        >  1:636507/1‑238 (MQ=255)
                                                                                                     cccGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATg                                                                                                                                      <  1:609341/238‑1 (MQ=255)
                                                                                                       cgcatggcgcagGTTATCGGTCATGCCTTGGAGCACGTCTTCTTTGCTGCGGCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTg                                                                                                                                   <  1:1082615/231‑1 (MQ=255)
                                                                                                       ccgcggcagctgTTTTTCCGTTATGCCTTTGAGTACTTCTGTTTTGGGGGGCCAGTCGTGTTCCAGTGTGAAGACCTCGGCGGTGACCAGCGGTTCGTCTTGTGGGGTCATGCGGAATTGCGCCAGCAGGTGTTCGACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTg                                                                                                                                   <  2:47365/234‑1 (MQ=255)
                                                                                                          cagcagcagGTTCTCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTgcgc                                                                                                                                <  1:359879/239‑1 (MQ=255)
                                                                                                               gcagGGTCTCGGTCTTTCCTTGGATCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGGCGTCGGGGGTGACCCGCGGTGCGTCTTGTGGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAAGTGATGGTGCGCTCATg                                                                                                                           <  2:1189647/239‑1 (MQ=255)
                                                                                                                 atgttgttggTCAGTCCTTTGAGCCCTGCTTGTTGGCTGCGCCAGTCGATTTGCAGTTGGGTGCACTCGGCGGTGACCAGCGGTGCGTCTTGTGGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTGCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGc                                                                                                                          <  2:713082/230‑1 (MQ=255)
                                                                                                                 agGTTACCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATg                                                                                                                           <  1:380699/237‑1 (MQ=255)
                                                                                                                     tatctggcgtGCCTTTGAGCACTTCTCCTTGGATGGGCCAGTCGCTTCCCAGTGTGATGCCCCCGGCGGTGTCCAGCGGTGCGTCTTGTGGGGTCATGCGGACTGTCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCAGTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGt                                                                                                                     <  2:172842/230‑1 (MQ=255)
                                                                                                                          gttgtgtcTTTGAGCTCTTTTTTTTTGCTGCGGCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATc                                                                                                                <  2:344112/232‑1 (MQ=255)
                                                                                                                           ttctGCTTTTGAGCACTTCTGCTTGGCTGGGCCAGTGGCTTTCCAGTGGTAGGCACGCGGCGGTGACCAGTGGTTCGTCTTGTTGGGTCATGCGGAATTGTGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATc                                                                                                                <  2:927894/235‑1 (MQ=255)
                                                                                                                           tCATGCCTTTGAGCACTTCTTCTTTGCTGCGGCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGAACAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATc                                                                                                                <  1:582919/238‑1 (MQ=255)
                                                                                                                             aTGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTgg                                                                                                                                                                                                          >  2:55026/1‑147 (MQ=255)
                                                                                                                             aTGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCt                                                                                                                               >  1:250276/1‑221 (MQ=255)
                                                                                                                                 cTTTGAGCACTGCTGCTTTGCTGCGCCAGTCGCTTGCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCAGGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCAc                                                                                                                                                                                    <  2:323417/165‑1 (MQ=255)
                                                                                                                                         cttttttttttGCTGCGCGTGTCGGTTTCCAGTGTGCTGCACGCTGCGGTGTCCCGCGGCGCGTTGTGTTGGGTCAGGCGGAATTGCCACAGCCGGTGGTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGAc                                                                                                 <  2:249076/231‑1 (MQ=255)
                                                                                                                                           ttcttcttTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCCGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAg                                                                                                                                                                         >  1:323417/1‑165 (MQ=255)
                                                                                                                                           ttcttcttTGCTGCGCCAGTCCCTTTCCAGTTTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAAt                                                                                                                                           <  1:941634/195‑1 (MQ=255)
                                                                                                                                             cttcttTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGc                                                                                             >  1:741786/1‑239 (MQ=255)
                                                                                                                                              ttcttTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGt                                                                                                                     <  1:857425/214‑1 (MQ=255)
                                                                                                                                                  ttGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGCCCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACgg                                                                                        <  1:600875/239‑1 (MQ=255)
                                                                                                                                                  ttGCGGGGCCAGTCGCTTTCCAGTGGGATGCCCTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACgg                                                                                        <  1:625067/239‑1 (MQ=255)
                                                                                                                                                   tgctgcGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCCGGCCATTGCCTGGGTGAGCCAATTGGCGCCCCCCCGGCCATCAAGGTTCGCCATCCGAGCTTTCGCCGCCCGCATTGGTTGTGCCCGCCTGCTGCGGTCCCCCACCCCCAAGCCCACCAACAACGgc                                                                                       >  2:1007366/1‑239 (MQ=255)
                                                                                                                                                        ggccaggtggtttCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCAc                                                                                  <  1:600683/229‑1 (MQ=255)
                                                                                                                                                        cgcCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGGGCCATCCGGGCTTTCGCCCCCGGAATTTATGGTGCGCCCCATCTGATTTATTCCCGCCCCAACCCGCCCAACAACGGGCGCAc                                                                                  >  2:949449/1‑239 (MQ=255)
                                                                                                                                                          ccAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGGGGCCATGCTGGGTTCATCCAGCCCTAAGCCGACCAGCAACggacgcccgg                                                                                >  1:19858/1‑231 (MQ=255)
                                                                                                                                                                  tttCCAAGGTGATGGCCTCGGCGGTGACCAGCGGTGCCTCTTTTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTc                                                                        <  1:14838/239‑1 (MQ=255)
                                                                                                                                                                         tgtgATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGGGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGGATTGATGGTGCGCTCATGCTGAGGTCATCCAGCCCTAAGCCGACCAGCAACGGCAGGACGGAACCTTtccctccc                                                                  >  1:820195/1‑238 (MQ=255)
                                                                                                                                                                           tgATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAg                                                                >  1:456687/1‑238 (MQ=255)
                                                                                                                                                                               tGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCAACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTc                                                             >  1:856743/1‑238 (MQ=255)
                                                                                                                                                                                 cACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTTAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGGCCAGCAACGGCAGCACGGAACCTTTCGCTCCCCGccccc                                                           >  1:1009570/1‑234 (MQ=255)
                                                                                                                                                                                    tCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGca                                                       <  1:205257/239‑1 (MQ=255)
                                                                                                                                                                                    tCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGca                                                       <  1:444139/239‑1 (MQ=255)
                                                                                                                                                                                    ccggcggTGACCAGCGGTGCGTGTTGTTGGGTCATGCGGAATTGCGCCAGGAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGca                                                       <  1:202100/238‑1 (MQ=255)
                                                                                                                                                                                     ctggcgtttgccACCGGTTCTTTTTTTTGGGTCCTGCGGAATTGGGCCTGCAGGTGTTCTCCTCCCAGCGTAGTACGGCAGGCCATTTCCTGGTTGGGCATCGTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGca                                                       <  2:402331/229‑1 (MQ=255)
                                                                                                                                                                                     cggggtggggccAGGGGTGCGTTTTGTGGGGTCACGGGGAATTGCGCCGTCAGGTGTTCTTCTGCCCGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGca                                                       <  2:741786/229‑1 (MQ=255)
                                                                                                                                                                                      ggcggTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGGGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCCCCCCCACCTCAcccccc                                                      >  1:691175/1‑234 (MQ=255)
                                                                                                                                                                                      ggcggTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTTCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCAcccccca                                                     >  1:777852/1‑234 (MQ=255)
                                                                                                                                                                                        cggTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTCAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAg                                                    >  1:1015229/1‑238 (MQ=255)
                                                                                                                                                                                                      tGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTg                                       >  1:929195/1‑238 (MQ=255)
                                                                                                                                                                                                         gTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGccc                                    <  1:764315/238‑1 (MQ=255)
                                                                                                                                                                                                          tCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCCCAGACCAATCCATTTGCCCt                                   >  1:411235/1‑238 (MQ=255)
                                                                                                                                                                                                          tCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCCCAGACCAATCCATTTGCCCTg                                  >  1:40459/1‑239 (MQ=255)
                                                                                                                                                                                                          tCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCATGGATGAGCAACTTGCGGCACTCACGGCAATCACGTTGCGCCATCCTAGCTTTCGCCGCCGGAATTGTGGGTGCGCCAATGCTGAGTTCCTCCAGCCCTAATCCGCCCACCCAACGCAGCCCGGAGACTTTCGCTCCCCGCTCACCGCACAACCCAACACATTTAccccg                                  >  2:364905/1‑237 (MQ=255)
                                                                                                                                                                                                                tgGGTTATGGGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCACCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATg                            <  1:307113/239‑1 (MQ=255)
                                                                                                                                                                                                                tgGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGAt                             >  1:69405/1‑238 (MQ=255)
                                                                                                                                                                                                                 ggtgcATGCGGAATTGCGCCGGCGGGGGTTCTACTTCCCACCAAGTGCGGAAGGCCATTGCTTGGGTGAGCAAGTGGCGGCCCTGAGGGCTATCAAGTTGCGCCAACCGAGCCTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGc                           <  2:281569/235‑1 (MQ=255)
                                                                                                                                                                                                                       tGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTa                                                                                                                                                                               >  2:260872/1‑86 (MQ=255)
                                                                                                                                                                                                                       tGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCCATGCACCGCt                     >  1:1219529/1‑239 (MQ=255)
                                                                                                                                                                                                                       tGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGCTTCGCCATCCGCGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGGTCATCCAGCCCTCAGCCGACCCGCAACCGCCGCCCCGCAGCTTTCCGCCCCCCCTCCCCGCCCCACCCAACCCATTTGCCCCCGCGATGCAcc                        >  2:918605/1‑236 (MQ=255)
                                                                                                                                                                                                                         tggACTTGCGCCAGCAGGGGTTCTCCTTCCAGCGAAGGCCGGCAGGCCATTGCCTGGTGTGGCAACCTGCGGCCCTTACGGCTATCAAGTTTCGCCATCCGGGCTTTCGCCGCCGGAATTGATGGTGCGCTCAGGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCGTTCGCCCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTg                   <  2:301431/238‑1 (MQ=255)
                                                                                                                                                                                                                         cGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCACCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCtt                    >  1:380853/1‑238 (MQ=255)
                                                                                                                                                                                                                          ggAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTg                   <  1:1064595/238‑1 (MQ=255)
                                                                                                                                                                                                                          ggAATTGCCCCAGCAGGTGTTCTAATTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTg                   <  1:1007366/238‑1 (MQ=255)
                                                                                                                                                                                                                             atggggcccAGCAGGTGTTCTTCTTCACGCGAAGTACGGCAGGCCAGTGCCTGGTTGAGCAACTTGCGGCACTCGCGGCTATCAAGTGGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTGCa                 <  2:411235/231‑1 (MQ=255)
                                                                                                                                                                                                                             aagtGCGCCAGCATGTGTTCTACTTCCAGGGAAGTACGGCAGGCCATTGCCTGGGTGAGCAAGTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGGGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTGCa                 <  1:559697/235‑1 (MQ=255)
                                                                                                                                                                                                                               ttGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCCCCGCCCCAACCAACCCATTTGCCCTGGCGAGGCCCCCCttcccagg              >  2:785552/1‑233 (MQ=255)
                                                                                                                                                                                                                                  cgcCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGCGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTGCACGGCGt           >  1:1205753/1‑239 (MQ=255)
                                                                                                                                                                                                                                  cgcCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCCCAGACCAATCCATTTGCCCTGGCGATGCACGGCTTGCACGGCGt           >  1:843519/1‑239 (MQ=255)
                                                                                                                                                                                                                                            tgtTCTACTTCCAGCGAAGTACGGCAGGCCTTTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTGCACGGCGTAATCGAgcg  <  1:852057/238‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
ATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCACGGCTATCAAGTTGCGCCATCCGAGCTTTCGCCGCCGGAATTGATGGTGCGCTCATGCTGAGTTCATCCAGCCCTAAGCCGACCAGCAACGGCAGCACGGAACCTTTCGCTCCCAGCTCACCGCACAGACCAATCCATTTGCCCTGGCGATGCACCGCTTGCACGGCGTAATCGAGCG  >  NC_000913/4139965‑4140435

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