Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,372,968 G→A D191N (GAT→AAT)  ydcO → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,372,9680GA97.4% 1403.4 / 13.7 393D191N (GAT→AAT) ydcOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base A (172/208);  total (177/216)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 393 total aligned reads displayed.

TGCCAGTCTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTGATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGATATT  >  DH1/2372820‑2373115
                                                                                                                                                    |                                                                                                                                                   
  ccAGTCTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAAtgt                                                                                                                                                 <  1:511229/151‑1 (MQ=255)
   cAGTCTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAAtgtt                                                                                                                                                <  1:1249960/151‑1 (MQ=255)
   cAGTCTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAAtgtt                                                                                                                                                <  7:1303441/151‑1 (MQ=255)
     gTCTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAAtgttgt                                                                                                                                              <  7:157348/151‑1 (MQ=255)
       cTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCt                                                                                                                                            >  3:554586/1‑151 (MQ=255)
          gACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTa                                                                                                                                         <  5:1002533/151‑1 (MQ=255)
           aCGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTa                                                                                                                                         >  6:708420/1‑150 (MQ=255)
             ggTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAAc                                                                                                                                      >  5:1574390/1‑151 (MQ=255)
              gTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAAcc                                                                                                                                     <  6:1273819/151‑1 (MQ=255)
               tCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAAccc                                                                                                                                    >  5:1089350/1‑151 (MQ=255)
                cAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCg                                                                                                                                   <  2:352154/151‑1 (MQ=255)
                cAATTTACGTTGGGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAAcccc                                                                                                                                   <  3:1734797/151‑2 (MQ=255)
                 aaTTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATCATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTATAAACCCGt                                                                                                                                  >  6:1754290/1‑151 (MQ=255)
                  aTTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGtt                                                                                                                                 >  3:1379173/1‑151 (MQ=255)
                       cGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTccc                                                                                                                            <  6:1116580/151‑1 (MQ=255)
                          tgtgtgGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCa                                                                                                                          <  5:1361350/150‑1 (MQ=255)
                          tgtgtgGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCAc                                                                                                                         >  7:1169530/1‑151 (MQ=255)
                           gtgtgGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCAc                                                                                                                         <  5:561711/150‑1 (MQ=255)
                             gtgGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACtt                                                                                                                       <  6:1333538/150‑1 (MQ=255)
                             gtgGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTta                                                                                                                      <  7:886960/151‑1 (MQ=255)
                               ggAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTtat                                                                                                                     <  7:1262296/150‑1 (MQ=255)
                               ggAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTtat                                                                                                                     <  4:1763369/150‑1 (MQ=255)
                               ggAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTtat                                                                                                                     >  4:1352197/1‑150 (MQ=255)
                               ggAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTtat                                                                                                                     >  3:939766/1‑150 (MQ=255)
                                 aaGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATAtt                                                                                                                  <  5:344225/151‑1 (MQ=255)
                                 aaGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATAtt                                                                                                                  >  1:1685826/1‑151 (MQ=255)
                                 aaGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATAtt                                                                                                                  <  8:560363/151‑1 (MQ=255)
                                   gTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTAc                                                                                                                <  8:894856/151‑1 (MQ=255)
                                       gTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTAcccc                                                                                                             <  8:632440/150‑1 (MQ=255)
                                         tGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTg                                                                                                           <  7:651420/150‑1 (MQ=255)
                                             ggTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGAtttt                                                                                                      >  3:69023/1‑151 (MQ=255)
                                             ggTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGAttt                                                                                                       <  5:868022/150‑1 (MQ=255)
                                                aTGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTc                                                                                                    <  5:496044/150‑1 (MQ=255)
                                                      ggcAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGcg                                                                                             >  6:372832/1‑150 (MQ=255)
                                                         aaCCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTcac                                                                                          <  2:1642232/151‑1 (MQ=255)
                                                          aCCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTcaca                                                                                         <  3:1451536/151‑1 (MQ=255)
                                                            cAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGc                                                                                       <  8:433186/151‑1 (MQ=255)
                                                                 ccGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTg                                                                                    >  6:1673901/1‑149 (MQ=255)
                                                                  cGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGc                                                                                 >  2:196517/1‑151 (MQ=255)
                                                                  cGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGc                                                                                 >  7:345169/1‑151 (MQ=255)
                                                                   gTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCg                                                                                <  1:228457/151‑1 (MQ=255)
                                                                     tGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGtt                                                                              >  6:64538/1‑151 (MQ=255)
                                                                      gcgcCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTg                                                                             <  3:1678140/151‑1 (MQ=255)
                                                                      gcgcCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTg                                                                             <  4:91798/151‑1 (MQ=255)
                                                                        gccgcGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGc                                                                            >  6:439924/1‑150 (MQ=255)
                                                                            cgcTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTcc                                                                       >  3:75525/1‑151 (MQ=255)
                                                                             gcTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTccc                                                                      >  1:94906/1‑151 (MQ=255)
                                                                              cTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTcccc                                                                     >  4:420908/1‑151 (MQ=255)
                                                                               tATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTAGCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTccccc                                                                    >  5:40891/1‑151 (MQ=255)
                                                                                aTGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTccccc                                                                    <  8:334461/150‑1 (MQ=255)
                                                                                    ggTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCttttt                                                               <  4:29027/151‑1 (MQ=255)
                                                                                       aaTTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCt                                                             <  7:714154/150‑1 (MQ=255)
                                                                                          tGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGt                                                          <  8:1662885/150‑1 (MQ=255)
                                                                                           gCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTg                                                         <  8:211442/150‑1 (MQ=255)
                                                                                            ccGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGCTGCACTCCCCCTTTTTCTGGTGAc                                                       <  2:121939/151‑1 (MQ=255)
                                                                                             cgcgATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTCATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACg                                                      >  1:1400885/1‑151 (MQ=255)
                                                                                               cgATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGAt                                                    >  7:201118/1‑151 (MQ=255)
                                                                                                gatgatTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCCTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATg                                                   >  6:1547302/1‑151 (MQ=255)
                                                                                                  tgatTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGc                                                 >  2:80021/1‑151 (MQ=255)
                                                                                                  tgatTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGc                                                 >  5:167716/1‑151 (MQ=255)
                                                                                                  tgatTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGc                                                 >  6:881683/1‑151 (MQ=255)
                                                                                                   gatTATTGGGATCGTGATCGTCATCGCTCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCa                                                <  1:1016149/151‑1 (MQ=255)
                                                                                                   gatTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCa                                                <  7:771019/151‑1 (MQ=255)
                                                                                                     ttattGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATc                                              >  5:117916/1‑151 (MQ=255)
                                                                                                         tGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCa                                           <  4:1499362/150‑1 (MQ=255)
                                                                                                          gggATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCaaa                                         <  7:417829/151‑1 (MQ=255)
                                                                                                          gggATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCaa                                          >  5:1815709/1‑150 (MQ=255)
                                                                                                            gatcgtgatcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAAc                                       >  4:162650/1‑151 (MQ=255)
                                                                                                              tcgtgatcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTTCCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACg                                      <  2:450463/150‑1 (MQ=255)
                                                                                                               cgtgatcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGc                                     <  7:573244/150‑1 (MQ=255)
                                                                                                                gtgatcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCAc                                   >  6:1450934/1‑151 (MQ=255)
                                                                                                                 tgatcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCAcc                                  >  3:495522/1‑151 (MQ=255)
                                                                                                                 tgatcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCAcc                                  >  2:896879/1‑151 (MQ=255)
                                                                                                                   atcgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCgg                                <  8:1291670/151‑1 (MQ=255)
                                                                                                                     cgtCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGt                              <  4:661513/151‑1 (MQ=255)
                                                                                                                       tCATCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCCCCGGGTAt                            >  1:374553/1‑151 (MQ=255)
                                                                                                                         aTCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCg                          <  8:736986/151‑1 (MQ=255)
                                                                                                                         aTCGCGCAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCg                          <  7:989664/151‑1 (MQ=255)
                                                                                                                            gcgcAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCgcag                       <  8:1527847/151‑1 (MQ=255)
                                                                                                                            gcgcAAGGTGACGATGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCgcac                       >  6:281995/1‑150 (MQ=255)
                                                                                                                              gcAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCgcagca                     <  3:358488/151‑1 (MQ=255)
                                                                                                                              gcAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCgcagc                      >  4:381527/1‑150 (MQ=255)
                                                                                                                               cAAGGTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCgcagca                     >  6:1002193/1‑150 (MQ=255)
                                                                                                                                aaGGTGACGTTGTCACAACTAATGTTGTCTTTACACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAAt                   <  8:1763663/151‑1 (MQ=255)
                                                                                                                                  ggTGACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGa                 <  6:1848735/151‑1 (MQ=255)
                                                                                                                                     gACGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGaaa               <  3:1617960/150‑1 (MQ=255)
                                                                                                                                       cGTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAagca            <  4:870115/151‑1 (MQ=255)
                                                                                                                                        gTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAagcag           <  1:1500787/151‑1 (MQ=255)
                                                                                                                                        gTTGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAagcag           <  5:708735/151‑1 (MQ=255)
                                                                                                                                          tGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCt         >  7:1746203/1‑151 (MQ=255)
                                                                                                                                          tGTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCt         >  3:207134/1‑151 (MQ=255)
                                                                                                                                           gTCACAACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGTAGCt         <  3:1021610/150‑1 (MQ=255)
                                                                                                                                             cacaACTGATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTgg       <  6:773557/150‑1 (MQ=255)
                                                                                                                                              acaACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGat     >  1:640023/1‑151 (MQ=255)
                                                                                                                                              acaACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGat     <  2:646370/151‑1 (MQ=255)
                                                                                                                                               caACTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGat     >  3:1193412/1‑150 (MQ=255)
                                                                                                                                               aaaCTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGata    <  2:23906/150‑1 (MQ=255)
                                                                                                                                                aaCTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGatat   <  4:1162237/151‑1 (MQ=255)
                                                                                                                                                aaCTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGatat   <  1:757270/151‑1 (MQ=255)
                                                                                                                                                aaCTAATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGatat   <  1:1123310/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCAGTCTGGACGGTCAATTTACGTTGTGTGGAAGTATGTTGCTGGTATGGCTGGCAACCAAGGCCGTTGCGCCGCGCTATGCGGTAATTGCCGCGATGATTATTGGGATCGTGATCGTCATCGCGCAAGGTGACGTTGTCACAACTGATGTTGTCTTTAAACCCGTTCTCCCCACTTATATTACCCCTGATTTTTCGTTTGCTCACAGCCTGAGCGTTGCACTCCCCCTTTTTCTGGTGACGATGGCATCGCAAAACGCACCGGGTATCGCAGCAATGAAAGCAGCTGGATATT  >  DH1/2372820‑2373115

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