Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,371,554 C→A L48F (TTG→TTT ydcP ← peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,371,5540CA97.4% 1127.4 / 12.9 307L48F (TTG→TTTydcPpeptidase
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base A (148/150);  total (155/152)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 309 total aligned reads displayed.

ACAATCAGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTCAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGATACGGTCACATTC  >  DH1/2371405‑2371701
                                                                                                                                                     |                                                                                                                                                   
aCAATCAGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTCa                                                                                                                                                    <  4:454914/151‑1 (MQ=255)
 aaaTCAGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTaaa                                                                                                                                                   <  8:1753919/150‑1 (MQ=255)
   aTCAGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGc                                                                                                                                                 <  3:1134235/151‑1 (MQ=255)
     cAGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTa                                                                                                                                               <  4:1894581/151‑1 (MQ=255)
      aGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTAt                                                                                                                                              >  5:1347875/1‑151 (MQ=255)
       gcgcATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTAtt                                                                                                                                             <  2:1521126/151‑1 (MQ=255)
        cgcATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTAtt                                                                                                                                             >  4:844656/1‑150 (MQ=255)
        cgcATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTa                                                                                                                                            <  3:1217600/151‑1 (MQ=255)
         gcATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTAc                                                                                                                                           >  8:1931592/1‑151 (MQ=255)
              gACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGc                                                                                                                                       >  2:1504955/1‑150 (MQ=255)
              gACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGc                                                                                                                                       >  4:656309/1‑150 (MQ=255)
                   cGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATg                                                                                                                                 <  6:121201/151‑1 (MQ=255)
                   cGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATg                                                                                                                                 <  1:1161886/151‑1 (MQ=255)
                    ggTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGTCATTATg                                                                                                                                 >  2:1737063/1‑150 (MQ=255)
                     gtgtGGTAGAGGTCAGTTATGAGCCGTTTCGCGGGTTCCAGCTGATCATCATGCAAAATGGTGGTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACGAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATg                                                                                                                                 <  6:1420841/146‑1 (MQ=255)
                     gTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGAc                                                                                                                               <  7:264371/151‑1 (MQ=255)
                         ggTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGc                                                                                                                           >  3:1762332/1‑151 (MQ=255)
                          gTAGATCTCACTAATCAGGAGTTGCTCCGGTTCCATCTCATCATCAAGCAAAATGTTGTTAATCGTGACTAAAATTTTTGCACCATAACGATTGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGc                                                                                                                           <  4:759270/150‑1 (MQ=255)
                          gTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCa                                                                                                                          <  7:1050444/151‑1 (MQ=255)
                            agagGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCAcc                                                                                                                        >  5:1194247/1‑151 (MQ=255)
                              agGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCa                                                                                                                       >  3:273489/1‑150 (MQ=255)
                                gTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCAGCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCtaa                                                                                                                     >  4:1292912/1‑150 (MQ=255)
                                gTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCaaa                                                                                                                     >  3:704680/1‑150 (MQ=255)
                                   aGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCa                                                                                                                 <  3:1792079/151‑1 (MQ=255)
                                    gTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCATAATGGTGTTAAGTGTGACGAAAATTTTTGCACCATAAGGATGGGCAAACGGCTCCAGGTCGGCAATATCTTTAAAGCTATTAGTGGGATTATGTCGGGCACCAAAACCAg                                                                                                                <  6:1559624/151‑1 (MQ=255)
                                     tAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAgg                                                                                                               >  6:309645/1‑151 (MQ=255)
                                      aaTCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAggg                                                                                                              >  2:1310740/1‑151 (MQ=255)
                                         cAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGgccg                                                                                                           >  7:1859934/1‑151 (MQ=255)
                                         cAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGgcc                                                                                                            <  8:1106574/150‑1 (MQ=255)
                                         cAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCAGTATGACGGGCACCAAAACCAGGgccg                                                                                                           >  2:1645217/1‑151 (MQ=255)
                                          aGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACAAGCTCGGCAATATCTTTAAAGCTATTACTGGAATTATGACGGGCACCAAAACCAGGgccgc                                                                                                          >  7:1703945/1‑151 (MQ=255)
                                           gCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGgccgcc                                                                                                         <  5:987985/151‑1 (MQ=255)
                                           gCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGgccgcc                                                                                                         <  7:1926501/151‑1 (MQ=255)
                                              gTTGCGCGGGTTCCAGCTCATCATAATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGat                                                                                                      >  5:275005/1‑151 (MQ=255)
                                                 gcgcGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGAtata                                                                                                   <  5:928202/151‑1 (MQ=255)
                                                   gcgGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATaaa                                                                                                 <  6:380119/151‑1 (MQ=255)
                                                    cgGGCTCCAGCTCAGCATCATGCAAAATGGTGTTAACCGTGACGAAAATTTTTGCACCATAACGATGGCCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATCACGGCCACCAAAACCAGGGCCGCCGATATAAAc                                                                                                <  2:1155594/151‑1 (MQ=255)
                                                     gggTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACa                                                                                               <  7:1147726/151‑1 (MQ=255)
                                                     gggTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACa                                                                                               <  6:651753/151‑1 (MQ=255)
                                                      ggTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAg                                                                                              <  1:1525684/151‑1 (MQ=255)
                                                      ggTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAg                                                                                              <  3:416752/151‑1 (MQ=255)
                                                         tCCAGCTCATCATCATGCAAATTGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCAt                                                                                           >  1:399969/1‑151 (MQ=255)
                                                          ccAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATc                                                                                          <  1:948367/151‑1 (MQ=255)
                                                           cAGCTCATCATCATGCAAAATGGGGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCg                                                                                         <  4:1695291/151‑1 (MQ=255)
                                                              cTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGc                                                                                       >  6:859618/1‑150 (MQ=255)
                                                                catcatcatGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCAc                                                                                     <  7:606815/150‑1 (MQ=255)
                                                                    atcatGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGGCGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGc                                                                                <  3:873143/151‑1 (MQ=255)
                                                                     tcatGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGc                                                                                >  7:1551758/1‑150 (MQ=255)
                                                                       atGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGTCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCaaa                                                                             <  8:1486485/151‑1 (MQ=255)
                                                                        tGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCaaa                                                                             >  3:993383/1‑150 (MQ=255)
                                                                        tGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCaaa                                                                             >  5:446733/1‑150 (MQ=255)
                                                                          cAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATa                                                                          <  2:300648/151‑1 (MQ=255)
                                                                          cAAAATGGTGCTAAGCGTGACGAAAATTATTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATa                                                                          <  4:1569698/151‑1 (MQ=255)
                                                                            aaaTGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGc                                                                        >  4:1711379/1‑151 (MQ=255)
                                                                            aaaTGGTGTTAAGCGTGACGAAAATTTTTGAACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACGGTGCAAAATAGc                                                                        >  4:1659743/1‑151 (MQ=255)
                                                                             aaTGGTGTTAAGCGGGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCt                                                                       <  8:1271579/151‑1 (MQ=255)
                                                                                ggTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTc                                                                     <  4:606309/150‑1 (MQ=255)
                                                                                 gtgtTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTcgc                                                                   <  6:759800/151‑1 (MQ=255)
                                                                                 gtgtTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTcgc                                                                   >  8:735003/1‑151 (MQ=255)
                                                                                   gtTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCgg                                                                 <  1:1613603/151‑1 (MQ=255)
                                                                                    ttAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCggg                                                                >  1:476570/1‑151 (MQ=255)
                                                                                     tAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTCAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGc                                                               >  7:1254745/1‑151 (MQ=255)
                                                                                     tAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGGACCAAACCCAGGGCCGCCGATATAAACAGCCTCGGCGGCGTGCAAAATAGCGTCGCGGGc                                                               >  6:1957226/1‑151 (MQ=255)
                                                                                      aaGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGCCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGc                                                               >  5:1454760/1‑150 (MQ=255)
                                                                                        gCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAAt                                                            <  1:1094319/151‑1 (MQ=255)
                                                                                         cGTGAAGAAAATTTTTGCACCATATCGATGGGCAAACGGCACCAGCTCGGCAATATATTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATg                                                           <  5:180845/151‑1 (MQ=255)
                                                                                          gTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATgg                                                          <  1:1751094/151‑1 (MQ=255)
                                                                                           tGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATggc                                                         >  7:1400837/1‑151 (MQ=255)
                                                                                              cGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATggcggc                                                      <  5:832639/151‑1 (MQ=255)
                                                                                                  aaTTTTTGCACCATAACGATGGGCAAACGTCACCAGCTCTGCAATATCTTTAAAGCTATTACTGGCATTATGACTGGCACCAAAACCATGTCCGCCGATATAAACAGCATCGACACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATcg                                                  <  8:1432886/151‑1 (MQ=255)
                                                                                                  aaTTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATcg                                                  <  8:1686458/151‑1 (MQ=255)
                                                                                                  aaTTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATcg                                                  >  4:1962510/1‑151 (MQ=255)
                                                                                                     ttttGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATAAAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGt                                               >  4:1338723/1‑151 (MQ=255)
                                                                                                        tGCACCATAACGATGGGCAAACGGCACAAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCCAAACCAGGGCCGCCGATATAAACAGCATAGGCACCGTGCAAAATAGCTTAGCGGGCAATGGCGGCATCGCGTGcc                                            >  1:1251685/1‑151 (MQ=255)
                                                                                                         gCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGcc                                            <  8:30040/150‑1 (MQ=255)
                                                                                                         gCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCg                                           <  5:1629293/151‑1 (MQ=255)
                                                                                                           aCCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCggg                                         <  3:1678786/151‑1 (MQ=255)
                                                                                                            ccATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGc                                        <  6:1821276/151‑1 (MQ=255)
                                                                                                            ccATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGc                                        <  3:1568743/151‑1 (MQ=255)
                                                                                                              aTAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCtt                                      >  5:564431/1‑151 (MQ=255)
                                                                                                              aTAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCtt                                      <  1:183278/151‑1 (MQ=255)
                                                                                                              aTAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCtt                                      <  5:1927578/151‑1 (MQ=255)
                                                                                                               tAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCtt                                      >  6:1664913/1‑150 (MQ=255)
                                                                                                                aaCGATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTa                                     <  2:18158/150‑1 (MQ=255)
                                                                                                                   gATGGGCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAg                                 <  8:1021110/151‑1 (MQ=255)
                                                                                                                       ggCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCCTCGCGTGCCGGGCTTAACAGTTc                              >  5:1537357/1‑150 (MQ=255)
                                                                                                                        gCAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCaa                            <  5:1906293/151‑1 (MQ=255)
                                                                                                                         cAAACGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAg                           <  8:453939/151‑1 (MQ=255)
                                                                                                                           aaCGGCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTc                         >  5:516022/1‑151 (MQ=255)
                                                                                                                               gCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATg                     <  3:1286234/151‑1 (MQ=255)
                                                                                                                               gCACCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACAAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATg                     >  1:1430377/1‑151 (MQ=255)
                                                                                                                                 aCCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATgag                   <  4:1127829/151‑1 (MQ=255)
                                                                                                                                 aCCAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATg                     >  1:1913292/1‑149 (MQ=255)
                                                                                                                                   cAGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATgaga                  >  1:780529/1‑150 (MQ=255)
                                                                                                                                    aGCTCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGagaag                <  4:972178/151‑1 (MQ=255)
                                                                                                                                       tCGGCAATATCTTTCAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGAt              >  1:287429/1‑150 (MQ=255)
                                                                                                                                       tCGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGAt              >  6:1001113/1‑150 (MQ=255)
                                                                                                                                        cGGCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGAt              <  7:1359623/149‑1 (MQ=255)
                                                                                                                                         ggCAATATCTTTAAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGATACg           >  4:368423/1‑151 (MQ=255)
                                                                                                                                                  tttCAAGCAATTACTGGCATTATGACGGGCACCAAAACAAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCCTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGATACGGTCACATTc  >  2:761308/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACAATCAGCGCATCGACACCGGTCTGGTAGAGGTCAGTAATCAGCCGTTGCGCGGGTTCCAGCTCATCATCATGCAAAATGGTGTTAAGCGTGACGAAAATTTTTGCACCATAACGATGGGCAAACGGCACCAGCTCGGCAATATCTTTCAAGCTATTACTGGCATTATGACGGGCACCAAAACCAGGGCCGCCGATATAAACAGCATCGGCACCGTGCAAAATAGCTTCGCGGGCAATGGCGGCATCGCGTGCCGGGCTTAACAGTTCAAGTCGATGAGAAGATACGGTCACATTC  >  DH1/2371405‑2371701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: