Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,573,476 C→A A146A (GCG→GCT oppF ← oligopeptide transport ATP‑binding protein OppF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,573,4760CA98.6% 1354.9 / 6.1 362A146A (GCG→GCToppFoligopeptide transport ATP‑binding protein OppF
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (164/192);  total (167/195)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 373 total aligned reads displayed.

CCAGCGCCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATCGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGCGAACGGCA  >  DH1/2573328‑2573623
                                                                                                                                                    |                                                                                                                                                   
ccAGCGCCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGc                                                                                                                                                   <  6:659765/151‑1 (MQ=255)
  aGCGCCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCt                                                                                                                                                 <  8:1808138/151‑1 (MQ=255)
    cgcCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTc                                                                                                                                               <  4:263047/151‑1 (MQ=255)
    cgcCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCACGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCtt                                                                                                                                                >  7:385697/1‑150 (MQ=255)
     gcCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTc                                                                                                                                               >  4:429053/1‑150 (MQ=255)
     gcCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTc                                                                                                                                               >  8:87220/1‑150 (MQ=255)
      ccGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCAc                                                                                                                                             <  6:774114/151‑1 (MQ=255)
       cgccacCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATCCGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCAcg                                                                                                                                            <  8:1736733/148‑1 (MQ=255)
       cGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGGATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACAAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCAcg                                                                                                                                            >  2:1350135/1‑151 (MQ=255)
       cGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCAcg                                                                                                                                            >  4:1569209/1‑151 (MQ=255)
       cGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCAcg                                                                                                                                            >  4:1154698/1‑151 (MQ=255)
         acacCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCAcgcg                                                                                                                                          <  7:934616/151‑1 (MQ=255)
           acCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCt                                                                                                                                        <  8:1683235/151‑1 (MQ=255)
                 tcatcaCAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTATCAGCATCATAGCCTTCACGCGCTCGCGac                                                                                                                                  >  2:239048/1‑150 (MQ=255)
                    tcaCAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACtt                                                                                                                               >  7:1507785/1‑151 (MQ=255)
                     cacaGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTc                                                                                                                              >  8:1580624/1‑151 (MQ=255)
                       caGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCt                                                                                                                            >  7:1770187/1‑151 (MQ=255)
                          aTAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGc                                                                                                                         >  5:1674334/1‑151 (MQ=255)
                            acgcAGCTTCGGTGCCAGAATAACAGCAAGCGCAATCCCGATACGCTGGCACTGCCCAGCAGAGAACTCATGCGGATAGCGGTTAATCAGCTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGTCGt                                                                                                                       <  8:571408/148‑1 (MQ=255)
                             aTCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTg                                                                                                                      >  5:427064/1‑151 (MQ=255)
                              tCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTg                                                                                                                      >  8:121571/1‑150 (MQ=255)
                                 gCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACAt                                                                                                                  >  1:860435/1‑151 (MQ=255)
                                 gCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTAACAt                                                                                                                  >  6:527959/1‑151 (MQ=255)
                                    tCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATCGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACAtt                                                                                                                 <  8:681665/149‑1 (MQ=255)
                                      ggTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCg                                                                                                             <  4:1732438/151‑1 (MQ=255)
                                          cAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGgatg                                                                                                         <  6:438043/151‑1 (MQ=255)
                                          cAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGgat                                                                                                          <  4:1465016/150‑1 (MQ=255)
                                          cAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGgat                                                                                                          <  1:1030057/150‑1 (MQ=255)
                                           aaGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGgatga                                                                                                        <  7:670967/151‑1 (MQ=255)
                                             gAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATa                                                                                                      <  3:287653/151‑1 (MQ=255)
                                             gAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATa                                                                                                      <  5:849909/151‑1 (MQ=255)
                                              aaTAAGAGCACGCGCACTCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATa                                                                                                      >  1:1089483/1‑150 (MQ=255)
                                               aTAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATa                                                                                                      >  8:934926/1‑149 (MQ=255)
                                                 aaGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGt                                                                                                   >  4:1028703/1‑150 (MQ=255)
                                                     gCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGc                                                                                               >  3:1731753/1‑150 (MQ=255)
                                                     gCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCa                                                                                              <  3:1085305/151‑1 (MQ=255)
                                                       aCGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGt                                                                                            >  8:1608985/1‑151 (MQ=255)
                                                         gcggAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCTTGCGGATAGCGGGTAATCAGGTTAGGCATTAACCGGACTTTCAGGATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTgg                                                                                          <  2:1003122/151‑1 (MQ=255)
                                                          cgcAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGt                                                                                         >  3:1437423/1‑151 (MQ=255)
                                                           gcAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGtt                                                                                        <  5:1253528/151‑1 (MQ=255)
                                                            cAATCGCGATACGCTGGCACTGCCCACCAGAGAAGTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTc                                                                                       <  2:628815/151‑1 (MQ=255)
                                                            cAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTc                                                                                       <  5:990172/151‑1 (MQ=255)
                                                            cAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTc                                                                                       >  1:68275/1‑151 (MQ=255)
                                                             aaTCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTc                                                                                       >  6:1194849/1‑150 (MQ=255)
                                                              aTCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTc                                                                                       >  5:1505389/1‑149 (MQ=255)
                                                              aTCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAg                                                                                     <  1:1352506/151‑1 (MQ=255)
                                                              aTCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAg                                                                                     <  7:1158853/151‑1 (MQ=255)
                                                              aTCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTATCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACACAGTGGTTc                                                                                       >  4:1417510/1‑149 (MQ=255)
                                                                cccGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCg                                                                                   >  6:1269894/1‑151 (MQ=255)
                                                                   gATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCgatg                                                                                >  1:1394607/1‑151 (MQ=255)
                                                                   gATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCgatg                                                                                <  5:729056/151‑1 (MQ=255)
                                                                      aCGCTGGCACTGCCCACCAGAGCACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCgatgatg                                                                             >  2:1562172/1‑150 (MQ=255)
                                                                        gCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATACCCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGAtctc                                                                           <  1:504007/151‑1 (MQ=255)
                                                                          tGGCCCTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGc                                                                         >  2:56252/1‑151 (MQ=255)
                                                                                  cccACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTc                                                                 <  3:725215/151‑1 (MQ=255)
                                                                                   ccaccaGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCa                                                                <  6:578983/151‑1 (MQ=255)
                                                                                     accaGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATa                                                              <  5:302809/151‑1 (MQ=255)
                                                                                         gagaACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATAcgcg                                                          <  2:1409893/151‑1 (MQ=255)
                                                                                          agaACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCgg                                                         >  3:1669927/1‑151 (MQ=255)
                                                                                            aaCTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGt                                                       >  8:780077/1‑151 (MQ=255)
                                                                                            aaCTCATGCGGATAGCGGGTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCgggg                                                       <  5:1045355/151‑2 (MQ=255)
                                                                                             aCTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGtt                                                      >  8:1258099/1‑151 (MQ=255)
                                                                                                cATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGttt                                                     <  4:82004/149‑1 (MQ=255)
                                                                                                    gggATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGAt                                               >  3:1627916/2‑151 (MQ=255)
                                                                                                    cGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGAt                                               >  5:635259/1‑151 (MQ=255)
                                                                                                      gATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGCTGc                                             >  6:271471/1‑151 (MQ=255)
                                                                                                       aTAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGcc                                            <  3:645494/151‑1 (MQ=255)
                                                                                                          gCGGTTATTCAGGTTAGGCAATATCCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACGTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAAc                                         <  2:369666/151‑1 (MQ=255)
                                                                                                          gCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCCTCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTGCAGCGATGATCTCGACGATGGTCATACGCGGGTTTAGCGATGCCAAc                                         <  8:254119/151‑1 (MQ=255)
                                                                                                          gCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAAc                                         <  2:810394/151‑1 (MQ=255)
                                                                                                           cGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACg                                        <  1:207622/151‑1 (MQ=255)
                                                                                                            ggTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACgg                                       <  1:155310/151‑1 (MQ=255)
                                                                                                            ggTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACgc                                       >  6:1445971/1‑150 (MQ=255)
                                                                                                               tAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATc                                    >  2:1641338/1‑151 (MQ=255)
                                                                                                               tAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATc                                    >  5:1309615/1‑151 (MQ=255)
                                                                                                                aaTCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATcc                                   <  4:741539/151‑1 (MQ=255)
                                                                                                                 aTCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCt                                  <  1:1696074/151‑1 (MQ=255)
                                                                                                                  tCAGGTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTg                                 >  6:347273/1‑151 (MQ=255)
                                                                                                                     ggTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGaa                              <  6:1704215/151‑1 (MQ=255)
                                                                                                                     ggTTAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGaa                              <  7:765145/151‑1 (MQ=255)
                                                                                                                       ttAGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGaaa                             >  7:1718887/1‑150 (MQ=255)
                                                                                                                         aGGCAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAat                           <  7:1451868/150‑1 (MQ=255)
                                                                                                                            cAATAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAatcatc                       >  3:889322/1‑151 (MQ=255)
                                                                                                                             aaTAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCt                      >  5:1154822/1‑151 (MQ=255)
                                                                                                                              aTAACCCGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTg                     <  7:1612894/151‑1 (MQ=255)
                                                                                                                                 aCCCGACTTTCAGCATCATAGCCTTCCCGCGCTCGCGCACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGGTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGCTGCCAACGGATCCTGGAAAATCATCTGAat                  >  8:471081/1‑151 (MQ=255)
                                                                                                                                 aCCCGACTTTCAGCATCTAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAata                 <  7:604778/151‑1 (MQ=255)
                                                                                                                                  cccGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAata                 <  1:781225/151‑1 (MQ=255)
                                                                                                                                   ccGACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAatat                <  2:1739770/151‑1 (MQ=255)
                                                                                                                                     gACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATc               >  8:1057676/1‑150 (MQ=255)
                                                                                                                                     gACTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGGGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCa              <  2:745039/151‑1 (MQ=255)
                                                                                                                                      aCTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCAc             <  1:1214389/151‑1 (MQ=255)
                                                                                                                                      aCTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCAc             >  4:374809/1‑151 (MQ=255)
                                                                                                                                       cTTTCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACt            >  2:1178829/1‑151 (MQ=255)
                                                                                                                                        tttCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTg           >  3:1231911/1‑151 (MQ=255)
                                                                                                                                          tCAGCATCATAGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGc          <  1:362714/150‑1 (MQ=255)
                                                                                                                                               atcataGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGCGAACgg    >  6:92680/1‑151 (MQ=255)
                                                                                                                                               atcataGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGCGAACgg    >  3:42734/1‑151 (MQ=255)
                                                                                                                                                 cataGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGCGAACGGc   <  5:1793177/150‑1 (MQ=255)
                                                                                                                                                 cataGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGCGAACGGCa  <  3:974240/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGCGCCGACACCGGCTCATCACAGATAATCAGCTTCGGTTCAAGAATAAGAGCACGCGCAATCCCGATACGCTGGCACTGCCCACCAGAGAACTCATGCGGATAGCGGTTAATCAGGTTAGGCAATAACCCGACTTTCAGCATCATCGCCTTCACGCGCTCGCGAACTTCCTGGCGTGACATTTTCGGATGATAAGTACGCAGTGGTTCAGCGATGATCTCGCCGATGGTCATACGCGGGTTTAGCGATGCCAACGGATCCTGGAAAATCATCTGAATATCACTGCGAACGGCA  >  DH1/2573328‑2573623

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: