Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 668,994 T→G E575A (GAA→GCA)  rpoD ← RNA polymerase sigma factor RpoD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1668,9940TG95.0% 1287.3 / 35.9 343E575A (GAA→GCA) rpoDRNA polymerase sigma factor RpoD
Reads supporting (aligned to +/- strand):  ref base T (5/12);  new base G (122/204);  total (127/216)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 350 total aligned reads displayed.

GTGCGGCGTAACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTCGAGGGTG  >  DH1/668845‑669140
                                                                                                                                                     |                                                                                                                                                  
gTGCGGCGTAACGCCTGATCCGGCCTACCGATTAATCGTAC‑GGAAGCTACGCAGCACTTCAGAACGTCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCActgct                                                                                                                                                  >  3:744196/1‑151 (MQ=255)
     gcgTAACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCa                                                                                                                                              <  6:328451/151‑1 (MQ=255)
      ggTAACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAg                                                                                                                                             <  3:946606/150‑1 (MQ=255)
      cgTCACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAg                                                                                                                                             <  8:1355064/151‑1 (MQ=255)
      cgTAACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAg                                                                                                                                             <  6:1464228/151‑1 (MQ=255)
       gTAACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGc                                                                                                                                            <  3:348412/151‑1 (MQ=255)
         accGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCgt                                                                                                                                          >  6:180474/3‑151 (MQ=255)
         aaCGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCgt                                                                                                                                          <  7:792095/151‑1 (MQ=255)
          aCGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCgtg                                                                                                                                         <  1:1057149/151‑1 (MQ=255)
             ccTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCgtgt                                                                                                                                        >  2:1626210/1‑149 (MQ=255)
             ccTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCCCGTAACGTCGAACTGTTTACCCACTGCTTCCAGCgtgt                                                                                                                                        >  4:442676/1‑149 (MQ=255)
                gATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCg                                                                                                                                   <  5:1855007/151‑1 (MQ=255)
                 aTCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCgg                                                                                                                                  <  6:787497/151‑1 (MQ=255)
                     ggCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGtt                                                                                                                              <  6:276875/151‑1 (MQ=255)
                     ggCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGtt                                                                                                                              <  4:657642/151‑1 (MQ=255)
                      gCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTc                                                                                                                             <  1:811866/151‑1 (MQ=255)
                          aCCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCatat                                                                                                                         <  8:757593/151‑1 (MQ=255)
                          aCCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCat                                                                                                                           <  2:1492173/149‑1 (MQ=255)
                            cGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCg                                                                                                                       <  3:1675919/151‑1 (MQ=255)
                            cGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCg                                                                                                                       <  2:755654/151‑1 (MQ=255)
                             gATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGa                                                                                                                      <  7:1665909/151‑1 (MQ=255)
                              aTTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGAt                                                                                                                     >  2:973198/1‑151 (MQ=255)
                                 aaTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATAcc                                                                                                                  <  3:1150903/151‑1 (MQ=255)
                                 aaTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATAcc                                                                                                                  <  2:376580/151‑1 (MQ=255)
                                 aaTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATAcc                                                                                                                  <  7:90643/151‑1 (MQ=255)
                                 aaTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATAc                                                                                                                   <  7:159785/150‑1 (MQ=255)
                                 aaTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATAc                                                                                                                   <  6:1037202/150‑1 (MQ=255)
                                  aTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCg                                                                                                                 <  8:367130/151‑1 (MQ=255)
                                  aTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCg                                                                                                                 >  8:1805930/1‑151 (MQ=255)
                                  aTCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCg                                                                                                                 <  2:1127908/151‑1 (MQ=255)
                                    cGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGaa                                                                                                               <  3:1095675/151‑1 (MQ=255)
                                     gTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGaaa                                                                                                              <  4:1531615/151‑1 (MQ=255)
                                     gTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGaa                                                                                                               >  8:1228375/1‑150 (MQ=255)
                                       ccAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACg                                                                                                            <  2:512251/151‑1 (MQ=255)
                                       ccAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACg                                                                                                            <  5:887273/151‑1 (MQ=255)
                                         aGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGc                                                                                                           >  6:92599/1‑150 (MQ=255)
                                         aGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCa                                                                                                          <  4:925177/151‑1 (MQ=255)
                                          ggAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCAt                                                                                                         <  2:942984/151‑1 (MQ=255)
                                            aaGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATAc                                                                                                       <  7:531203/151‑1 (MQ=255)
                                                    cagcaCTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGAAACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACt                                                                                               >  2:703182/1‑151 (MQ=255)
                                                     agcaCTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACtt                                                                                              <  1:1421097/151‑1 (MQ=255)
                                                      gcaCTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACttt                                                                                             <  3:1501635/151‑1 (MQ=255)
                                                        aCTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTg                                                                                           >  1:641955/1‑151 (MQ=255)
                                                        aCTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTg                                                                                           >  5:119990/1‑151 (MQ=255)
                                                           tCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGc                                                                                          <  3:1348019/149‑1 (MQ=255)
                                                            ggatctggCTCGGGTGACGCAGTTTGCGCATCGACTTCGCTTCGATCTTCCGGATACGTTCGCGGGTAACGTCGAACTTGGTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTc                                                                                       <  1:1609623/145‑1 (MQ=255)
                                                            cAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTc                                                                                       >  6:1221922/1‑151 (MQ=255)
                                                               aaCGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCAcg                                                                                    >  4:1607603/1‑151 (MQ=255)
                                                                 cGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCAcgcg                                                                                  <  5:526898/151‑1 (MQ=255)
                                                                 cGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCAcgcg                                                                                  >  7:27511/1‑151 (MQ=255)
                                                                     tCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGt                                                                              >  7:1137351/1‑151 (MQ=255)
                                                                          ggACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGgc                                                                         <  5:369471/150‑1 (MQ=255)
                                                                           gACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTGTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGgcc                                                                        <  5:92394/151‑1 (MQ=255)
                                                                               cAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGgccagcc                                                                    >  6:53204/1‑151 (MQ=255)
                                                                                 gTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGccagccag                                                                  <  3:247556/151‑1 (MQ=255)
                                                                                   ttGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCcagccagc                                                                 <  4:1696038/150‑1 (MQ=255)
                                                                                   ttGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCcagccagc                                                                 >  1:1316079/1‑150 (MQ=255)
                                                                                    tGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCAc                                                               >  5:1753483/1‑151 (MQ=255)
                                                                                     gcgcAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCAc                                                               <  6:401127/150‑1 (MQ=255)
                                                                                     gcgcAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCAc                                                               <  5:1369061/150‑1 (MQ=255)
                                                                                       gcAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCAcgtc                                                            <  7:1709145/151‑1 (MQ=255)
                                                                                         aGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCAcgtcgt                                                          <  2:844740/151‑1 (MQ=255)
                                                                                          gcgcCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTg                                                         <  5:613979/151‑1 (MQ=255)
                                                                                          gcgcCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACCCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTg                                                         <  8:1742504/151‑1 (MQ=255)
                                                                                            gcCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCg                                                       >  3:499511/1‑151 (MQ=255)
                                                                                                tcgcttcgATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGtt                                                     <  7:1529446/149‑1 (MQ=255)
                                                                                                    ttcgATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCcgc                                                >  1:981379/1‑150 (MQ=255)
                                                                                                     tcgATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCcgca                                               <  8:1493101/150‑1 (MQ=255)
                                                                                                       gATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCcgcacgc                                            >  5:158001/1‑151 (MQ=255)
                                                                                                        aTCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCcgcacgc                                            >  8:422548/1‑150 (MQ=255)
                                                                                                          cTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAgg                                         >  3:483710/1‑151 (MQ=255)
                                                                                                            gACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCcgcacgca                                           <  4:210757/147‑1 (MQ=255)
                                                                                                            gACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCcgcacgca                                           <  7:667828/147‑1 (MQ=255)
                                                                                                            gACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCt                                       <  4:741918/151‑1 (MQ=255)
                                                                                                            gACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCt                                       <  1:223793/151‑1 (MQ=255)
                                                                                                            gACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCt                                       <  5:38012/151‑1 (MQ=255)
                                                                                                             aCGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCtt                                      >  2:261608/1‑151 (MQ=255)
                                                                                                             aCGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCtt                                      <  4:1210875/151‑1 (MQ=255)
                                                                                                              cGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCttt                                     <  6:1555482/151‑1 (MQ=255)
                                                                                                              cGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCttt                                     <  6:1121643/151‑1 (MQ=255)
                                                                                                              cGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCtat                                     >  4:601277/1‑149 (MQ=255)
                                                                                                              cGGATACGGTCGCGGGTCACGTCGAACTGGTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATAGCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACCTCGTGCGTTGCCGCACGCAGGCttt                                     <  1:453646/151‑1 (MQ=255)
                                                                                                                gATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCg                                   <  1:1133681/151‑1 (MQ=255)
                                                                                                                       tCGCGGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCg                            >  8:1060789/1‑151 (MQ=255)
                                                                                                                        cgcgGGTAACGTCGAACTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCg                            <  6:1638520/150‑1 (MQ=255)
                                                                                                                          cgGGTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAg                         <  4:398368/151‑1 (MQ=255)
                                                                                                                           gggTAACGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCAGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGa                        <  1:701286/151‑1 (MQ=255)
                                                                                                                               aaCGTCGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATcc                    <  2:669953/151‑1 (MQ=255)
                                                                                                                                    cGAACTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCgg               <  2:1237462/151‑1 (MQ=255)
                                                                                                                                      aaCTGTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCa             <  1:122531/151‑1 (MQ=255)
                                                                                                                                      aaCTGTGTACCCACTGCTTCCAGCGTGTAGTCGGGGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCa             <  7:1296632/151‑1 (MQ=255)
                                                                                                                                        cTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGc           >  1:563416/1‑151 (MQ=255)
                                                                                                                                        cTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGc           >  7:108274/1‑151 (MQ=255)
                                                                                                                                          ggttACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTc         <  8:629205/149‑1 (MQ=255)
                                                                                                                                          gTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGc           <  3:1441323/149‑1 (MQ=255)
                                                                                                                                          gTTTACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTc         >  1:369098/1‑151 (MQ=255)
                                                                                                                                           tttACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTCg        >  1:17358/1‑151 (MQ=255)
                                                                                                                                            ttACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTCg        <  5:112716/150‑1 (MQ=255)
                                                                                                                                             tACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTCg        >  8:101799/1‑149 (MQ=255)
                                                                                                                                             tACCCACTGCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTCGa       <  4:1412636/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGCGGCGTAACGCCTGATCCGGCCTACCGATTAATCGTCCAGGAAGCTACGCAGCACTTCAGAACGGCTCGGGTGACGCAGTTTGCGCAGCGCCTTCGCTTCGATCTGACGGATACGTTCGCGGGTAACGTCGAACTGTTTACCCACTTCTTCCAGCGTGTAGTCGGTGTTCATATCGATACCGAAACGCATACGCAGAACTTTTGCTTCACGCGCGGTCAGGCCAGCCAGCACGTCGTGCGTTGCCGCACGCAGGCTTTCGGTGGTCGCAGAATCCAGCGGCAGCTCGAGGGTG  >  DH1/668845‑669140

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