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,3710CT98.3% 1938.0 / 10.5 525Q337* (CAG→TAG) phnMalpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Reads supporting (aligned to +/- strand):  ref base C (6/3);  new base T (251/265);  total (257/268)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 528 total aligned reads displayed.

GAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGATTTG  >  DH1/4198224‑4198518
                                                                                                                                                   |                                                                                                                                                   
gAACTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAg                                                                                                                                                   <  3:475246/150‑1 (MQ=255)
 aaCTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGgc                                                                                                                                                 >  7:66870/1‑151 (MQ=255)
  aCTGGCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGgcg                                                                                                                                                <  4:1630978/151‑1 (MQ=255)
     ggCGCAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCtt                                                                                                                                             <  6:813703/151‑1 (MQ=255)
         cAGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAAtc                                                                                                                                         <  5:1129317/151‑1 (MQ=255)
          aGCTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCAGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAAtct                                                                                                                                        >  1:204617/1‑151 (MQ=255)
            cTTGGCCTGCTGGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTcc                                                                                                                                      >  5:713578/1‑151 (MQ=255)
                    gctgGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATcgc                                                                                                                              >  3:843352/1‑151 (MQ=255)
                    gctgGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATcg                                                                                                                               <  3:1271104/150‑1 (MQ=255)
                    gctgGATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTATACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATc                                                                                                                                <  5:1396848/149‑1 (MQ=255)
                        gATATCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCgggg                                                                                                                          <  3:1744188/151‑1 (MQ=255)
                         atatCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGt                                                                                                                         <  8:1457446/151‑1 (MQ=255)
                          tatCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTg                                                                                                                        >  2:1095849/1‑151 (MQ=255)
                           atCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTg                                                                                                                        >  7:389426/1‑150 (MQ=255)
                            tCCTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGAt                                                                                                                      >  6:1070915/1‑151 (MQ=255)
                             ccTCTCTTCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGAt                                                                                                                      <  1:971798/150‑1 (MQ=255)
                                  attCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTGGGGTCTGCGATGACCAGAGCAACCGCTTTACTCTGACGCAGGCGGTGAAGCTTGTGACTAAAAATCCAGCGTAGGCGCTTAATCTAAAGGATCGCGGGGTGATTGGc                                                                                                                  <  4:464702/148‑1 (MQ=255)
                                   ttCCGACTACTACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCg                                                                                                                 >  4:1747766/1‑149 (MQ=255)
                                        actactACCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCaa                                                                                                          <  1:1183926/151‑1 (MQ=255)
                                         ctactaCCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCaa                                                                                                          >  5:1535831/1‑150 (MQ=255)
                                           actacCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAAcg                                                                                                       <  6:88347/151‑1 (MQ=255)
                                           actacCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAAc                                                                                                        <  8:1658972/150‑1 (MQ=255)
                                           actacCCCGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAAc                                                                                                        <  2:589348/150‑1 (MQ=255)
                                                cccGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCg                                                                                                  >  1:1648017/1‑151 (MQ=255)
                                                 ccGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGa                                                                                                 <  5:1815917/151‑1 (MQ=255)
                                                  cGCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCCTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGAc                                                                                                <  1:710590/151‑1 (MQ=255)
                                                   gCCAGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGAcc                                                                                               >  3:372507/1‑151 (MQ=255)
                                                      aGCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTgg                                                                                            >  1:292330/1‑151 (MQ=255)
                                                      aGCCTGCGCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTgg                                                                                            <  4:848602/151‑1 (MQ=255)
                                                       gCCTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGt                                                                                           <  6:1337841/151‑1 (MQ=255)
                                                         cTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGc                                                                                         >  8:1031561/1‑151 (MQ=255)
                                                         cTGCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGc                                                                                         >  2:1245879/1‑151 (MQ=255)
                                                           gCTCGTTGCGGCATGTCGCGTCGCCGATGACGAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTg                                                                                       <  2:1617849/151‑1 (MQ=255)
                                                           gCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGc                                                                                         <  4:115733/149‑1 (MQ=255)
                                                           gCTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTg                                                                                       <  7:171343/151‑1 (MQ=255)
                                                            cTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTgg                                                                                      >  4:697923/1‑151 (MQ=255)
                                                            cTCGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTgg                                                                                      >  3:63515/1‑151 (MQ=255)
                                                              cGATGCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGgc                                                                                     <  5:881599/150‑1 (MQ=255)
                                                                  gCGGCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATc                                                                                <  5:174225/151‑1 (MQ=255)
                                                                    ggCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGc                                                                              >  1:1020928/1‑151 (MQ=255)
                                                                     gCATTTCGCGTCGCCGATGACCAGAGCACCCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCa                                                                             <  7:542127/151‑1 (MQ=255)
                                                                     gCATTTCGCGTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGc                                                                              >  4:1483829/1‑150 (MQ=255)
                                                                           cgcgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGAc                                                                        >  7:484772/1‑150 (MQ=255)
                                                                           cgcgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGAc                                                                        >  1:578081/1‑150 (MQ=255)
                                                                            gcgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACaa                                                                      <  8:1324569/151‑1 (MQ=255)
                                                                             cgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAt                                                                     <  3:1810963/151‑1 (MQ=255)
                                                                             cgTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAt                                                                     >  5:832793/1‑151 (MQ=255)
                                                                              gTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtc                                                                    >  2:1767083/1‑151 (MQ=255)
                                                                              gTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtc                                                                    <  8:288444/151‑1 (MQ=255)
                                                                              gTCGCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtc                                                                    >  2:455328/1‑151 (MQ=255)
                                                                                 gCCGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcata                                                                 <  7:1305551/151‑1 (MQ=255)
                                                                                  ccGATGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatat                                                                >  8:1556336/1‑151 (MQ=255)
                                                                                     aTGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatattca                                                             <  7:949046/151‑1 (MQ=255)
                                                                                     aTGACCAGAGCAACAGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGAGTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATTGCAAGGACAAtcatattc                                                              >  3:1117344/1‑150 (MQ=255)
                                                                                      tGACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatattcat                                                            >  4:372257/1‑151 (MQ=255)
                                                                                       gACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcttattcata                                                           <  1:488058/151‑1 (MQ=255)
                                                                                       gACCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatattcat                                                            >  7:838116/1‑150 (MQ=255)
                                                                                        aCCAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAAtcatattcatat                                                          <  5:1624539/151‑1 (MQ=255)
                                                                                          cAGAGCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCg                                                        <  5:1285134/151‑1 (MQ=255)
                                                                                             agCAACCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGAcc                                                     <  8:1181664/151‑1 (MQ=255)
                                                                                                aaCCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACg                                                  <  6:519671/151‑1 (MQ=255)
                                                                                                 aCCGCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGt                                                 >  2:978761/1‑151 (MQ=255)
                                                                                                    gCTTTACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTg                                              >  6:44683/1‑151 (MQ=255)
                                                                                                       ttACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGc                                            <  1:1103248/150‑1 (MQ=255)
                                                                                                        ttcGCTGCCGCAGGGGGTGTAGCTGGTGACTAAAAATCGAGCGTTGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTGATATCGACCACGTCTGGCg                                           <  8:880315/148‑1 (MQ=255)
                                                                                                        tACGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGt                                          >  1:1706152/1‑151 (MQ=255)
                                                                                                         aCGCTGCCGCAGGCGGTGAAGGTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGGGGGTTGATTGGCGAGGGCAACCGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTc                                         <  7:1549901/151‑1 (MQ=255)
                                                                                                         aCGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATCTTCATATCGACCACGTCTGGCGTc                                         >  1:792107/1‑151 (MQ=255)
                                                                                                         aCGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTc                                         >  5:978601/1‑151 (MQ=255)
                                                                                                          cGCTGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTATGCGCTTAATCTCCAGGATCGCGGGGTGATTTGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCa                                        <  4:1300302/151‑1 (MQ=255)
                                                                                                             tGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAggg                                     >  1:1326829/1‑151 (MQ=255)
                                                                                                             tGCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAggg                                     >  3:275669/1‑151 (MQ=255)
                                                                                                              gcggCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAggg                                     <  6:254334/147‑1 (MQ=255)
                                                                                                              gCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAggg                                     <  2:141165/150‑1 (MQ=255)
                                                                                                              gCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGt                                    >  5:617479/1‑151 (MQ=255)
                                                                                                              gCCGCAGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGt                                    >  2:1353245/1‑151 (MQ=255)
                                                                                                                   aGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTaaa                                 >  7:1390559/1‑149 (MQ=255)
                                                                                                                   aGGCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTaaa                                 >  2:401318/1‑149 (MQ=255)
                                                                                                                     gCGGTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGCCCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAggg                             >  6:1748496/1‑151 (MQ=255)
                                                                                                                      gggTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGgt                            <  1:853767/150‑1 (MQ=255)
                                                                                                                      cGGTGAAGCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGgt                            >  5:388235/1‑151 (MQ=255)
                                                                                                                        gTGAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGgtgt                          <  8:443065/151‑1 (MQ=255)
                                                                                                                          gAAGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAggg                                     <  6:1509440/138‑1 (MQ=255)
                                                                                                                           aaGGTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGACGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCt                       <  6:931521/151‑1 (MQ=255)
                                                                                                                            aGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtg                      >  3:1438579/1‑151 (MQ=255)
                                                                                                                            aGCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtg                      <  3:1238056/151‑1 (MQ=255)
                                                                                                                             gCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtga                     <  2:206003/151‑1 (MQ=255)
                                                                                                                             gCTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtg                      >  6:387274/1‑150 (MQ=255)
                                                                                                                             gCTGGTGACTAAAAATCCAGCGCAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtg                      >  5:1031247/1‑150 (MQ=255)
                                                                                                                              cTGGTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtgat                    >  7:1782465/1‑151 (MQ=255)
                                                                                                                                ggTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCtgatg                   <  2:764789/150‑1 (MQ=255)
                                                                                                                                 gTGACTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTgatgat                 <  7:1177867/151‑1 (MQ=255)
                                                                                                                                    aCTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATggg              >  8:67798/1‑151 (MQ=255)
                                                                                                                                     cTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATggg              >  5:1190842/1‑150 (MQ=255)
                                                                                                                                     cTAAAAATCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATggg              >  5:487724/1‑150 (MQ=255)
                                                                                                                                       aaaaaTCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGaaa           <  8:1772088/151‑1 (MQ=255)
                                                                                                                                          aaTCCAGCGTAGGCTCTTATTCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGTTAAAAGGGTGTTCTGATGATGGGAAAACt        <  6:609444/151‑1 (MQ=255)
                                                                                                                                          aaTCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACt        <  4:1089512/151‑1 (MQ=255)
                                                                                                                                           aTCCAGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTg       >  7:257337/1‑151 (MQ=255)
                                                                                                                                               aGCGTAGGCGCTTAATCTCCAGGATCGCGGGGTGATTGGCGAGGGCAAACGCGCCGACCTGGTGCTGGCGCATCGCAAGGACAATCATATTCATATCGACCACGTCTGGCGTCAGGGTAAAAGGGTGTTCTGATGATGGGAAAACTGAttt   >  8:74313/1‑151 (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: