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,7950AC98.2% 1889.0 / 5.4 497Q156P (CAG→CCG) ptsPphosphoenolpyruvate‑protein phosphotransferase PtsP
Reads supporting (aligned to +/- strand):  ref base A (7/2);  new base C (241/247);  total (248/249)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 497 total aligned reads displayed.

TCCGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTGGATCC  >  DH1/915648‑915942
                                                                                                                                                   |                                                                                                                                                   
tCCGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTc                                                                                                                                                  >  5:1828594/1‑151 (MQ=255)
 ccGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTc                                                                                                                                                  <  2:1736936/150‑1 (MQ=255)
   gcgcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGc                                                                                                                                                >  6:234422/1‑150 (MQ=255)
     gcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGt                                                                                                                                             <  6:681057/151‑1 (MQ=255)
     gcgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGt                                                                                                                                             >  8:284646/1‑151 (MQ=255)
      cgTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTCCGCCAGTATGACGAAAGTGAAGAATCCTGCCTGGTGACGCTTGCCCCCCAGATGGCAGCTATTCATTCTCCGTCGCAGtg                                                                                                                                            >  8:935282/1‑150 (MQ=255)
          tttAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGAc                                                                                                                                         <  3:563045/150‑1 (MQ=255)
            tAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGc                                                                                                                                      >  8:13986/1‑151 (MQ=255)
              ggCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCt                                                                                                                                    >  8:122635/1‑151 (MQ=255)
                cGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTg                                                                                                                                  <  3:128505/151‑1 (MQ=255)
                cGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTg                                                                                                                                  <  2:670994/151‑1 (MQ=255)
                   aCCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGttt                                                                                                                               >  1:1231215/1‑151 (MQ=255)
                   aCCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGtt                                                                                                                                <  1:884101/150‑1 (MQ=255)
                    ccAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTAGAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTg                                                                                                                              <  2:1513157/151‑1 (MQ=255)
                    ccAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTg                                                                                                                              <  7:752893/151‑1 (MQ=255)
                     cAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTgg                                                                                                                             >  1:1508150/1‑151 (MQ=255)
                      aaTTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTggg                                                                                                                            >  1:1159052/1‑151 (MQ=255)
                       attattCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGc                                                                                                                           >  4:1542746/1‑151 (MQ=255)
                       attattCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGc                                                                                                                           >  7:341250/1‑151 (MQ=255)
                            tCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAata                                                                                                                       <  3:879585/150‑1 (MQ=255)
                             cAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATc                                                                                                                     <  7:1091860/151‑1 (MQ=255)
                                 gTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGcc                                                                                                                  >  1:1727135/1‑150 (MQ=255)
                                 gTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCa                                                                                                                 <  4:193682/151‑1 (MQ=255)
                                   cGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGcc                                                                                                                  >  4:1376394/1‑148 (MQ=255)
                                      cAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGAcgc                                                                                                            >  1:596988/1‑151 (MQ=255)
                                       aGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGGCTGCCTTGTTTGGGCAATATCGCCAGAcgcg                                                                                                           >  5:1089186/1‑151 (MQ=255)
                                        gTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGa                                                                                                          <  2:725562/151‑1 (MQ=255)
                                         ttGCTTCGTGTACTGGTGGTACAGCAACGAGAGTTGCCCCAGTATGACGAAACTCAAGAAACCTACCTGCAGACTCTAGCCACCCAGATGCCAGCTATTCTTTATCCGTCGCAGTTCACTACCTGGTATGCGCAATATCGCAAGACGCGa                                                                                                          >  2:935339/1‑150 (MQ=255)
                                              tGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATcc                                                                                                      <  5:587087/149‑1 (MQ=255)
                                                 tgtACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCgcgc                                                                                                  >  8:1081538/1‑150 (MQ=255)
                                                    aCTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCAtt                                                                                               <  7:1391749/150‑1 (MQ=255)
                                                    aCTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTa                                                                                              <  7:168218/151‑1 (MQ=255)
                                                       ggtggtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTAc                                                                                             <  7:448238/149‑1 (MQ=255)
                                                       ggtggtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCg                                                                                           <  5:1760462/151‑1 (MQ=255)
                                                          ggtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGgc                                                                                         >  6:456847/1‑150 (MQ=255)
                                                           gtACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGgcag                                                                                       <  4:96622/151‑1 (MQ=255)
                                                             ccAGCAACGCGCGGTGCGCCAGTATGACGAACGTGACGAATCCTTCCTGGTGACGCTTGCCACCCAGATGTCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGCCGCGAATCCGCGCATTACCGgcagc                                                                                      <  5:1278453/149‑1 (MQ=255)
                                                               aGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCAcc                                                                                   <  4:145667/151‑1 (MQ=255)
                                                                gCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCt                                                                                  <  2:1176249/151‑1 (MQ=255)
                                                                gCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCt                                                                                  <  2:246187/151‑1 (MQ=255)
                                                                 cAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTg                                                                                 <  4:1189666/151‑1 (MQ=255)
                                                                  aaCGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTgg                                                                                <  5:256526/151‑1 (MQ=255)
                                                                    cGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGgtg                                                                              <  2:185714/151‑1 (MQ=255)
                                                                     gagaGTTGCGCCAGTATGACTAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGgtgt                                                                             <  4:27317/151‑1 (MQ=255)
                                                                     gagaGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGgtgt                                                                             >  8:1549894/1‑151 (MQ=255)
                                                                     gagaGTTGCGCCAGTATGACGAAAGTGAAGAATCCTGCCTGGTGACGCTTGCCACCCAGATCGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGgtgt                                                                             <  7:1210068/151‑1 (MQ=255)
                                                                      agagTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGtgtg                                                                            <  4:739277/151‑1 (MQ=255)
                                                                      agagTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATCGCAGCTATTCTTTCTCCCTCGCAGTTGACTGCCTTGTTTGGGCAATATAGCCAGACGCCACTCCGCGCATTACCGGCAGCACCTGGtgtg                                                                            >  5:1406529/1‑151 (MQ=255)
                                                                       gagTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTgg                                                                           >  7:585129/1‑151 (MQ=255)
                                                                       gagTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGT‑ACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGc                                                                          <  7:151956/151‑1 (MQ=255)
                                                                        agTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGc                                                                          >  8:1272163/1‑151 (MQ=255)
                                                                           tGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGAt                                                                       >  8:582039/1‑151 (MQ=255)
                                                                            gcgcCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGAtt                                                                      <  3:1781798/151‑1 (MQ=255)
                                                                            gcgcCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGAtt                                                                      >  6:774303/1‑151 (MQ=255)
                                                                               ccAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCCCCTGGTGTGGCGATTGcc                                                                   >  3:1534784/1‑151 (MQ=255)
                                                                                cAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCg                                                                  >  4:1645117/1‑151 (MQ=255)
                                                                                  gTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGaa                                                                <  8:1490779/151‑1 (MQ=255)
                                                                                    aTGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAgg                                                              <  3:1645787/151‑1 (MQ=255)
                                                                                    aTGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAgg                                                              >  1:1587207/1‑151 (MQ=255)
                                                                                    aTGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAgg                                                              <  1:977680/151‑1 (MQ=255)
                                                                                      gACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGc                                                             <  3:373078/150‑1 (MQ=255)
                                                                                      gACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGc                                                             <  3:236471/150‑1 (MQ=255)
                                                                                      gACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCt                                                            <  1:471510/151‑1 (MQ=255)
                                                                                       aCGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTg                                                           <  2:169758/151‑1 (MQ=255)
                                                                                        cGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCCGCAGCACCTGGTGTGGCGATTGCCGAAGGCTgg                                                          >  2:484156/1‑151 (MQ=255)
                                                                                          aaaGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGc                                                         >  1:1001296/1‑150 (MQ=255)
                                                                                            aGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAgg                                                      <  4:604915/151‑1 (MQ=255)
                                                                                            aGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAgg                                                      <  5:711005/151‑1 (MQ=255)
                                                                                                 agaaTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGCTGcc                                                 >  6:1034391/1‑151 (MQ=255)
                                                                                                   aaTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCAc                                               <  2:427474/151‑1 (MQ=255)
                                                                                                   aaTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCAc                                               <  8:739070/151‑1 (MQ=255)
                                                                                                    aTCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACg                                              <  4:1737817/151‑1 (MQ=255)
                                                                                                       cttcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGtt                                            >  6:1748373/1‑150 (MQ=255)
                                                                                                         tcctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTAc                                          <  8:1345599/150‑1 (MQ=255)
                                                                                                          cctGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCt                                        <  3:719018/151‑1 (MQ=255)
                                                                                                             ggTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCttt                                      >  6:979435/1‑150 (MQ=255)
                                                                                                               tGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAAt                                   >  6:623060/1‑151 (MQ=255)
                                                                                                                gACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATg                                  <  1:30499/151‑1 (MQ=255)
                                                                                                                gACGCGTGCCACCCAGATGGCAGCTATGCTTTCTCCGTCGCAGTTGACTGCCTTGTTGGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATg                                  <  2:1260251/151‑1 (MQ=255)
                                                                                                                      tGCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAAAATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAAc                              >  3:1021300/1‑149 (MQ=255)
                                                                                                                       gCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAgg                           <  8:1821365/151‑1 (MQ=255)
                                                                                                                       gCCACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAgg                           <  3:650560/151‑1 (MQ=255)
                                                                                                                         cACCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGgtg                         <  7:1181684/151‑1 (MQ=255)
                                                                                                                          aCCCAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGgtgt                        <  3:1389500/151‑1 (MQ=255)
                                                                                                                           cccAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTa                       <  7:1475707/151‑1 (MQ=255)
                                                                                                                            ccAGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTAt                      >  7:187497/1‑151 (MQ=255)
                                                                                                                              aGATGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATc                     >  8:103878/1‑150 (MQ=255)
                                                                                                                                aTGGCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAgg                  <  6:1355221/151‑1 (MQ=255)
                                                                                                                                  ggCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGc                 <  1:757566/150‑1 (MQ=255)
                                                                                                                                  ggCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGc                 >  8:1624077/1‑150 (MQ=255)
                                                                                                                                  ggCAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGc                 <  4:455439/150‑1 (MQ=255)
                                                                                                                                    cAGCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATc              <  3:1178659/151‑1 (MQ=255)
                                                                                                                                      gCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCaa            <  7:9627/151‑1 (MQ=255)
                                                                                                                                      gCTATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCaa            <  3:1708242/151‑1 (MQ=255)
                                                                                                                                        tATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACg          >  8:451372/1‑151 (MQ=255)
                                                                                                                                        tATTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGACGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCACCg          >  7:1657764/1‑151 (MQ=255)
                                                                                                                                         aTTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACTGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACg          <  5:480711/150‑1 (MQ=255)
                                                                                                                                         aTTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGc         >  5:1643485/1‑151 (MQ=255)
                                                                                                                                         aTTCTTTCTCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATCCCACGTTACCTTTAATCGTAGAGGAGTATCAGGCATCAACg          <  6:1206257/150‑1 (MQ=255)
                                                                                                                                             tttctCCGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTgg      >  2:1327890/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCCGCGCGTTTTTAGGCGTACCAATTATTCAACGTCGCCAGTTGCTTGGTGTACTGGTGGTACAGCAACGAGAGTTGCGCCAGTATGACGAAAGTGAAGAATCCTTCCTGGTGACGCTTGCCACCCAGATGGCAGCTATTCTTTCTCAGTCGCAGTTGACTGCCTTGTTTGGGCAATATCGCCAGACGCGAATCCGCGCATTACCGGCAGCACCTGGTGTGGCGATTGCCGAAGGCTGGCAGGATGCCACGTTACCTTTAATGGAACAGGTGTATCAGGCATCAACGCTGGATCC  >  DH1/915648‑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: ‑