Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 876,146 G→A A560T (GCC→ACC)  unknown → Overcoming lysogenization defect protein Old

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1876,1460GA97.8% 116.2 / NA 46A560T (GCC→ACC) unknownOvercoming lysogenization defect protein Old
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (25/20);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAGGAAATATACCTTTGTACTTGAAATTATTAAATAAATTTAAGATAAAATATATCACTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGGCCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGAGCGGGAGACGTTGGTAGTGGAGCAGGAAAAGAAGAAAGTGGCGCAACTG  >  1/875907‑876384
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
cAAAAGGAAATATACCTTTGTACTTGAAATTATTAAATAAATTTAAGATAAAATATATCACTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAatat                                                                                                                                                                                                                                      <  1:343521/251‑1 (MQ=255)
                                                    atatatCACTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAAt                                                                                                                                                                                                                               >  1:690860/1‑206 (MQ=255)
                                                    atatatCACTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGaa                                                                                                                                                                                  >  1:316246/1‑251 (MQ=255)
                                                         tCACTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAata                                                                                                                                                                             >  1:158416/1‑251 (MQ=255)
                                                            cTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAatatat                                                                                                                                                                                                                                    >  1:622177/1‑193 (MQ=255)
                                                            cTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATc                                                                                                                                                                                                                              <  1:543547/199‑1 (MQ=255)
                                                            cTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                <  1:696056/213‑1 (MQ=255)
                                                            cTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                <  1:90207/213‑1 (MQ=255)
                                                            cTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                >  1:205532/1‑213 (MQ=255)
                                                              atatatGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                >  1:567158/1‑211 (MQ=255)
                                                              atatatGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCa                                                                                                                                                                        >  1:78526/1‑251 (MQ=255)
                                                                atatGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATc                                                                                                                                                                                                                              <  1:111277/195‑1 (MQ=255)
                                                                     aTCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAg                                                                                                                                                                                                      <  1:839557/214‑1 (MQ=255)
                                                                                                    tGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTg                                                                                                                                                     <  1:255613/232‑1 (MQ=255)
                                                                                                    tGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGa                                                                                                                                  >  1:228817/1‑251 (MQ=255)
                                                                                                     gATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATc                                                                                                                                                                                                                              >  1:528477/1‑158 (MQ=255)
                                                                                                                  ccACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                >  1:228990/1‑159 (MQ=255)
                                                                                                                                     aaTTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCa                                                                                                 >  1:470947/1‑251 (MQ=255)
                                                                                                                                                ctATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAatatat                                                                                                                                                                                                                                    >  1:877469/1‑109 (MQ=255)
                                                                                                                                                ctATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAAt                                                                                                                                                                                                                               <  1:198419/114‑1 (MQ=255)
                                                                                                                                                         aagaGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATca‑c                                                                              <  1:178285/249‑1 (MQ=255)
                                                                                                                                                         aagaGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATCATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATca‑c                                                                              <  1:492038/249‑1 (MQ=255)
                                                                                                                                                           gagTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAAt                                                                                                                                                                                                                               >  1:663870/1‑103 (MQ=255)
                                                                                                                                                           gagTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATc                                                                                                                                                                                                                              >  1:51010/1‑104 (MQ=255)
                                                                                                                                                           gagTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAAc                                                                                                     <  1:738560/225‑1 (MQ=255)
                                                                                                                                                                      ttCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGGCCTTGCAATATATaaaa                                                                                                                                                                                                                                <  1:43821/91‑1 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAatat                                                                                                                                                                                                                                      >  1:701114/1‑83 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATc                                                                                                                                                                                                                              <  1:699601/91‑1 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATc                                                                                                                                                                                                                              >  1:772040/1‑91 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                >  1:246233/1‑105 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATaa                                                                                                                                                        >  1:92203/1‑161 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATaa                                                                                                                                                        >  1:581275/1‑161 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATca‑c                                                                              >  1:864164/1‑234 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATca‑c                                                                              <  1:126565/234‑1 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATc                                                                                 >  1:662375/1‑232 (MQ=255)
                                                                                                                                                                        cTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTAt                                                              >  1:855876/1‑250 (MQ=255)
                                                                                                                                                                                   ataAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAac                                                                                             <  1:308631/209‑1 (MQ=255)
                                                                                                                                                                                                         ttGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGAGCGGGAGACGTTGGTAGTGGAGc                            <  1:58605/251‑1 (MQ=255)
                                                                                                                                                                                                           ggTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATca‑ccac                                                                           <  1:205227/202‑1 (MQ=255)
                                                                                                                                                                                                            gTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAAc                                                                                                                                                                                                                >  1:368965/1‑69 (MQ=255)
                                                                                                                                                                                                                     gtgGTAAATCAGGGAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATc                                                                                                                                                                         <  1:248519/99‑1 (MQ=255)
                                                                                                                                                                                                                                 ggAAGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGAGCGGGAGACGTTGGTAGTGGAGCAGGAAAAGAAGAAAGTGgcgc       <  1:827268/248‑1 (MQ=255)
                                                                                                                                                                                                                                    aGCCTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGAGCGGGAGACGTTGGTAGTGGAGCAGGAAAAGAAGAAAGTGGCGCAACTg  >  1:869434/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cTTTTCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGa                                                   >  1:252756/1‑198 (MQ=255)
                                                                                                                                                                                                                                        ttttCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCACCCACTCAGGGCAGTTATCTGCATAAGGAGCGGGAGACGTTGGt                                    <  1:27033/213‑1 (MQ=255)
                                                                                                                                                                                                                                        ttttCAGACCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGAGCGGgag                                            <  1:413455/204‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
CAAAAGGAAATATACCTTTGTACTTGAAATTATTAAATAAATTTAAGATAAAATATATCACTATATATGATCAAGATCATCAAGATAGAAAATCTGAAAATGATAAAGCTGTTGCCACTAAACTGACTAATAAAATTTTATCTTCTATTGAAGAAGAGTTAGGTCATTCTATTATGCTTATAAATGACATTGAAGAAGAGCTTGGTTATGAATGTGGTAAATCAGGGAAGCCTTTTCAGGCCTTGCAATATATAAAATCTGATAGTTTTAAACTATCCCCAAGTTTTGAAGAAAAAATCAGAAAAATATATCAATAATCAAAATTATAATTGCTAGTTAATTGAAGTAAGAAAAATTTTAGGAGCCACCCATGCGGAAACTGCAAAACACGCTCTATATCA‑CCACTCAGGGCAGTTATCTGCATAAGGAGCGGGAGACGTTGGTAGTGGAGCAGGAAAAGAAGAAAGTGGCGCAACTG  >  1/875907‑876384

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