Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,581,314 C→T E237K (GAA→AAA)  unknown ← Guanosine‑3',5'‑bis(diphosphate) 3'‑pyrophosphohydrolase (EC 3.1.7.2) / GTP pyrophosphokinase (EC 2.7.6.5), (p)ppGpp synthetase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,581,3140CT100.0% 219.8 / NA 74E237K (GAA→AAA) unknownGuanosine‑3',5'‑bis(diphosphate) 3'‑pyrophosphohydrolase (EC 3.1.7.2) / GTP pyrophosphokinase (EC 2.7.6.5), (p)ppGpp synthetase II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (46/28);  total (46/28)

GTGCGGACCAATCATCGACGTATGTAAGGACTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAGAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTCACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTGGCGATACGACGG  >  1/1581081‑1581558
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
gTGCGGACCAATCATCGACGTATGTAAGGACTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCcgcg                                                                                                                                                                                                                                     >  1:719484/1‑251 (MQ=255)
                     aTGTAAGGACTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCa                                                                                                                                                                                                                >  1:476598/1‑251 (MQ=255)
                      tGTAAGGACTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACAt                                                                                                                                                                                                                      <  1:700758/244‑1 (MQ=255)
                      tGTAAGGACTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCaa                                                                                                                                                                                                               >  1:462924/1‑251 (MQ=255)
                             aCTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGcc                                                                                                                                                                                                          >  1:99793/1‑249 (MQ=255)
                              cTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCa                                                                                                                                                                                                                <  1:704302/242‑1 (MQ=255)
                                  tAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGAt                                                                                                                                                                                                   <  1:885516/251‑1 (MQ=255)
                                          ggCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCgg                                                                                                                                                                                           >  1:850191/1‑251 (MQ=255)
                                                       gcaATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATAc                                                                                                                                                                                                                            >  1:121963/1‑205 (MQ=255)
                                                       gcaATATAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACAt                                                                                                                                                                                                                      >  1:346955/1‑211 (MQ=255)
                                                          atatAATCTTTTACCCGACCCGGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATACTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCAt                                                                                                                                                                           >  1:435034/1‑251 (MQ=255)
                                                                           acccgGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTc                                                                                                                                                                        <  1:176756/237‑1 (MQ=255)
                                                                           acccgGGCGCGGTTTATAAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACcg                                                                                                                                                           >  1:550091/1‑250 (MQ=255)
                                                                                         tataAAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACg                                                                                                                                            >  1:597052/1‑251 (MQ=255)
                                                                                               ggCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACtttt                                                                                                                                      >  1:574226/1‑251 (MQ=255)
                                                                                                gCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTc                                                                                                                                                                        >  1:518560/1‑216 (MQ=255)
                                                                                                gCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTAc                                                                                                                                                             <  1:566737/227‑1 (MQ=255)
                                                                                                     gCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACtt                                                                                                                                        >  1:600424/1‑243 (MQ=255)
                                                                                                     gCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCaaa                                                                                                                                >  1:81370/1‑251 (MQ=255)
                                                                                                     gCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTCTTGGCAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCTTTCAATCATGTCTTTGCGATTACCCCCCGCAATCTGACCGCCTTTTGGCaaa                                                                                                                                >  1:685426/1‑251 (MQ=255)
                                                                                                       aTTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAAc                                                                                                                              <  1:565457/251‑1 (MQ=255)
                                                                                                                 gCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAAc                                                                                                                      <  1:332959/249‑1 (MQ=255)
                                                                                                                   acacGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCa                                                                                                                                                                                                                >  1:610913/1‑157 (MQ=255)
                                                                                                                    cacGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATcc                                                                                                                                                                                                                    >  1:731225/1‑152 (MQ=255)
                                                                                                                    cacGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTg                                                                                                                  >  1:55834/1‑250 (MQ=255)
                                                                                                                    cacGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTg                                                                                                                  >  1:156753/1‑250 (MQ=255)
                                                                                                                      cGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTAc                                                                                                                                                             <  1:597928/205‑1 (MQ=255)
                                                                                                                      cGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTa                                                                                                                                              <  1:328831/220‑1 (MQ=255)
                                                                                                                      cGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCgg                                                                                                               >  1:795504/1‑251 (MQ=255)
                                                                                                                              ggTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATgg                                                                                                       >  1:129726/1‑251 (MQ=255)
                                                                                                                               gTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATc                                                                                                                                                  <  1:830558/207‑1 (MQ=255)
                                                                                                                               gTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATgg                                                                                                       <  1:43460/250‑1 (MQ=255)
                                                                                                                                  tCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGAt                                                                                                                                                                >  1:415989/1‑190 (MQ=255)
                                                                                                                                          gttCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATAc                                                                                                                                                                                                                            >  1:292577/1‑122 (MQ=255)
                                                                                                                                          gttCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTAc                                                                                                                                                             <  1:20843/185‑1 (MQ=255)
                                                                                                                                          gttCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCaa                                                                                           >  1:62591/1‑251 (MQ=255)
                                                                                                                                            tCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGa                                                                                         >  1:16496/1‑251 (MQ=255)
                                                                                                                                                  cGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTcc                                                                                   <  1:253913/251‑1 (MQ=255)
                                                                                                                                                    gTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGa                                                                                                              >  1:471004/1‑222 (MQ=255)
                                                                                                                                                    gTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAg                                                                                 >  1:400510/1‑251 (MQ=255)
                                                                                                                                                     tGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATGCATGGCTTCAAAACCAAGATCTTCCAg                                                                                 <  1:741625/250‑1 (MQ=255)
                                                                                                                                                            cgcgTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTc                                                                                                                                                                        >  1:846325/1‑156 (MQ=255)
                                                                                                                                                             gcgTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATc                                                                                       >  1:110513/1‑236 (MQ=255)
                                                                                                                                                             gcgTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGttttt                                                                        >  1:80084/1‑251 (MQ=255)
                                                                                                                                                             gcgTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGttttt                                                                        >  1:856662/1‑251 (MQ=255)
                                                                                                                                                              cgTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACa                                                                                                          >  1:56064/1‑216 (MQ=255)
                                                                                                                                                              cgTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGtt                                                                               >  1:398731/1‑243 (MQ=255)
                                                                                                                                                                tAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGAc                                                                                                                                          >  1:830507/1‑182 (MQ=255)
                                                                                                                                                                tAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGAt                                                                     >  1:466963/1‑251 (MQ=255)
                                                                                                                                                                     atCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATTCCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATc                                                                                                                                                  <  1:118447/169‑1 (MQ=255)
                                                                                                                                                                       ccATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTa                                                                                                                                              <  1:662275/171‑1 (MQ=255)
                                                                                                                                                                       ccATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGt                                                                                                                    <  1:853690/197‑1 (MQ=255)
                                                                                                                                                                          tAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAg                                                                                 >  1:215444/1‑229 (MQ=255)
                                                                                                                                                                            aTGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATcc                                                          >  1:706129/1‑250 (MQ=255)
                                                                                                                                                                              ggAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCa                                                       >  1:257088/1‑251 (MQ=255)
                                                                                                                                                                               gAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCaa                                                      >  1:674996/1‑251 (MQ=255)
                                                                                                                                                                               gAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCaa                                                      >  1:357349/1‑251 (MQ=255)
                                                                                                                                                                               gAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGC‑AACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGAtg                                                                    >  1:336593/1‑236 (MQ=255)
                                                                                                                                                                                aaTGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAg                                                     <  1:48880/251‑1 (MQ=255)
                                                                                                                                                                                aaTGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAg                                                     <  1:306005/251‑1 (MQ=255)
                                                                                                                                                                                 aTGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAg                                                     >  1:251912/1‑250 (MQ=255)
                                                                                                                                                                                 aTGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGc                                                    <  1:322139/251‑1 (MQ=255)
                                                                                                                                                                                             aTCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAAc                                         <  1:416900/250‑1 (MQ=255)
                                                                                                                                                                                              tCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACgg                                       <  1:409065/251‑1 (MQ=255)
                                                                                                                                                                                                 tttAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACggcg                                     >  1:901830/1‑248 (MQ=255)
                                                                                                                                                                                                       cgcATTTTTTGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCCTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATaaa                               <  1:451715/250‑1 (MQ=255)
                                                                                                                                                                                                             ttttGGTAGATGGCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAg                                                     <  1:222740/222‑1 (MQ=255)
                                                                                                                                                                                                                         gCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTgg             >  1:732024/1‑250 (MQ=255)
                                                                                                                                                                                                                         gCATACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTgg             >  1:224220/1‑250 (MQ=255)
                                                                                                                                                                                                                           aTACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCaaaa                                                                                               <  1:726862/166‑1 (MQ=255)
                                                                                                                                                                                                                           aTACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGAtg                                                                    <  1:636516/193‑1 (MQ=255)
                                                                                                                                                                                                                            tACAAATGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTGGCGAt         >  1:524674/1‑251 (MQ=255)
                                                                                                                                                                                                                                aaTGCTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTGGCGATacga     <  1:371695/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cTTTTTACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTGGCGATACGACgg  <  1:559403/250‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
GTGCGGACCAATCATCGACGTATGTAAGGACTGATAGCCATTGGCTTTCGGCACGGCAATATAATCTTTTACCCGACCCGGGCGCGGTTTATAGAGGCTGTGCATTTGTCCAAGCACACGATAGCAGGTATCCACGTTGTTCACCACGGTGCGGAACGCGTAAATATCCATAATGGAATGAAACTGCTGATCTTTTAAGCGCATTTTTTGGTAGATGGCATACAAATGCTTTTCACGCCCGAATACCCGCGCTTGAATACCCACATCCTGCAAACGCCCTTTGATTTCGTCGGAAATTCGTTCAATCATGTCTTTGCGATTACCGCGCGCAATCTGTACGACTTTTTGCAAAACCGCATAACGTTGCGGATACATGGCTTCAAAACCAAGATCTTCCAGTTCGTTTTTGATGTGCTCGATCCCCAAGCGATGCGCCAACGGTGCATAAATTTCCAGGGTTTCTTTGGCGATACGACGG  >  1/1581081‑1581558

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑