Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,198,371 C→T Q337* (CAG→TAG)  phnM → alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,198,3710CT99.2% 2397.2 / 1.9 638Q337* (CAG→TAG) phnMalpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (311/321);  total (313/325)
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.51e-01

Only 100 of 645 total aligned reads displayed.

GAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGATTTG  >  DH1/4198224‑4198518
                                                                                                                                                   |                                                                                                                                                   
gAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCTGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAg                                                                                                                                                   <  8:304116/150‑1 (MQ=255)
gAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAg                                                                                                                                                   <  6:1229748/150‑1 (MQ=255)
gAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAg                                                                                                                                                   <  7:1349435/150‑1 (MQ=255)
 aaCTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGgc                                                                                                                                                 >  3:1009389/1‑151 (MQ=255)
     ggCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGgcgc                                                                                                                                               <  7:1488726/149‑1 (MQ=255)
     ggCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCtt                                                                                                                                             >  3:1709787/1‑151 (MQ=255)
      gcgcAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTa                                                                                                                                            <  2:290884/151‑1 (MQ=255)
       cgcAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTaa                                                                                                                                           <  1:1225172/151‑1 (MQ=255)
          aGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAAtct                                                                                                                                        >  7:801033/1‑151 (MQ=255)
          aGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAAtct                                                                                                                                        >  1:167233/1‑151 (MQ=255)
             ttGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTcc                                                                                                                                      >  4:1732065/1‑150 (MQ=255)
             ttGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTcc                                                                                                                                      >  2:1652562/1‑150 (MQ=255)
                gCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAgg                                                                                                                                   >  7:688611/1‑150 (MQ=255)
                gCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAgg                                                                                                                                   >  5:149185/1‑150 (MQ=255)
                 ccTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGAt                                                                                                                                 <  1:1013010/151‑1 (MQ=255)
                   tgctgGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATcg                                                                                                                               >  1:505712/1‑151 (MQ=255)
                    gctgGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTTCTCGATTCGGATTTTCGATTAGCAGATGAAAAGAGCAACCTATTTACGCTTCCGCAGGAGGTGAAGCTTGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATc                                                                                                                                <  6:1853951/149‑1 (MQ=255)
                    gctgGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATc                                                                                                                                <  5:1960816/149‑1 (MQ=255)
                       ggATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATcgcg                                                                                                                             <  8:1506323/149‑1 (MQ=255)
                        gATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCgggg                                                                                                                          <  5:1006017/151‑1 (MQ=255)
                           atCCTCTCTTCCGACTACTACCCCGCCAGCCTGCCCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTg                                                                                                                        >  7:640934/1‑150 (MQ=255)
                             ccTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGAt                                                                                                                      <  1:1333777/150‑1 (MQ=255)
                              cactctTCCGACTCCTACCCCGCCACCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTCCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTg                                                                                                                    <  8:1322830/149‑1 (MQ=255)
                                ctctTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGc                                                                                                                  >  4:1043083/1‑151 (MQ=255)
                                  ctTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGGCGCCGATGACCAGCGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTgg                                                                                                                   <  8:1928252/148‑1 (MQ=255)
                                    tCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCg                                                                                                                 <  7:541881/148‑1 (MQ=255)
                                     ccGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAgg                                                                                                              <  2:716123/150‑1 (MQ=255)
                                      cGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGc                                                                                                            >  7:453832/1‑151 (MQ=255)
                                           actacCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCaaa                                                                                                         <  7:1552014/149‑1 (MQ=255)
                                           actacCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCaaa                                                                                                         <  1:608202/149‑1 (MQ=255)
                                           actacCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAAcg                                                                                                       >  3:931977/1‑151 (MQ=255)
                                              acCCCGCCAGCCTGCTCGATGCGGGATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATGCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAAcgc                                                                                                      <  2:1018273/149‑1 (MQ=255)
                                               ccccGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGcc                                                                                                   >  4:1877680/1‑151 (MQ=255)
                                                    ccAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCt                                                                                              >  1:1626640/1‑151 (MQ=255)
                                                     cAGCCTGCTCGAGGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTg                                                                                             <  3:1235509/151‑1 (MQ=255)
                                                        ccTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTg                                                                                          >  8:116593/1‑151 (MQ=255)
                                                          tGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCt                                                                                        >  1:1774889/1‑151 (MQ=255)
                                                           gCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCt                                                                                        <  2:149754/150‑1 (MQ=255)
                                                           gCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTg                                                                                       >  5:1704907/1‑151 (MQ=255)
                                                           gCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTg                                                                                       >  7:1830226/1‑151 (MQ=255)
                                                           gCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTg                                                                                       <  7:1739455/151‑1 (MQ=255)
                                                              cGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGgcg                                                                                    >  8:169536/1‑151 (MQ=255)
                                                              cGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGgcg                                                                                    >  7:1136832/1‑151 (MQ=255)
                                                              cGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGgc                                                                                     <  2:193852/150‑1 (MQ=255)
                                                                  gCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATc                                                                                <  7:289541/151‑1 (MQ=255)
                                                                   cGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCg                                                                               <  1:904405/151‑1 (MQ=255)
                                                                   cGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCg                                                                               <  7:953794/151‑1 (MQ=255)
                                                                       aTTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCa                                                                             >  3:738377/1‑149 (MQ=255)
                                                                           cgcgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGAc                                                                        >  2:1839873/1‑150 (MQ=255)
                                                                            gcgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACaa                                                                      <  3:658764/151‑1 (MQ=255)
                                                                             cgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAt                                                                     >  3:462121/1‑151 (MQ=255)
                                                                              gTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtc                                                                    <  3:303368/151‑1 (MQ=255)
                                                                                cGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcat                                                                  <  3:1103/151‑1 (MQ=255)
                                                                                 gCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGAGTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCACGGACAAtcat                                                                  >  8:846177/1‑150 (MQ=255)
                                                                                 gCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcat                                                                  >  7:225357/1‑150 (MQ=255)
                                                                                  ccGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatat                                                                <  5:1177879/151‑1 (MQ=255)
                                                                                  ccGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatat                                                                >  6:1073136/1‑151 (MQ=255)
                                                                                   cGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatatt                                                               >  4:733735/1‑151 (MQ=255)
                                                                                    gATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatattc                                                              >  3:89581/1‑151 (MQ=255)
                                                                                            gagCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGAc                                                      <  3:1405508/151‑1 (MQ=255)
                                                                                             agCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGAcc                                                     <  6:1611874/151‑1 (MQ=255)
                                                                                              gCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCa                                                    <  8:560362/151‑1 (MQ=255)
                                                                                               cAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCAc                                                   <  8:688575/151‑1 (MQ=255)
                                                                                                aaCCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACg                                                  <  5:136697/151‑1 (MQ=255)
                                                                                                aaCCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACg                                                  <  4:597132/151‑1 (MQ=255)
                                                                                                   cGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCt                                               <  5:974020/151‑1 (MQ=255)
                                                                                                    gCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTg                                              >  6:920155/1‑151 (MQ=255)
                                                                                                        tACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCg                                           <  4:420671/150‑1 (MQ=255)
                                                                                                         aCGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGt                                          <  6:644376/150‑1 (MQ=255)
                                                                                                              gCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAgg                                      <  1:1558946/149‑1 (MQ=255)
                                                                                                              gCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGt                                    <  3:1932441/151‑1 (MQ=255)
                                                                                                                   aGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTaaa                                 >  1:1464250/1‑149 (MQ=255)
                                                                                                                        gTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGgtgt                          >  1:19781/1‑151 (MQ=255)
                                                                                                                          gAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTc                        <  7:1242890/151‑1 (MQ=255)
                                                                                                                           aaGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGAGCACGTCTGGCGTCAGGGTAAAAGGGTGTTCt                       >  3:162375/1‑151 (MQ=255)
                                                                                                                           aaGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCt                       <  5:805849/151‑1 (MQ=255)
                                                                                                                           aaGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCt                       >  2:1251066/1‑151 (MQ=255)
                                                                                                                           aaGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCt                       <  6:426933/151‑1 (MQ=255)
                                                                                                                            aGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtg                      >  5:317546/1‑151 (MQ=255)
                                                                                                                            aGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtg                      <  5:584479/151‑1 (MQ=255)
                                                                                                                               tGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtgatg                   >  2:493088/1‑151 (MQ=255)
                                                                                                                                ggTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtgat                    <  8:727367/149‑1 (MQ=255)
                                                                                                                                 gTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTgatgat                 <  4:1513505/151‑1 (MQ=255)
                                                                                                                                  tGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGatgatg                >  2:856446/1‑151 (MQ=255)
                                                                                                                                   gACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATgg               >  3:1462702/1‑151 (MQ=255)
                                                                                                                                    aCTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATggg              <  4:823105/151‑1 (MQ=255)
                                                                                                                                    aCTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATggg              >  5:704315/1‑151 (MQ=255)
                                                                                                                                        aaaaTCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGaaa           >  6:349440/1‑150 (MQ=255)
                                                                                                                                         aaaTCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGaaac          >  4:1201268/1‑149 (MQ=255)
                                                                                                                                           aTCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGAGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTg       >  7:1910479/1‑151 (MQ=255)
                                                                                                                                             ccAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAACCTGAt     <  6:1795625/151‑1 (MQ=255)
                                                                                                                                             ccAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGa      <  5:497505/150‑1 (MQ=255)
                                                                                                                                             ccAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGAt     <  4:1377671/151‑1 (MQ=255)
                                                                                                                                             ccAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGAt     <  8:379774/151‑1 (MQ=255)
                                                                                                                                             ccAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGAt     <  5:1836468/151‑1 (MQ=255)
                                                                                                                                              cAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGAtt    >  6:145663/1‑151 (MQ=255)
                                                                                                                                              cAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGAtt    >  7:141466/1‑151 (MQ=255)
                                                                                                                                                gcgtAGGCGCTTAATCTCCATGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGATTTg  <  5:1066217/151‑1 (MQ=255)
                                                                                                                                                gcgtAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGATTTg  >  7:100647/1‑151 (MQ=255)
                                                                                                                                                gcgtAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGATTTg  <  4:200708/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGATTTG  >  DH1/4198224‑4198518

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