Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 444,535 A→C 16.2% L130V (TTG→GTG)  hdfR_1 ← HTH‑type transcriptional regulator HdfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122444,5350AC16.2% 114.9 / 3.5 67L130V (TTG→GTG) hdfR_1HTH‑type transcriptional regulator HdfR
Reads supporting (aligned to +/- strand):  ref base A (39/15);  new base C (2/11);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 2.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAGCGACTGGCCCTTG  >  ACB122/444406‑444668
                                                                                                                                 |                                                                                                                                     
gggtATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  <  8:70404/135‑1 (MQ=255)
gggtATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  <  3:335093/135‑1 (MQ=255)
gggtATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  >  4:211702/1‑135 (MQ=255)
        ggTATGACCAGTTGCAGCTCTTCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCg                                                                                                                          <  3:267153/135‑1 (MQ=255)
           aTGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTaa                                                                                                                       >  7:268321/1‑135 (MQ=255)
               ccAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGCTCTCCAACCCCGACAACTCGTAAATCt                                                                                                                   >  5:263895/1‑135 (MQ=255)
                cAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTg                                                                                                                  >  1:22384/1‑135 (MQ=255)
                  gTTGCAGCTCATCTTCGAACAACAGCAGGCCATGCAGCCCTGGGTGGCGCGGCGGAAAATAGCTGATGCCGATGTCCTGGTTGCCTTTGAGCAGGCGCCGTTCTATCACCAACCCCGACAACTCGTAAATCTCCa                                                                                                                <  6:30862/135‑1 (MQ=255)
                     gCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACACCTCGTAAATCTGCACGa                                                                                                             >  3:289454/1‑135 (MQ=255)
                      cAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCCACCCCGACAACTCGTAAATCTGCACGa                                                                                                             >  8:305295/1‑134 (MQ=255)
                      cAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCCACCCCGACAACTCGTAAATCTGCACGa                                                                                                             >  6:102725/1‑134 (MQ=255)
                      cAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGa                                                                                                             >  4:121032/1‑134 (MQ=255)
                      cAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGa                                                                                                             >  5:403584/1‑134 (MQ=255)
                        tcTCATCTTCGTACATACGCAGGCCATGCAGCCCTGGCTGGCGCGGCGTCAAATAGCTTATGCCTATTTCCAGGTTTCCTTTTAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGAACGACCa                                                                                                          <  8:332713/134‑1 (MQ=255)
                                tCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGt                                                                                                  >  7:227607/1‑135 (MQ=255)
                                 cGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGt                                                                                                  >  8:305607/1‑134 (MQ=255)
                                  gTACACCAGCCGGCCAGGCAGCCGTGGCTGGCGCGGCGGAAAATCGCTGAGGCCGAGGTCCTGGTGGCCGTGGAGCAGGCGCGGTTCGAGCTCCACCCCCCACAACTCGTAAATCGGCACGACCAGGTGAGGtt                                                                                                 <  2:247028/134‑1 (MQ=255)
                                      aacaGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCACCCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGcc                                                                                            >  1:133832/1‑135 (MQ=255)
                                           cagGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTgc                                                                                        >  8:391437/1‑134 (MQ=255)
                                              gCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAg                                                                                    >  7:460075/1‑135 (MQ=255)
                                              gCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAg                                                                                    >  3:62228/1‑135 (MQ=255)
                                              gCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAg                                                                                    >  8:259723/1‑135 (MQ=255)
                                              gCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAg                                                                                    >  1:24136/1‑135 (MQ=255)
                                                cATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  <  3:167210/87‑1 (MQ=255)
                                                cATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  <  3:159774/87‑1 (MQ=255)
                                                cATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  >  4:159774/1‑87 (MQ=255)
                                                cATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCg                                                                                                                                  >  4:167210/1‑87 (MQ=255)
                                                 aTGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCg                                                                                 <  7:259723/135‑1 (MQ=255)
                                                   gCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTTCAGGTAGCCGTTGAGCAGGCGCAGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGAGAGGTTGAGCCTTGCGCAGACGt                                                                                >  5:148139/1‑134 (MQ=255)
                                                     aGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGtt                                                                               >  7:84302/1‑133 (MQ=255)
                                                           ggctggcGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAagccgc                                                                       >  2:238720/1‑135 (MQ=255)
                                                            gctggcGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCt                                                                      >  1:155892/1‑135 (MQ=255)
                                                               ggcGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCACCCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCg                                                                   >  3:263081/1‑135 (MQ=255)
                                                               ggcGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCg                                                                   >  1:88951/1‑135 (MQ=255)
                                                                gcgcGGCGGCAAATCGCTGAGGCCGATGCCCAGGTGGCCGTTGAGCAGGCGCCGTTCGATCCCCACCCCCCACAACTCGTAAATCGGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCgg                                                                  <  8:122072/135‑1 (MQ=255)
                                                                 cgcgGCGCAAAATCGCTGATGCCGATGTCCAGGTGGCGGTTGAGCTGGCGCCGTTCGATCTCCAACCCCGACACCTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCggg                                                                 <  4:396248/135‑1 (MQ=255)
                                                                    ggcggcAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGact                                                              >  4:287636/1‑133 (MQ=255)
                                                                    ggcggcAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGa                                                                >  3:355890/1‑133 (MQ=255)
                                                                    ggcggcAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGa                                                                >  6:351973/1‑133 (MQ=255)
                                                                      cggaAAATCGCTGAGGCCGTTGTCCTGGTTGCCGTTGAGCAGCCGCCGTGCGATCTCCACCCCCGACAACTCGAAAATCTCCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAg                                                            <  6:403584/135‑1 (MQ=255)
                                                                         cAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCt                                                         >  8:400835/1‑135 (MQ=255)
                                                                         cAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGAGCAGGCt                                                         >  2:200954/1‑135 (MQ=255)
                                                                            aTAGCTGATGCCGAGGTCCAGGTTGCCGTGGAGCCGGCGCCGTGCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTgg                                                      <  1:39522/135‑1 (MQ=255)
                                                                              aGCTGATGCGGTTGTCCTGGTTGCGTTTGAGCCGGCGCGTTGGGATCTCCAACCCCCACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGcc                                                    <  4:454743/135‑1 (MQ=255)
                                                                              aGCTGATGCCGTTGTCCTGGGTGCCGTTGCGCCGGCGCCGTTCGATCCCCACCCCCGACAACTCGTAAATCGCCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGcc                                                    <  1:195144/135‑1 (MQ=255)
                                                                              aGCTGATGCCGATGTCCAGGGTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGcc                                                    <  4:214472/135‑1 (MQ=255)
                                                                              aGCTGAGGCCGAGGCCCCGGTGGCCGTTGCGCTGGCGCCGGGCTACCCCAACCCCCGACAACTCGAAACTCTCCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGcc                                                    <  4:83183/135‑1 (MQ=255)
                                                                              aGCTCATGCCGAGGGCCCGGTTGCCTTTGAGCAGGCGCCGTTCGATCTCCAACCCCCACACCTCGTAAAGCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGcc                                                    <  3:121032/135‑1 (MQ=255)
                                                                                 tGATGCCGATGGCCCGGGGGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCCACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGc                                                  <  2:15092/134‑1 (MQ=255)
                                                                                    tGCCGATGGCCCGGTGGCCGTTGAGCCGGCGCCGGTCGATCTCCACCCCCGCCAACGCGAAAATCCGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCAc                                                <  4:79356/133‑1 (MQ=255)
                                                                                    tGCCGAGGTCCAGGTGGCCGTTGAGCCGGCGCCGTTCGAGCACCACCCCCGACACCTCGAAAATCGGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCAc                                                <  7:78016/133‑1 (MQ=255)
                                                                                    tGCCGAGGTCCAGGTGGCCGTTGAGCAGGCGCGGTTCGATCCCCACCCCCGACAACTCGAAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCAc                                                <  1:52251/133‑1 (MQ=255)
                                                                                      ccaaTGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTc                                            >  5:109622/4‑135 (MQ=255)
                                                                                      ccGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTc                                            >  8:67878/1‑135 (MQ=255)
                                                                                      ccGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTc                                            >  6:271959/1‑135 (MQ=255)
                                                                                      ccGATGTCCAGCTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCAGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTc                                            >  4:421393/1‑135 (MQ=255)
                                                                                          tgccctggTGGCCGTTGCGCCGGCGCCGTGCGAGCACAACCCCCCACAACTCGAAAATCGCCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCa                                        <  5:102725/129‑1 (MQ=255)
                                                                                             ccAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAAt                                      >  5:99238/1‑134 (MQ=255)
                                                                                                       gTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAg                            >  6:174171/1‑134 (MQ=255)
                                                                                                       gTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAg                            >  4:16953/1‑134 (MQ=255)
                                                                                                        ttGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGcc                          <  6:157543/135‑1 (MQ=255)
                                                                                                         tGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGc                           >  5:87063/1‑133 (MQ=255)
                                                                                                         tGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGc                           >  4:396660/1‑133 (MQ=255)
                                                                                                                   ccGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAg               >  8:205631/1‑135 (MQ=255)
                                                                                                                      ttCGATCTCCACCCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGGTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAGCGa            <  6:202349/135‑1 (MQ=255)
                                                                                                                        cGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAGCGACt          >  7:386740/1‑135 (MQ=255)
                                                                                                                           tccccACCCCCGACAACTCGTAAAGCTGCACGACCAGGTGAGGGTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAGCGACTGGc       <  4:263081/132‑1 (MQ=255)
                                                                                                                                caCCCCCCCAACTCCTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAGCGACTGGCCCTTg  <  3:38979/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
GGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGGCGCGGCGGCAAATAGCTGATGCCGATGTCCAGGTTGCCGTTGAGCAGGCGCCGTTCGATCTCCAACCCCGACAACTCGTAAATCTGCACGACCAGGTGAGGTTGAGCCTTGCGCAGCCGTTCAAGCAGCTGCGGGATCAGGCTTGGCCGCACGGTCTGCAATACACCGATAGCCAAGGTGCGCAGCGACTGGCCCTTG  >  ACB122/444406‑444668

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