Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 915,795 A→C Q156P (CAG→CCG)  ptsP → phosphoenolpyruvate‑protein phosphotransferase PtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1915,7950AC97.9% 2131.5 / 9.9 566Q156P (CAG→CCG) ptsPphosphoenolpyruvate‑protein phosphotransferase PtsP
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (288/265);  total (296/270)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 566 total aligned reads displayed.

TTCCGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTGGATCC  >  DH1/915647‑915942
                                                                                                                                                    |                                                                                                                                                   
  ccGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTc                                                                                                                                                  <  5:88856/150‑1 (MQ=255)
    gcgcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGc                                                                                                                                                >  4:782538/1‑150 (MQ=255)
    gcgcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCa                                                                                                                                               <  8:834237/151‑1 (MQ=255)
    gcgcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGAGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGc                                                                                                                                                >  8:483441/1‑150 (MQ=255)
     cgcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTACGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAg                                                                                                                                              <  6:1011543/151‑1 (MQ=255)
       cgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGtt                                                                                                                                            <  7:1588472/151‑1 (MQ=255)
       cgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGtt                                                                                                                                            >  1:644040/1‑151 (MQ=255)
           tttAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGCTGGCAGCTATTCTTTCTCCGTCGCAGTTGACt                                                                                                                                        >  6:1167904/1‑151 (MQ=255)
            ttAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTg                                                                                                                                       >  4:1713387/1‑151 (MQ=255)
                gCGTACCATTTATTCAACGTCGCCAGTTGCTTGGGGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCtt                                                                                                                                   <  6:1559784/151‑1 (MQ=255)
                     ccAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTg                                                                                                                              <  7:1722268/151‑1 (MQ=255)
                      cAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTgg                                                                                                                             <  8:1453529/151‑1 (MQ=255)
                       aaTTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTggg                                                                                                                            <  7:137275/151‑1 (MQ=255)
                       aaTTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTggg                                                                                                                            >  6:598722/1‑151 (MQ=255)
                        attattCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGc                                                                                                                           <  1:1781027/151‑1 (MQ=255)
                             tCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAata                                                                                                                       <  2:607745/150‑1 (MQ=255)
                               aaCGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCg                                                                                                                    >  5:1590867/1‑151 (MQ=255)
                                aCGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGc                                                                                                                   >  6:1747433/1‑151 (MQ=255)
                                 cGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGcc                                                                                                                  >  8:75202/1‑151 (MQ=255)
                                  gTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGcc                                                                                                                  >  1:1124346/1‑150 (MQ=255)
                                  gTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGcc                                                                                                                  >  1:541158/1‑150 (MQ=255)
                                    cGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGcc                                                                                                                  >  3:497024/1‑148 (MQ=255)
                                    cGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAg                                                                                                                <  8:358277/150‑1 (MQ=255)
                                    cGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCGTGTTTGGGCAATATCGcc                                                                                                                  >  1:1234086/1‑148 (MQ=255)
                                    cGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGGGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCan                                                                                                                <  6:1311402/150‑2 (MQ=255)
                                     gCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGACTTGTTTGGGCAATATCGCCAGAc                                                                                                              >  5:1181820/1‑151 (MQ=255)
                                      ccAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCGGAcg                                                                                                             >  3:751060/1‑151 (MQ=255)
                                      ccAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGAcg                                                                                                             >  7:186559/1‑151 (MQ=255)
                                      ccAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGAcg                                                                                                             >  5:1794754/1‑151 (MQ=255)
                                       cAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCTCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGAcgc                                                                                                            <  2:1659181/151‑1 (MQ=255)
                                        aGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGAcgcg                                                                                                           <  1:1395162/151‑1 (MQ=255)
                                         gTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGa                                                                                                          <  2:280027/151‑1 (MQ=255)
                                            gCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATc                                                                                                       <  4:678817/151‑1 (MQ=255)
                                             cTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATcc                                                                                                      >  4:1855316/1‑151 (MQ=255)
                                              ttGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCcg                                                                                                     >  3:1830437/1‑151 (MQ=255)
                                               tGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATcc                                                                                                      <  5:1953260/149‑1 (MQ=255)
                                                   gtACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCa                                                                                                 <  3:444032/150‑1 (MQ=255)
                                                   gtACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCa                                                                                                 <  1:1948557/150‑1 (MQ=255)
                                                    tACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCG‑ATCCGGGCAtt                                                                                               >  3:1060709/1‑150 (MQ=255)
                                                     aCTGGTGGTAAAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTa                                                                                              <  8:1842252/151‑1 (MQ=255)
                                                      cTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTAc                                                                                             >  2:819894/1‑151 (MQ=255)
                                                        ggtggtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCg                                                                                           <  4:384587/151‑1 (MQ=255)
                                                         gtggtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCgg                                                                                          >  7:429352/1‑151 (MQ=255)
                                                           ggtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGgca                                                                                        <  7:865508/151‑1 (MQ=255)
                                                            gtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGgcag                                                                                       <  3:994413/151‑1 (MQ=255)
                                                             tACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGgcagc                                                                                      >  3:1509345/1‑151 (MQ=255)
                                                             tACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGgcagc                                                                                      >  4:537070/1‑151 (MQ=255)
                                                                aGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCAcc                                                                                   >  8:627998/1‑151 (MQ=255)
                                                                 gCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCt                                                                                  <  5:419824/151‑1 (MQ=255)
                                                                      gagaGTTGCGCCAGTATGACGAAAGTGAAGATTCCTTCCTGCTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTGGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGgtgt                                                                             <  1:1854936/151‑1 (MQ=255)
                                                                      gagaGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGgtgt                                                                             >  4:923071/1‑151 (MQ=255)
                                                                       agagTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGtgtg                                                                            <  1:1052373/151‑1 (MQ=255)
                                                                              cgcCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTg                                                                     >  1:803901/1‑151 (MQ=255)
                                                                               gcCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGCCGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGc                                                                    >  6:758945/1‑151 (MQ=255)
                                                                               gcCAGTATGACGAAAGTGAAGAATCCTTCCTGGAGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCAtt                                                                                               >  4:1132683/1‑124 (MQ=255)
                                                                                ccAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGcc                                                                   <  5:58418/151‑1 (MQ=255)
                                                                                ccAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGAAGCACCTGGTGTGGCGATTGcc                                                                   >  8:130899/1‑151 (MQ=255)
                                                                                  aGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGa                                                                 <  4:237770/151‑1 (MQ=255)
                                                                                   gTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGaa                                                                <  5:1447762/151‑1 (MQ=255)
                                                                                    tATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGa                                                                 >  2:514569/1‑149 (MQ=255)
                                                                                          gAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGc                                                         >  4:98824/1‑151 (MQ=255)
                                                                                          gAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGc                                                         >  6:195387/1‑151 (MQ=255)
                                                                                           aaaGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCa                                                        <  5:782059/151‑1 (MQ=255)
                                                                                            aaGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAg                                                       <  8:109206/151‑1 (MQ=255)
                                                                                                gaagaaTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGAt                                                    <  2:1544552/150‑1 (MQ=255)
                                                                                                  agaaTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGcc                                                 >  5:1093539/1‑151 (MQ=255)
                                                                                                     aTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACg                                              >  8:1657987/1‑151 (MQ=255)
                                                                                                      tccttcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGt                                             >  5:1278197/1‑151 (MQ=255)
                                                                                                        cttcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGtt                                            >  2:1463159/1‑150 (MQ=255)
                                                                                                        cttcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGtt                                            >  7:1184166/1‑150 (MQ=255)
                                                                                                        cttcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATAACGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGtt                                            >  3:557728/1‑150 (MQ=255)
                                                                                                         ttcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTAc                                          >  1:1688344/1‑151 (MQ=255)
                                                                                                          tcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTAcc                                         >  3:209037/1‑151 (MQ=255)
                                                                                                          tcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTAc                                          <  1:738203/150‑1 (MQ=255)
                                                                                                           cctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTAcc                                         <  5:1673255/150‑1 (MQ=255)
                                                                                                               gTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTaa                                    <  3:1745561/151‑1 (MQ=255)
                                                                                                                  aCGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATgg                                 <  3:1411082/151‑1 (MQ=255)
                                                                                                                  aCGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATgg                                 <  6:228654/151‑1 (MQ=255)
                                                                                                                   cGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGa                                <  1:1623810/151‑1 (MQ=255)
                                                                                                                   cGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGa                                <  1:51798/151‑1 (MQ=255)
                                                                                                                    gCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGa                                >  4:1210076/1‑150 (MQ=255)
                                                                                                                    gCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGa                                >  5:293525/1‑150 (MQ=255)
                                                                                                                          cACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGgtg                         >  7:1028129/1‑151 (MQ=255)
                                                                                                                          cACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGgtg                         >  6:1433943/1‑151 (MQ=255)
                                                                                                                            cccAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTa                       <  6:409538/151‑1 (MQ=255)
                                                                                                                              cAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATc                     >  7:311305/1‑151 (MQ=255)
                                                                                                                               aGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCa                    <  5:1871535/151‑1 (MQ=255)
                                                                                                                                gATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATc                     >  7:560674/1‑149 (MQ=255)
                                                                                                                                   ggCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCa                <  3:1288868/151‑1 (MQ=255)
                                                                                                                                    gCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCAt               <  3:665421/151‑1 (MQ=255)
                                                                                                                                    gCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCAt               <  2:1108522/151‑1 (MQ=255)
                                                                                                                                        cTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAAc           >  2:1865704/1‑151 (MQ=255)
                                                                                                                                           ttctttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCt        >  7:71435/1‑151 (MQ=255)
                                                                                                                                            tctttctCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGGCGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCt        <  3:810894/150‑1 (MQ=255)
                                                                                                                                             ctttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTgg      >  4:1442172/1‑151 (MQ=255)
                                                                                                                                               ttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTgg      <  3:1655144/149‑1 (MQ=255)
                                                                                                                                               ttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTg       <  2:1330280/148‑1 (MQ=255)
                                                                                                                                               ttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTg       <  1:374844/148‑1 (MQ=255)
                                                                                                                                               ttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTGGAt    >  2:697869/1‑151 (MQ=255)
                                                                                                                                                 ctcCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTGGATcc  >  3:114867/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCCGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTGGATCC  >  DH1/915647‑915942

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: ‑