Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,103,606 T→C 13.4% S457G (AGT→GGT)  cpxA ← sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,103,6060TC13.4% 304.1 / 3.7 125S457G (AGT→GGT) cpxAsensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base T (48/56);  new base C (5/16);  total (53/72)
Fisher's exact test for biased strand distribution p-value = 8.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TGATGACGATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATAC  >  NC_000913/4103369‑4103839
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
ttagtaccgtCAATTTCCCCGTTTACTCTCCCGTGCCGGGTTATCGAAAAGTGGGACGCGGACGATGCAAAAGCGCCGAATGGGCGGAAAGGTCGAAGATTTTCCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTc                                                                                                                                                                                                                                               <  2:300844/230‑1 (MQ=255)
   agaCGATAAATATCCCCGTTATACCCCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGc                                                                                                                                                                                                                                            <  1:461680/238‑1 (MQ=255)
        aTCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTCCCTGCAAATGCGAAGTTTAACCCCCCTtata                                                                                                                                                                                                                                       >  2:605664/1‑239 (MQ=255)
             tatCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCg                                                                                                                                                                                                                                  <  1:630116/239‑1 (MQ=255)
                   cGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTCCCTGCAAATCCGAAGTTTAACCCCCCTTTTACAACGGCAAcc                                                                                                                                                                                                                            >  2:966878/1‑239 (MQ=255)
                   cGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGAGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAAACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGAAGAGGGGGTAAACGCCTTATCCTGCCTGCAAAACCGAAGTTTACCTCCGCTtataagccgccaaca                                                                                                                                                                                                                            >  2:763031/1‑228 (MQ=255)
                     tttACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGCAGTTTAACTCCGCTTATACAGCGcaaacaaa                                                                                                                                                                                                                         >  2:283075/1‑239 (MQ=255)
                     tttACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAAACGCTATCAGCGTCGCATCAGGCATCCTGCTCAACTGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATCCGAAGTTTAACTCCGCTTTAACAACGGAAACaaa                                                                                                                                                                                                                          >  2:162036/1‑239 (MQ=255)
                        tACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAAACGCTATAAGCGTCGCATCAGGCACCCTGCTCAACTGCCGAAGCCGCGGTAAAGGCCTTATCCTGTCTGCAAATGCGAAGTTTAAATCCGCTTATAACACGCCAACCAAAt                                                                                                                                                                                                                        >  2:701158/1‑239 (MQ=255)
                           ctgccgccgccgGTTTAGCAAAAAGGGGGACGGGGACGGTGAAAAAGCTCCTACTGAGCCGCAAAATCGAAAATTTTTCCCGCCGGTTGGGGGGAAAATAACCCCCGGAGTGGAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCTTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAc                                                                                                                                                                                                                     <  2:207762/236‑1 (MQ=255)
                           ctcccGCGGCGGGGTTTACGAAAAGCGGGGCGCGGACGGTGCAAAAGCGCAGACTGAGCGGCAAGCTCGAAGCTTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGGGTGTAGGCCTGATAAGCCGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAc                                                                                                                                                                                                                     <  2:909085/236‑1 (MQ=255)
                           ccgccGCTGCCGGGTTATCGAAACGCGGGACGCGGACGGTGCAAAAGCGCCGCACGAGCCGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAc                                                                                                                                                                                                                     <  1:963771/236‑1 (MQ=39)
                              ccGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTTTCCTGCCTGCAAATGCGAAGTTTAACCCCGCTTATAAAACGGGAACaaaaaaccca                                                                                                                                                                                                                   >  2:345014/1‑232 (MQ=255)
                              ccGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGGGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCCTCCTGCTCAAATGCCGGATCGCGCGTAAACGCCTTATCCTCCTTCCAAATGCGAAGTTTAACCCCCCTTTTAAACAGGCAACaaaacaccca                                                                                                                                                                                                                   >  2:479349/1‑232 (MQ=255)
                                  tGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGTGGGTAGGCCTGCTCAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCCAATGCCAGATGCGGCGACAAGGCCCTATTCTGCCTGCACATGCGCAGTTTAACTCCGCTTAACAAGGGCCAACCaaaaaaccccccgg                                                                                                                                                                                                              >  2:694303/1‑228 (MQ=255)
                                      ggtttATCGAAAATGTGGCCCCGGGCGTTGAAAGGGGCCGACTGAGCGGCAGAGTGGAAGATTTTCCCCGCCGGTTGTGGGGAAAAGAACCCCCGGGGGGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAAcc                                                                                                                                                                                                         <  2:484643/236‑1 (MQ=37)
                                            tCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATAAAGCGGCAACCAAAACACCACCCGAAAACCGCCCa                                                                                                                                                                                                    >  2:1082443/1‑239 (MQ=255)
                                            tCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCCCTTATACAGCGGCAACCAAATCACCAGCCGCAAACCGccc                                                                                                                                                                                                     >  1:98258/1‑238 (MQ=255)
                                            ccaaaagGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCa                                                                                                                                                                                                    <  1:728662/236‑1 (MQ=255)
                                                 aggtgggACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGTGCGAAGATTTTTCCCGCCGGTTGGGGGGAAAATAACCCCCGGAGTTTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                                                                                                                                                                                               <  1:902504/234‑1 (MQ=255)
                                                 aaGCGGGGCGGGGACGAGCAAAAAGCCCCGTCTGAGGGCCACGACCGAGATTTTTCCCCGCCGGTTGGGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAGGCCGCTATCAGCTGCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                                                                                                                                                                                               <  2:323150/239‑1 (MQ=255)
                                                  acctgaccgcgGACAAGGAAAAAGGCCCAACGGGGGGGCAAGATCGAAGATTTTCCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGc                                                                                                                                                                                              <  1:193309/232‑1 (MQ=255)
                                                  aGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGc                                                                                                                                                                                              <  1:345014/239‑1 (MQ=255)
                                                  aGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAAATCCCCTTATACAGCGGCAACCAAAACACCAACCGTAAACCGCCCACCGGGc                                                                                                                                                                                              >  2:910958/1‑239 (MQ=255)
                                                      ggACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCgg                                                                                                                                                                                                >  1:323150/1‑233 (MQ=255)
                                                         ccggggcggtgaaaagGCGCCGACTGGGGGGCAAGTTCGAAGTTTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtctt                                                                                                                                                                                        <  2:111312/226‑1 (MQ=39)
                                                              aCGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTTGCTTCGGcc                                                                                                                                                                                   >  1:936612/1‑238 (MQ=255)
                                                                          cgcCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGGTTAACTCCGCTTATACAGCGGCAACCAAATCACCACCCGTAAACCGCCCACCGGGCTGGCTCCTGCCTTCACCCAAccc                                                                                                                                                                       >  1:840940/1‑237 (MQ=255)
                                                                           gcgggcTGGGCGGCAGGAGCGAAGTTTTTTCCCGCCGGTTGGGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCGGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCa                                                                                                                                                                       <  2:261463/232‑1 (MQ=255)
                                                                            cgcgtggAGCGGCAATATCGGAGTTTTTTCCAGCCGGTTGGGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGCCGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTACCCTGCCTGCAAATGCGAAGTTTAACTCCGCTGATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                                                                                    <  2:411536/233‑1 (MQ=255)
                                                                            ccGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTACACCCCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                                                                                    >  1:1054681/1‑239 (MQ=255)
                                                                            ccGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACg                                                                                                                                                                     >  1:261463/1‑238 (MQ=255)
                                                                            ccGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAACGGCAACCAAATAACCAGCCGTAAACAGCCCACCGGGGTGTCTTCTTCCTTAACCCAACCAcc                                                                                                                                                                     >  2:685373/1‑237 (MQ=39)
                                                                            ccGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGACTTATCCTGCCTGCACAGGCGAAGTTTAACTCCGCTTATAACGCCGCCACCAAATCACCCGCCGTAAACCGCACACCGGCCTGTCTTTTGCCTTccccacaccaca                                                                                                                                                                     >  2:812512/1‑230 (MQ=255)
                                                                             ggACTGAGGGGCAAATGCGAAGATTTTCCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                                                                                    <  1:334844/237‑1 (MQ=255)
                                                                             cGACTGGGCGGCAATGTCGAAGTTTTTTCCCGCCGGTTGGGGGGAAAATAACCCCCGGAGGGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                                                                                    <  1:414982/238‑1 (MQ=255)
                                                                                 taggggacAAGACCGAAGTTTTTGCCCGCCGGTTGGGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGAtgctg                                                                                                                                                               <  1:694303/232‑1 (MQ=255)
                                                                                                    tttttCCGGCGGGTGGGGGGGAAAAAAAGCCCCGGGGTGGAGGCCGGATACGGCCGTAAGGCCGTTGGATCAGGCACGCAGCTAAAATGCCGGAGCCGGGTTAACCGCCTTACCCTGCCTGCAAATGCGAAGTGTAACCCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGGGGGCTGTCGTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTcaaca                                                                                                                                           <  2:676579/239‑1 (MQ=255)
                                                                                                    tttttCCGGCCGGTTGGGGGGAAAACTACCCCCGGAGTGTAGGCATGAAAAGCCGCTATCATCGTCGCATCAGGCTTCCTGCTCAAATGCCGGATGCGGCGAAACCGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTcaaca                                                                                                                                           <  1:1072034/239‑1 (MQ=255)
                                                                                                      tttCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCAAGCCACGATGCTGCTGAATGGCGGGTTtaaaaa                                                                                                                                          >  2:59780/1‑238 (MQ=255)
                                                                                                       ttCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGCAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCACCGGCCGGTCTTCTGCCTTCCCCCACCCACAATGTGGCTGAAGGGGGGGTTTcaaaac                                                                                                                                          >  2:673477/1‑235 (MQ=38)
                                                                                                           ccgccgGTTGTGGGAAAAATAACCCCCGGAGGTGAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGcc                                                                                                                                     <  1:81751/239‑1 (MQ=255)
                                                                                                                ggTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCACCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTTACCAAACCCCAGcc                                                                                                                                >  1:616874/1‑238 (MQ=255)
                                                                                                                ggTTGGGGGAAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGccc                                                                                                                               <  1:599506/239‑1 (MQ=255)
                                                                                                                            aaTAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGGTTTAAAAAACCCCAGCCCCAAAcccttcc                                                                                                                   >  2:209414/1‑233 (MQ=255)
                                                                                                                             ttaccccccGGAGTGTGGGCTTGAAACGGCGTTTTCGGGGTCGCATCAGGCTTCCTGCTAAAGTGCGGGATGCGGCGTAACCCCTTTACCCTGCCTGAAATTGCCAAGTTTAACCCCCTTTAGACAGCGGCAACCAAACCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCg                                                                                                                <  2:1045245/235‑1 (MQ=255)
                                                                                                                               aaCCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGcc                                                                                                              <  1:98604/239‑1 (MQ=255)
                                                                                                                               aaCCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAAGGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGcc                                                                                                              <  1:605664/239‑1 (MQ=255)
                                                                                                                               aaCCCCCGGAGTGTAGGCCTGATAAGACGCTACCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGcc                                                                                                              <  1:528877/239‑1 (MQ=255)
                                                                                                                               aaCCCCCGGAGTGCAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAACCGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGcc                                                                                                              <  1:222692/239‑1 (MQ=255)
                                                                                                                               aaCCCCCGGAGGGTAGGCCTGATAAGACGCTAGCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGcc                                                                                                              <  1:77139/239‑1 (MQ=255)
                                                                                                                                      ggtgtgtAGGCATGATAAGACGCTAACAGCGTCCCCGCAGGCCTCCTCCTCAAATGCCGGATGCGGCGTAAACGCCTTACCCTGCCTGCAAATGCGAATTTTAACCCCGTTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCa                                                                                                       <  2:713972/236‑1 (MQ=255)
                                                                                                                                       gAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTTAACAAACGCCCGCCCCAAACCTGTACCCCCAGATTCAc                                                                                                      >  1:389439/1‑239 (MQ=255)
                                                                                                                                          tgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                                                                                                                                                                                               <  1:1035183/151‑1 (MQ=255)
                                                                                                                                                 cTGATAAGCCGCTATCGGCGTCGCATCAGGCATCCTGCCCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACCCCGCTTATCCAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTg                                                                                             <  1:359713/238‑1 (MQ=255)
                                                                                                                                                  tGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGGTTCAACAAACGCCAGCCCCAAACCCGTCACGCCCGGTTACCGCTcgcgggtt                                                                                           >  1:26268/1‑235 (MQ=255)
                                                                                                                                                    aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGACTGGCGGGTTCAAAAAACGCCACCCCCCAAACACGACCGCCCCATTCACGATccgcgcgct                                                                                          >  2:38474/1‑230 (MQ=255)
                                                                                                                                                    aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCACCAATCCCCCGCCCCCAACCCGCACCCCCAAATTCACGATccccagctcc                                                                                         >  2:83261/1‑230 (MQ=255)
                                                                                                                                                    aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCCAACCTGTAACCCCAGATTAAccactccgtgct                                                                                           >  1:630377/1‑226 (MQ=255)
                                                                                                                                                    aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCCAAACTGTAACCCCAGATTCACGAACCCGGGCttt                                                                                         >  1:685308/1‑238 (MQ=255)
                                                                                                                                                     aaaGACGTTATCACGGCCGCATCAGGCATCCTGCCCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTc                                                                                         <  1:389254/237‑1 (MQ=255)
                                                                                                                                                      aaGACGCTAACAGCGTCGCATCAGGAATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCAt                                                                                       <  1:703986/239‑1 (MQ=255)
                                                                                                                                                         aCGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCCGCCCCAAACCCGTAACCGCCAATTCACGATGCCGGGGTTCATCgg                                                                                    >  1:644749/1‑239 (MQ=255)
                                                                                                                                                         aCGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCGGCCTTAACCCAGCCACGCTGCTGCTGAATGGGGGTTTCCACAACCCCCAGCCCCCAACCCGTCCCCCCAAGTTTACCATcgcgggtctccccgg                                                                                    >  2:985940/1‑228 (MQ=255)
                                                                                                                                                          agCTATCAGCTTCGCAGTAGGCATCCTGCGCAAAGGCCGGAAGGGGCCTAAACGCCTTATCCTGCCTGCAAAGGCGAAGTTTAACTCCGCTTCTACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCgg                                                                                    <  1:1071091/237‑1 (MQ=255)
                                                                                                                                                              aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                                                                                                                                                                                               >  2:436969/1‑131 (MQ=255)
                                                                                                                                                              aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                                                                                                                                                                                               <  1:436969/131‑1 (MQ=255)
                                                                                                                                                              aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGGTTCAACAATCGCCAGCCCCACACCTCTAACGACACCTTCACGATCGCGCGGTTccactggagga                                                                               >  2:83843/1‑229 (MQ=255)
                                                                                                                                                                      cGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCCAACCTGTAACGCCCGATTCACGATCGCGGGCTTCATCGGTCCGAAAGCacgg                                                                       >  1:830334/1‑239 (MQ=255)
                                                                                                                                                                          tCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGGGCTTCATCGGGACGAGAGAacggacgg                                                                   >  1:153535/1‑239 (MQ=255)
                                                                                                                                                                                  ccTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAACCCACACACCCAAACCCTGTCCCCCACGATCACGACACGGTGCTTCATCAGTACGAAAAAAAGCACGGaaaa                                                               >  2:404018/1‑235 (MQ=255)
                                                                                                                                                                                  ccTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGGTTCAACAATTCCCAGCCCCCAAACCTTACCCCCCCACTCACCATACCGTGGTTCTTCGGGACgcaaaaaacgccgagaaaacta                                                           >  2:935505/1‑218 (MQ=255)
                                                                                                                                                                                   cTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTTACCAAACCCCCCCCCCCACCCCGTACCGCCCCATTCACTACCCCCTGCTTCATCGGTCCAATAAAACGGACCGAAAAACttt                                                          >  1:906047/1‑237 (MQ=255)
                                                                                                                                                                                   cTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCCGATTCACGATCGCCTGCTTCATCGGAACGATAGAACGGACGGGAAATCTGt                                                          >  1:688921/1‑239 (MQ=255)
                                                                                                                                                                                   cTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAACTGCCCCCCCCAAACCCGTACCCCCCAAATCCCAAGCGCGCGTTTCACCGGGGCGCTAAAACAGAAGGaaaaaactt                                                          >  2:730696/1‑234 (MQ=255)
                                                                                                                                                                                   cTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACCCCGCTTATACAGCGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGt                                                          <  1:562843/239‑1 (MQ=255)
                                                                                                                                                                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCCGAACCGCCAGATTCACCATCGCGTCCTTCCTCGGTCCGAACGACCCGACGGGAAACCtttt                                                         >  1:881376/1‑238 (MQ=255)
                                                                                                                                                                                              ccggtggcgcgtaaaaCCCCTTACCCTTCCGGCAAAGGGGAAGTTTACCCCCGCTTTAACCGCGGACACCAAACCACCAGCCGTAACCCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTc                                              <  2:936612/226‑1 (MQ=255)
                                                                                                                                                                                               cggttccgcggTAAACCCCTTCCCCTCCCGCAAAATGCAAAGTTTAACCCCGCTTATACAGCGGCAACCAAATCCCCAGCCTTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTcc                                             <  2:1054681/229‑1 (MQ=255)
                                                                                                                                                                                                ggATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCa                                                                                                                             <  2:207327/161‑1 (MQ=255)
                                                                                                                                                                                                ggATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCa                                                                                                                             >  1:207327/1‑161 (MQ=255)
                                                                                                                                                                                                    gcggggTAAACGCCCTATCCTGCCGGCAAATGCGAAGTTTAACTCCGTCTAGACGGCGGCAACAAAATCACCAGCCGTAAACCGCCAAGCGGGCTGTCTTCGGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCt                                        <  1:870707/239‑1 (MQ=255)
                                                                                                                                                                                                     cggcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTGAACTCCGCTTATACAGCGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCCCCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCt                                        <  1:468741/238‑1 (MQ=255)
                                                                                                                                                                                                               gCCTTACCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATCCAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGAc                             <  1:730696/239‑1 (MQ=255)
                                                                                                                                                                                                                  ttACCCTGCCTGAAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCAt                          <  1:874192/239‑1 (MQ=255)
                                                                                                                                                                                                                    tcccttcctgCAAATGCGAAGTTTACCTCCCCTGATACAGCGGCAACACAAGCACCAGCCGAAAACCGCCAGGGGGGCTGTTTTCTGCCTGCACCCCGCCACGATGATGCGGAATGGCGGGTTCAACAATCGCCAGCCCCAAACCTGGACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGCCCAtcg                        <  2:916578/237‑1 (MQ=255)
                                                                                                                                                                                                                    aTCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCCACAATCTCCCGCCCCCAACCTATACCGCCAGCTTAACAGTCGCGGGCTTCATCCGGACGATCTAACAGGCGGAAAAACTTGTCCCGATCTTACCGGGCAACCCCAGGACCAacg                        >  2:892384/1‑239 (MQ=255)
                                                                                                                                                                                                                         gtcttaaattgggAAGTTATACCCCGCTTTTAGGGCGGCACACAAAACGCCAGCCGTAATCCGCCCAGCGGGCTTTCTTCTGCCTCCCACCAGCCACGATGCTGCTGGAGGGCGGTTTCAACACTCGCCTGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTTCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCgtcgtc                   <  2:612888/227‑1 (MQ=255)
                                                                                                                                                                                                                         gccttcAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCCTCCGTCACCCAGCCACGATGCTGCTGAATGGCGGTTGCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCgtcgtc                   <  1:529189/239‑1 (MQ=255)
                                                                                                                                                                                                                         gccttcAAATGCGAAGTTTAACTCCGCTTAAACCGGGGCAACCAAATCCCCAGCCGTAACCCGCCCAGCGGGCTGTCCTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCgtcgtc                   <  1:200946/239‑1 (MQ=255)
                                                                                                                                                                                                                         gcctccAAATGCGGATTTTACCCCCCGTTTAACAGGGCCCACCAAATCCCCGGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTCCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCgtcgtc                   <  2:644749/239‑1 (MQ=255)
                                                                                                                                                                                                                         gccgtcAAATGCGAAGTTTAACCCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTTTTCTGCCTTCCCCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCgtcgtc                   <  1:989264/234‑1 (MQ=255)
                                                                                                                                                                                                                           ctgcAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCCAACCTGTACCGCCAAATTCACGATCCCGGGCTTCATCGGTACGATAGAACGGACGGAAAAACTGTTCCCGATCTTCCGGGCTAACGCCCGGACCacccccgttc                  >  1:258065/1‑229 (MQ=255)
                                                                                                                                                                                                                           cctcAAATGCCAAGTTTAACCCCTCTTATACAGCGGCAACCAAACCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCCTCCCCCAGCCACGATGCGGCTGAATGGCGGTTTAAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCa                 <  2:534260/236‑1 (MQ=255)
                                                                                                                                                                                                                               aaatcgaaaGTTTAACTCCGTTTAGACAGGGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGt             <  1:83261/232‑1 (MQ=255)
                                                                                                                                                                                                                               aaaggggAAGTTTAACCCCGCTTACACAGGGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTTTTCTTCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGt             <  2:616874/233‑1 (MQ=255)
                                                                                                                                                                                                                               aaaTGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGt             <  1:59780/239‑1 (MQ=255)
                                                                                                                                                                                                                               aaaTGCGAAGTTTAACCCCGCTTATACAGGGGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGt             <  1:910958/239‑1 (MQ=255)
                                                                                                                                                                                                                                          ttAACTCCGCTTATACAGCGGCAACCAAACCCCCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATa   <  1:673477/238‑1 (MQ=255)
                                                                                                                                                                                                                                          ttAACCCCGTTTATACAGGCGCAACCAAATCCCCAGCCGTAAACCGCCCAGCGGGCTGGCCTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATAc  <  1:573306/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
TGATGACGATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGATCTTCCGGGCTAACGCCAGGACCATCGTCGTCCACCGTAATGGTGATAC  >  NC_000913/4103369‑4103839

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