Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,980,009 A→G *86* (TAA→TAG yjhX → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,980,0090AG99.7% 2258.7 / ‑3.8 584*86* (TAA→TAGyjhXhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (272/307);  total (276/308)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 588 total aligned reads displayed.

CAGCCGCGAAGGGCTGTTGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAAGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATGAAAGCGCGC  >  DH1/3979859‑3980156
                                                                                                                                                      |                                                                                                                                                   
  gCCGCGAAGGGCTGTTGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAggg                                                                                                                                                   <  5:439052/151‑1 (MQ=255)
   ccGCGAAGGGCTGTTGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGga                                                                                                                                                  <  1:318270/151‑1 (MQ=255)
      cgAAGGGCTGTTGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGgag                                                                                                                                                 >  7:901446/1‑149 (MQ=255)
               gTTGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTAt                                                                                                                                        <  5:1858420/149‑1 (MQ=255)
                ttGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGa                                                                                                                                      <  3:823658/150‑1 (MQ=255)
                  gCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACa                                                                                                                                    <  7:188980/150‑1 (MQ=255)
                  gCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAg                                                                                                                                   <  3:1552750/151‑1 (MQ=255)
                     gACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGt                                                                                                                                  <  5:136091/149‑1 (MQ=255)
                     gACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTac                                                                                                                                <  1:1603852/151‑1 (MQ=255)
                        cGACTGCACGCTCGCCGTCTTCAAAAAACTCAACACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTAcacc                                                                                                                             >  1:115829/1‑151 (MQ=255)
                        cGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTAcacc                                                                                                                             >  7:730138/1‑151 (MQ=255)
                             gCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACcacttc                                                                                                                        <  7:550831/151‑1 (MQ=255)
                             gCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACcactt                                                                                                                         <  2:613562/150‑1 (MQ=255)
                                 gCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACcacttcactt                                                                                                                    <  2:865690/151‑1 (MQ=255)
                                  cTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACttt                                                                                                                   >  4:1337817/1‑151 (MQ=255)
                                    cGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACtttt                                                                                                                  <  2:356885/150‑1 (MQ=255)
                                    cGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTc                                                                                                                 <  5:368163/151‑1 (MQ=255)
                                     gCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTc                                                                                                                 >  2:1719860/1‑150 (MQ=255)
                                      ccGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACAACTTCACTTTTCat                                                                                                               >  5:1957031/1‑151 (MQ=255)
                                       cGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCat                                                                                                               >  2:1370456/1‑150 (MQ=255)
                                       cGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCat                                                                                                               >  2:712306/1‑150 (MQ=255)
                                         tCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATc                                                                                                            >  4:1722064/1‑151 (MQ=255)
                                             cAAAAAACTGAAAACCCAAAAACATATGAACTCCGTCAAATGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATcaca                                                                                                         <  5:643158/150‑1 (MQ=255)
                                             cAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATcac                                                                                                          >  6:297523/1‑149 (MQ=255)
                                             cAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAg                                                                                                        <  6:1818557/151‑1 (MQ=255)
                                              aaaaaaCTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATcac                                                                                                          >  6:622647/1‑148 (MQ=255)
                                               aaaaaCTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGAt                                                                                                      <  1:892259/151‑1 (MQ=255)
                                               aaaaaCTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGGTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGAt                                                                                                      <  1:752925/151‑1 (MQ=255)
                                                aaaaCTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGAt                                                                                                      >  3:929858/1‑150 (MQ=255)
                                                aaaaCTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATa                                                                                                     <  6:166041/151‑1 (MQ=255)
                                                 aaaCTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATaa                                                                                                    <  8:1510372/151‑1 (MQ=255)
                                                      cAAACCCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCGCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAg                                                                                                 <  8:1056553/149‑1 (MQ=255)
                                                      cAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGt                                                                                                <  1:1016375/150‑1 (MQ=255)
                                                       aaaaCCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGa                                                                                              <  6:867701/151‑1 (MQ=255)
                                                        aaaCCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGAt                                                                                             >  4:778408/1‑151 (MQ=255)
                                                          aCCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGc                                                                                           >  8:702879/1‑151 (MQ=255)
                                                           ccAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGcc                                                                                          >  6:1149315/1‑151 (MQ=255)
                                                            cAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCa                                                                                         <  4:1887890/151‑1 (MQ=255)
                                                             aaaaaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAg                                                                                        <  1:219217/151‑1 (MQ=255)
                                                             aaaaaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAg                                                                                        <  4:1228261/151‑1 (MQ=255)
                                                              aaaaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGt                                                                                       >  8:1822516/1‑151 (MQ=255)
                                                              aaaaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGt                                                                                       <  4:1409419/151‑1 (MQ=255)
                                                               aaaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTg                                                                                      >  7:472730/1‑151 (MQ=255)
                                                                aaaCTTATCAAGTCCTTCAATGGCGAGCCCTATCGATTCATCAGCACCGAGCTGAATAAAGTTCGCGCTCAGGTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGa                                                                                     <  2:45664/151‑1 (MQ=255)
                                                                aaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTg                                                                                      >  3:319516/1‑150 (MQ=255)
                                                                aaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTg                                                                                      >  3:1852608/1‑150 (MQ=255)
                                                                aaaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGa                                                                                     <  7:1072548/151‑1 (MQ=255)
                                                                 aaCTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGat                                                                                    <  5:946796/151‑1 (MQ=255)
                                                                    ttATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATc                                                                                 >  5:1064831/1‑151 (MQ=255)
                                                                         aaGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGa                                                                            <  2:1526814/151‑1 (MQ=255)
                                                                             ccGTCACTGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGt                                                                         <  7:1474036/150‑1 (MQ=255)
                                                                              cGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCg                                                                       >  7:204341/1‑151 (MQ=255)
                                                                               gTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCCCTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCg                                                                       >  7:1217170/1‑150 (MQ=255)
                                                                                 cAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGaaa                                                                    <  8:204212/151‑1 (MQ=255)
                                                                                 cAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGaa                                                                     >  7:1596082/1‑150 (MQ=255)
                                                                                     ggCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGACACCcg                                                                >  3:1639137/1‑151 (MQ=255)
                                                                                       ccAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCcgcg                                                              >  6:987882/1‑151 (MQ=255)
                                                                                       ccAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCcgcg                                                              <  2:114438/151‑1 (MQ=255)
                                                                                          gCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGGTAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCAt                                                           <  8:518937/151‑1 (MQ=255)
                                                                                          gCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCAt                                                           >  3:1153457/1‑151 (MQ=255)
                                                                                          gCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCAt                                                           <  3:489823/151‑1 (MQ=255)
                                                                                               aTCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTg                                                        <  7:1321070/149‑1 (MQ=255)
                                                                                                 cGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCtt                                                    >  1:1435549/1‑151 (MQ=255)
                                                                                                  gAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTc                                                   >  7:1421183/1‑151 (MQ=255)
                                                                                                    aTCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTc                                                   >  7:347043/1‑149 (MQ=255)
                                                                                                       aaCACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCAGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCa                                               >  1:494010/1‑150 (MQ=255)
                                                                                                         caccacCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCaa                                              >  6:1023713/1‑149 (MQ=255)
                                                                                                            caccGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAaagaag                                         <  4:1484280/151‑1 (MQ=255)
                                                                                                            caccGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAaaga                                           >  4:1613552/1‑149 (MQ=255)
                                                                                                             accGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGc                                        >  3:1050474/1‑151 (MQ=255)
                                                                                                              ccGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGcc                                       <  2:777715/151‑1 (MQ=255)
                                                                                                               cGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCg                                      <  4:396721/151‑1 (MQ=255)
                                                                                                                gAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGa                                     <  6:764642/151‑1 (MQ=255)
                                                                                                                gAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGa                                     <  8:1195453/151‑1 (MQ=255)
                                                                                                                 aGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGAt                                    <  2:1040473/151‑1 (MQ=255)
                                                                                                                   cgtaaTAAAGTTCGGGCTCAGCTCGATAAACGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGAttt                                  <  1:1837768/148‑1 (MQ=255)
                                                                                                                    tGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGAttt                                  >  4:1956371/1‑150 (MQ=255)
                                                                                                                      aaTAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAGCCCGCGCATTTGGCTTCAGCAAAGACGCCGATGTAGt                               >  8:1347597/1‑151 (MQ=255)
                                                                                                                      aaTAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGt                               <  2:1908586/151‑1 (MQ=255)
                                                                                                                       aTAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTg                              >  5:1260337/1‑151 (MQ=255)
                                                                                                                       aTAAAGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTg                              <  1:1715125/151‑1 (MQ=255)
                                                                                                                         aaaGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGc                            >  3:1256072/1‑151 (MQ=255)
                                                                                                                         aaaGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGc                            >  6:1227956/1‑151 (MQ=255)
                                                                                                                          aaGTTCGCGCTCAGCTCGATAATCGCTAAGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGcc                           >  3:623671/1‑151 (MQ=255)
                                                                                                                           aGTTCGCGCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCct                          <  4:946827/151‑1 (MQ=255)
                                                                                                                               cgcgCTCAGCTCGATAATGCGTAGGGAGAGTTTCTATGACAGTACTGCAGTTCAGTTTTCATATCACTGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGGATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTGtt                       <  8:864131/150‑1 (MQ=255)
                                                                                                                               cgcgCTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCttt                      <  6:1587263/151‑1 (MQ=255)
                                                                                                                                gcgcTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTa                     <  5:827830/151‑1 (MQ=255)
                                                                                                                                 cgcTCAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTAc                    >  8:1770559/1‑151 (MQ=255)
                                                                                                                                     cAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTg                 <  5:1486481/150‑1 (MQ=255)
                                                                                                                                     cAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTg                 <  3:137264/150‑1 (MQ=255)
                                                                                                                                     cAGCTCGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTg                 <  4:645164/150‑1 (MQ=255)
                                                                                                                                         gcGATACTCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGAt            <  6:534415/150‑1 (MQ=255)
                                                                                                                                          cGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGAt            <  8:1404646/150‑1 (MQ=255)
                                                                                                                                          cGATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATg           >  3:2563/1‑151 (MQ=255)
                                                                                                                                           gATAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATg           >  8:398000/1‑150 (MQ=255)
                                                                                                                                            aTAATCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATGaa         <  1:1362730/151‑1 (MQ=255)
                                                                                                                                              aaTCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATGAAAg       <  3:1331136/151‑1 (MQ=255)
                                                                                                                                               aTCGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCAAAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCTCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATGAAAgc      <  4:1123922/151‑1 (MQ=255)
                                                                                                                                                 cGCTAGGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGATACCCGCGCATTTGGCTTCAGCAAAGACGCCGATTTAGTGGCCTCTTTACTTGAGGATGAAAgcgg    >  2:1751479/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CAGCCGCGAAGGGCTGTTGCTGACCGACTGCACGCTCGCCGTCTTCAAAAAACTCAAAACCAAAAAACTTATCAAGTCCGTCAATGGCCAGCCCTATCGAATCAACACCACCGAGCTGAATAAAGTTCGCGCTCAGCTCGATAATCGCTAAGGAGAGTTTCTATGACAGTACACCACTTCACTTTTCATATCACAGATAAAAGTGATGCCAGTGATATCCGGGAAGTCGAAACCCGCGCATTTGGCTTCAGCAAAGAAGCCGATTTAGTGGCCTCTTTACTTGAGGATGAAAGCGCGC  >  DH1/3979859‑3980156

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