Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,530,175 T→A 12.4% I24F (ATC→TTC)  sgcC ← putative PTS system EIIC permease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,1750TA12.4% 244.6 / 3.0 88I24F (ATC→TTC) sgcCputative PTS system EIIC permease component
Reads supporting (aligned to +/- strand):  ref base T (51/24);  new base A (2/10);  total (53/35)
Fisher's exact test for biased strand distribution p-value = 1.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAT  >  NC_000913/4529937‑4530408
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           
ttGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGCCAAAACCAATCCCAACCGTTACCCCTGCCTTAATAGCCTGCAGCCAGGGGATGCGGAAGAt                                                                                                                                                                                                                                             >  2:812410/1‑239 (MQ=255)
 tGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGTTa                                                                                                                                                                                                                                            >  2:240757/1‑239 (MQ=255)
        aaaaaTAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCATAGACATGCAGGTCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGCTCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATACGGAAGCTAAGGCCa                                                                                                                                                                                                                                      >  2:774049/1‑238 (MQ=255)
        aaaaaTAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATAACGATCACCAGTCCCCAACCGCACAAACCAATCCCAACCGTTCACCCAGCCTTAATCCGCCGAACCCAGGGGATGGCGAAAGAAAAgcaa                                                                                                                                                                                                                                       >  1:81315/1‑235 (MQ=255)
           aagatcTGGGATAATGTTGGCGCCAAGGGGTGGCGAATAACCAACGCCGGAGGCAGGGCCCGCACCAACGGCAAAGACATGCAGGGTAAGGCCGAAACGCTCAATCAACACTTTGATAGGGGGACTGAGGCTTTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAAGCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCa                                                                                                                                                                                                                                  <  1:812556/234‑1 (MQ=255)
           aaTGACTGGGATAATCATGGCGCCAATGGCTTTCGCATAACCAAGGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCa                                                                                                                                                                                                                                  <  2:18000/239‑1 (MQ=255)
             ataagggaaaggccagggCGCCAAGGGCTGTCCCATAACCACCGGCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCCCGTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACg                                                                                                                                                                                                                                <  1:669267/223‑1 (MQ=255)
              accggggATAATCATGGGCCCAATGGCTGACCCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGGCAAAGGCATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACg                                                                                                                                                                                                                                <  1:420537/234‑1 (MQ=255)
               aCTGGGATAATCAGGGCACCAATGGCTGGCGCAGAACCAACGCCGGAGGCAGGGCCTGCGCCAGCGTCAAAGACAGGCAGGGCTAGGCCGAAACGCTCAATCATCGCTGTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCTCCAGTCCCATACCGACATATCCAATCCCAACCGTTTCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGAt                                                                                                                                                                                                                              <  2:300557/239‑1 (MQ=255)
                  gggATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGAATGAGGCTGTCGATCGCCCTTACGATACCCAGTCCCATCCCGACAAAACCAACCCCAACCGTTACCCCAGCCTTATTCGCCTCCCACACGGGGTATCGGAAGATAAGGcacaaccacaacaa                                                                                                                                                                                                                           >  1:192511/1‑226 (MQ=255)
                    gATAACCATGGCCCCAATGCGTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGGCAAAGACATGCAGGTGTAGGCCGAAACGCTCAAGCAGCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCatg                                                                                                                                                                                                                          <  1:293301/238‑1 (MQ=255)
                       aaTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCACTCACGATCCtgatga                                                                                                                                                                                                                     >  2:420537/1‑239 (MQ=255)
                         tCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATa                                                                                                                                                                                                                   <  2:72296/239‑1 (MQ=255)
                             ggCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCCTGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCCCGATCATGTTGATAggg                                                                                                                                                                                                                >  2:783113/1‑238 (MQ=255)
                             ggCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCCCGATCATGCTGATAggg                                                                                                                                                                                                                >  2:345936/1‑238 (MQ=255)
                             ggCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCCTGATGATAggg                                                                                                                                                                                                                >  2:1011026/1‑238 (MQ=255)
                              gcgcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAgg                                                                                                                                                                                                                                         >  2:789968/1‑213 (MQ=255)
                                   aagggCTGGCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTAAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGaaaa                                                                                                                                                                                                         <  2:183175/236‑1 (MQ=255)
                                   aaTGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCCATCCCATACCGCCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGTTAAGGCCAATCACGATCAtgatga                                                                                                                                                                                                                     >  2:669267/1‑227 (MQ=255)
                                         cgtccatcaaCCCACCCGGGGGGCGGGGCCTGCCCCAAGGCAAAAGGCATGGAGGGTTGGGCCAAACCCCCCAATCACCCTTTTGATAGGGGGACTGGGGTTGCCAACCGCCATGCCGATCCCCAGCCCAATCCCGAAAAAACCAATCCAAACCTTTCCCCCTGCCTTAACCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTAc                                                                                                                                                                                                  <  1:977775/231‑1 (MQ=255)
                                                 aaCCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCcaga                                                                                                                                                                                          >  2:769775/1‑239 (MQ=255)
                                                       cGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTTCCGCCAGGGGATGCGGAAGATAAGGCCAACCACGATCATAATGATAGGGCCGAAACCGGTACCCCCGGAAGACaag                                                                                                                                                                                    >  1:426738/1‑237 (MQ=255)
                                                             aGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCTATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCACCCAGGGGAGACGGCAGTAAAGACCAAACACGATCATGGTTATAGGGCCGAAACCAGTACCGCCAGAAGACAgaggataa                                                                                                                                                                              >  1:877247/1‑232 (MQ=255)
                                                                       tGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCCCGATCATGATGATAGGGACGAAAACGGGACCGCCGAGAGACAGGATGTAATCAAACATaa                                                                                                                                                                    >  2:886580/1‑239 (MQ=255)
                                                                          cccAACGTAAAAGACATGCAGGGCTGGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCAATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGc                                                                                                                                                                 <  2:532611/238‑1 (MQ=255)
                                                                             aaCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCCCGATCATGTTGCTAGGGACGAAAACGGAACCGCCGGGAGACAGGAGGTAAACAAAAAAAAAGccc                                                                                                                                                               >  1:778162/1‑237 (MQ=255)
                                                                             aaCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATAACGATCATGATGTTAGGGACGAAAACGGTACCGCCGGGTGACAGGATGTAATAAAACATAATGcct                                                                                                                                                               >  2:890005/1‑238 (MQ=255)
                                                                                    aaGACATGAGGGGCTAGGCCGAAAGGCTCAATCATCACTTTGATAGGGGGACGGAGGCTGTCGATCGCAATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTAt                                                                                                                                                       <  2:767379/239‑1 (MQ=255)
                                                                                    aaGACAGGCAGGGCTAGGCCGAACCGCTAAACCATCCCTGTGATAGGGGGACTGAGGCTGCCGATCGCCATGACGATCACCAGTCCCATACCGACAAACCCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTAt                                                                                                                                                       <  2:426738/239‑1 (MQ=255)
                                                                                           gCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCAc                                                                                                                                                                                                                                 >  2:122424/1‑160 (MQ=255)
                                                                                           gCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCAc                                                                                                                                                                                                                                 <  1:122424/160‑1 (MQ=255)
                                                                                                 cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATCAGGCCAATCCCGATCATGTGGATAGGGCCGAAAACGGTACCCCCGCGAGACAGGAGGTAATCAAACATAATGCCACCTTTTATTGTGTTaaacag                                                                                                                                           >  1:198043/1‑234 (MQ=255)
                                                                                                    ggCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGAAATCAAAAATAAAGCCACCCTTTATTGAGTTaaaaaag                                                                                                                                          >  1:471078/1‑231 (MQ=255)
                                                                                                    ggCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATCATGATAGGGACGAAAACGGAACCGCCGAGAGCAAGGATGTAATCAAACATAATACCTCCTTTTATTCAGTTAACAAAGCCt                                                                                                                                       >  2:345946/1‑239 (MQ=255)
                                                                                                         aaaCGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGACGATAGGGACGAAAACGGTACCGCCGAGCGACAGGATGAAATCAAACATAAAGCCTCCCTTTATTGAGTTAACAGAGGCttagtt                                                                                                                                  >  1:716353/1‑235 (MQ=255)
                                                                                                         aaaCGCTCAATCATCACTTTGATAGCGGGACTGCGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGAAGCCAGGGCATGCGGAAGATAAGGCCAATACTGATCATGCTGCTAGGGCAGAAAACGGTTACGCCGAGGGACAGGATGTAATCAACACTAATGCCTCCTTTTATTGAGTGAACAGACACTTAAtt                                                                                                                                  >  1:237467/1‑239 (MQ=255)
                                                                                                           ccGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGACCCCCAGTCCCATACCGAAAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGACGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATttg                                                                                                                                <  2:716353/238‑1 (MQ=255)
                                                                                                                       caCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAAcgc                                                                                                                    >  2:281521/1‑239 (MQ=255)
                                                                                                                             gATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAAt                                                                                                                                                                                                                                    <  2:760311/123‑1 (MQ=255)
                                                                                                                             gATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAAt                                                                                                                                                                                                                                    >  1:760311/1‑123 (MQ=255)
                                                                                                                                  ggggACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAAAGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTg                                                                                                          >  2:97156/1‑238 (MQ=255)
                                                                                                                                   ggaaCTGAGGCTTTCGATCGCAAGAACGATCCCCAGTCCCATACCGACAAACCCAATCCCAACCGTTCCCCTTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGa                                                                                                         <  2:377599/235‑1 (MQ=255)
                                                                                                                                    ggACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACa                                                                                                                                             >  2:57142/1‑201 (MQ=255)
                                                                                                                                    ggACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGAt                                                                                                        >  2:861258/1‑238 (MQ=255)
                                                                                                                                    ggACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGCAATCAAACATAATGACTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAATGCGTCGTCGTTGAt                                                                                                        >  2:702973/1‑238 (MQ=255)
                                                                                                                                      aCTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGAACCGCCGAGAGACAGGATGTAATCAACCATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGCTTTGTTGTTTTAACGGCTCGCCGTTGATc                                                                                                       >  2:485688/1‑237 (MQ=255)
                                                                                                                                        tGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCCTTTATTGAGTTAAAAGAGCCTTGATTTTTTGTTTTAACGCGGCGTCGTTGATACCtt                                                                                                   >  1:377599/1‑238 (MQ=255)
                                                                                                                                          aGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGAATGCGGAAAATAAGGCCAATCACGATCAGGTGGAAAGGGACGAAAAGGGTACCGCGGCGACACAGACTGTAATCAAACAAAATGCCCCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTTTTTTAACGCGTCGCCGTTAATCCCtgtg                                                                                                 >  1:269325/1‑239 (MQ=255)
                                                                                                                                            gCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGAt                                                                                                                                                                                                                              <  1:262637/114‑1 (MQ=255)
                                                                                                                                            gCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGAt                                                                                                                                                                                                                              >  2:262637/1‑114 (MQ=255)
                                                                                                                                            gCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAAAAGAGCCTTGATTTGTTGTTTTAAAGCCTCTTCGTTTATCCCtgtggc                                                                                               >  1:640630/1‑237 (MQ=255)
                                                                                                                                               gTCGATCGCCATGACGATCACCAGTCCCATACCGGCAAAACCAATCCCAACCTTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATAATGATGATAGGGACGAAAACGGTACCGCCgagaga                                                                                                                                                                                        <  2:738330/147‑1 (MQ=255)
                                                                                                                                               gTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCgagaga                                                                                                                                                                                        >  1:738330/1‑147 (MQ=255)
                                                                                                                                                          tGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACAATCATGATGATAGGGACGAAAACGGTACCGCCGAGGGACCGGACGTAAAAAAACATAATGCATCCTTTTATTGAGTTAAAAGAGCCTTGATTTGGTGTTTTAACGCGTAGTGGTTGATCCCTGTGAGCAGGGGGGGGCCAt                                                                                 >  1:451233/1‑239 (MQ=255)
                                                                                                                                                           gccgccccccACTCCCAGACCGACAAAACCCAGCCCACACGTCCCCCCTGCCTTGCTCGCCGGCAGCCGGGGGATGGGGAAGAAAAGGCTAATCACGTGCATGAGTCTAGTAACGAAACCGGGACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCAtt                                                                                <  1:1001681/231‑1 (MQ=255)
                                                                                                                                                           aaaagggcgttgTCCCATGCCGAAAAACACATGCCCCAGCGTCGCCCACGCCATAATCCTCTGAAACAAGGGGGAGCTGACGACACGGCCAACCACGAGCATGATGGTAGGGACGAAGAGGGGACCGGAGAGAGACAGGATGTAATCAAACATAATGCCGCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCAtt                                                                                <  1:201770/228‑1 (MQ=255)
                                                                                                                                                                 ccccAGTCCCATCCCGAAAAACCCAACCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAggg                                                                                          <  2:979805/221‑1 (MQ=255)
                                                                                                                                                                      gTCCCATCCCGAAAAAACCAATCCCAACCGTTACCCCTCTCTTAATCGCCTGCGGCCCGGGGATGGGGAAGATAAGGCCAATCCCGATAATGATGGTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGa                                                                     <  2:269325/239‑1 (MQ=255)
                                                                                                                                                                      gTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACTCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGa                                                                     >  2:1017982/1‑239 (MQ=255)
                                                                                                                                                                      gTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAGGCGTAGGGa                                                                     >  2:672221/1‑239 (MQ=255)
                                                                                                                                                                         cgtctccaaaaaagCCGATCCCAACCTGTACCCGTGCTAGAAGCGCCTAGGGTCAGGGTATGCGGAAGAAAAGGCAAATCACGTTCATGATTATAGGGAGAAAAACGGTCCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGAGTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGACGc                                                                  <  1:801701/229‑1 (MQ=255)
                                                                                                                                                                          cATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGcc                                                                 >  2:687319/1‑239 (MQ=255)
                                                                                                                                                                           aTACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGCCGTCGTTGTCCCCTGTGGGCAGGGTGGCGCCATTacttgcccttaaaaac                                                                >  1:979805/1‑224 (MQ=255)
                                                                                                                                                                               ccaaaaaccccacccaccccGTGACCCCGTCCTTAATCGCCTGCAGCCAGGGGATGCAGAAGAAAATGTCAATCAGGATCATGATGATAGGGACGAAACGGGTCCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAt                                                            <  1:769775/221‑1 (MQ=255)
                                                                                                                                                                               caaaaaaaCCAATCCAAACCGTTCCCCCTGCCTTAACCCCCTGCAGCCGGGGAATGGGAAATATAAGGCAAATCAAGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAt                                                            <  2:876335/235‑1 (MQ=255)
                                                                                                                                                                                 aCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTAtt                                                                                                                                                      >  2:958426/1‑147 (MQ=255)
                                                                                                                                                                                 aCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTTTTTTTTTAACGCGTCGTGGTTTTTCCCTGTGAGCAGGGCGGGGCCATTAAGCGTAGGGTGGCCGAAAACa                                                          >  1:971229/1‑239 (MQ=255)
                                                                                                                                                                                 aCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATc                                                           >  2:369621/1‑238 (MQ=255)
                                                                                                                                                                                 aCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATc                                                           >  2:207851/1‑238 (MQ=255)
                                                                                                                                                                                 aCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGCCGGCGCCATTAAGCGTAggg                                                                      >  2:201770/1‑227 (MQ=255)
                                                                                                                                                                                         aaTCCCAACCGTTACCCCGGCCTAAAGCCCCGGCAGCCAGGGGATGCGAAAGAAAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGt                                                  <  2:442391/239‑1 (MQ=255)
                                                                                                                                                                                                ccctttcccccTGTCTTGACCGTCGGAAGCCAGGGAAGGCGGAAAGAAAGGACAACCGCGATCAAGATAAAAGGAACGAAACGGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATg                                           <  1:281521/232‑1 (MQ=255)
                                                                                                                                                                                                aCCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGCTTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGGAGGGATGCCGTAATAACTGTTGGGTCGCATg                                           >  1:24873/1‑239 (MQ=255)
                                                                                                                                                                                                  ctttcccctTGCGTTATTCGCGTGAGCCCGGGGATGGGGAAAGCAAAGGCCAATCAGGATCATGATTAGGGGGACGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGa                                         <  1:369621/234‑1 (MQ=255)
                                                                                                                                                                                                     tACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTGGCGTGGAGGt                                      >  2:267487/1‑239 (MQ=255)
                                                                                                                                                                                                      aCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGtt                                     >  2:912644/1‑239 (MQ=255)
                                                                                                                                                                                                           tGCCTAAACCGCCTGCAGCCGGGGGAGGGGGAAGAAAAGGCAAATCACGATCAGGATGATAGGGACGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAt                                <  2:198043/239‑1 (MQ=255)
                                                                                                                                                                                                           ctccTTAAGCGCTGGGAGCGGGGGGATGCGGAAAAAAAGGCAAACCAGGATATGTATGAAAGGGACAAAAAGGGTACCGACGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAt                                <  1:861258/237‑1 (MQ=255)
                                                                                                                                                                                                            gCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGTCGACGCCATTAAGCGTAGGGATGCCGCAATCCCTGTTGGTTCGCGTGGGGGTTACaaaa                               >  1:675648/1‑237 (MQ=255)
                                                                                                                                                                                                            gCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGGGGTTACAATa                               >  2:751203/1‑239 (MQ=255)
                                                                                                                                                                                                            gCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTAATCCCTGTAAGCAGGGCGGCGCCATTAAGCGTAGGGATGCGGAAATCCCTGTTGGTTCCCATGGGGGtaaaaaaa                               >  1:827538/1‑232 (MQ=255)
                                                                                                                                                                                                            gCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAACCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATa                               >  2:47482/1‑239 (MQ=255)
                                                                                                                                                                                                                      cctccggCCTGGGGATGGGGAAAAAAAGGCCAAACACGAACATTATGGTAGGGACGAAAACGCTACCGCCGAGAGACAGGATGTAAAAAAACATAATGCCTCCTTTTATTGAGTCAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATAc                     <  1:1017982/233‑1 (MQ=255)
                                                                                                                                                                                                                       ctgaagcagggggAGGCGGAAGATAAGGCCAATCGCGATAATGTTGATAGGGACAAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATAc                     <  2:192511/229‑1 (MQ=255)
                                                                                                                                                                                                                       cTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAGCAGAGCCTTGATTTGGTGTTTTAGCGCGTCGTCGTTGATCCCTGTGAGCAGGGGGGCGCCATTAAGCGTAGGGATGCCGTAATCCCTTTTGGTTCGCATGGGGGTTAAAAAAAAGACCATAcc                    >  1:270376/1‑239 (MQ=255)
                                                                                                                                                                                                                                            aataTAAGGACAATCCCGATCATGATGAGAGGGACAAAACGGGACCCGCCTGGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGATTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt  <  2:451233/235‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGAAAAGGCCAATCCCGATAAGGATAAAGGGGAAGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt  <  2:471078/238‑1 (MQ=255)
                                                                                                                                                                                                                                              gATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGa                                                                                                         >  1:336047/1‑133 (MQ=255)
                                                                                                                                                                                                                                              gATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGa                                                                                                         <  2:336047/133‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           
TTGAGAAGAAAAATAACTGGGATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAT  >  NC_000913/4529937‑4530408

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