Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,021,185:1 +T 5.8% coding (122/330 nt) yxiH ← conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,021,1851.T 5.8% 228.6 / 11.2 69coding (122/330 nt)yxiHconserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (39/26);  new base T (2/2);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AAAACCTTTACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGTTATGCTC  >  BS168_reference_genome/4021047‑4021317
                                                                                                                                           |                                                                                                                                    
aaaaCCTTTACCTTCGTCTCTTACGATTCGTAATAAACTAAGACATAGCCATCTGCACATTGTTGAAAGAGAAGAATATAACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGa                                                                                                                                       <  5:83350/139‑1 (MQ=255)
aaaaCCTTTACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGa                                                                                                                                       <  2:224211/139‑1 (MQ=255)
aaaaCCTTTACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGa                                                                                                                                       <  5:186059/139‑1 (MQ=255)
aaaaCCTTTACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGa                                                                                                                                       <  2:47900/139‑1 (MQ=255)
        tACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAAt                                                                                                                               <  5:279832/138‑1 (MQ=255)
        tACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATg                                                                                                                              >  4:208710/1‑139 (MQ=255)
               ggCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTTTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATATATTGTAATGGTGATa                                                                                                                        >  5:354821/1‑138 (MQ=255)
               ggCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATac                                                                                                                       >  6:32470/1‑139 (MQ=255)
               ggCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATac                                                                                                                       <  5:42789/139‑1 (MQ=255)
               ggCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATac                                                                                                                       <  2:55379/139‑1 (MQ=255)
               ggCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATac                                                                                                                       <  4:338896/139‑1 (MQ=255)
                gCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATacc                                                                                                                      <  6:307281/139‑1 (MQ=255)
                     tACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCCCCtt                                                                                                                 >  5:371367/1‑139 (MQ=255)
                     tACGATTCGTAATAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTTGTAATGTTGATACCACCtt                                                                                                                 >  3:30977/1‑139 (MQ=255)
                     tACGATTCGTAATAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTTGTAATGGTGATACCACCtt                                                                                                                 >  3:36464/1‑139 (MQ=255)
                       cGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTc                                                                                                               <  6:188467/139‑1 (MQ=255)
                         aTTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCaa                                                                                                             <  5:282232/139‑1 (MQ=255)
                               aataaACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCttt                                                                                                       >  6:288553/1‑139 (MQ=255)
                                ataaACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCtttt                                                                                                      <  8:578/139‑1 (MQ=255)
                                            ctcgCCATCTCCACATTGTCGAAAGAGCAGAATATCACTGTAGGAAATCTGCTATTTCACTCTAAAACCCCCAACATCACCGGAAGGATCATAGATTGTAATGGTCATACCACCTTTCACCCCTTTTCCAATGTAGcca                                                                                          <  6:100216/136‑1 (MQ=255)
                                            aTAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATATATTGTAATGGTGTTACCACCTTTCAACCCTTTTCCAATGTAGaca                                                                                          <  3:208710/139‑1 (MQ=255)
                                               gCCATCTGCAATTTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGACGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCAAATGTAGACACGc                                                                                       <  6:93571/139‑1 (MQ=255)
                                               cccATCTGCACATTGTTGAAACAGCAGAATATCACTGTATCAAATCTGGTCTTTCACTCTAAAATCCCGAACATCACGGGATGGATCATAGATTGTAATGGTTATACCACGTTTCAACCCTTTACCATTGTAGACACGc                                                                                       <  4:180069/138‑1 (MQ=255)
                                                    cTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCtg                                                                                  >  1:301226/1‑139 (MQ=255)
                                                          aTTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATc                                                                            >  3:180069/1‑139 (MQ=255)
                                                          aTTGTTGAAAGAGCACAATATCACTGTATGAAATCTGCTATTTCACTATAAAATCCAGAACATCACCGGAAGGATCATAGATTGTAATGGTTATACCACCTTTCAACCCTTTTCAAATGTAGACACGCCCCTGTGGATc                                                                            <  7:288214/139‑1 (MQ=255)
                                                            tgttgAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAg                                                                          >  3:3057/1‑139 (MQ=255)
                                                              ttgAGAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCa                                                                        >  4:331774/1‑139 (MQ=255)
                                                                 aaaaAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATg                                                                      >  5:333406/1‑138 (MQ=255)
                                                                 aaaGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTc                                                                                                               >  3:375764/1‑97 (MQ=255)
                                                                 aaaGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTc                                                                                                               <  4:375764/97‑1 (MQ=255)
                                                                 aaaGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATg                                                                      >  2:256002/1‑138 (MQ=255)
                                                                      gCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTTGACACGCCCCTGTGGATCAGCATGGaaaa                                                                 >  5:301407/1‑138 (MQ=255)
                                                                        aGAATTACACTGTATTAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGACGGATCATACATGCTAATGGTTATCCCACCTTTCAACCCTTTTCCAATGACCACACGCCCCTGTGGATCAGCATGGAAAAAAc                                                              <  7:330299/139‑1 (MQ=255)
                                                                               cACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGCACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGCCACGCCCCTGTGGATCAGCAGGGAAAAAACACTCTTg                                                       >  4:226247/1‑139 (MQ=255)
                                                                               cACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGaaaaaa                                                               >  8:222089/1‑131 (MQ=255)
                                                                               cACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGaaaaaa                                                               <  7:222089/131‑1 (MQ=255)
                                                                                 cTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTggg                                                     >  6:43519/1‑139 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATacc                                                                                                                      <  8:303809/71‑1 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATacc                                                                                                                      >  7:303809/1‑71 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCACCCCTTTTCCAATGTAGACCCGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGgtgt                                                  >  5:291124/1‑139 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGGCACGCCCCTGTGGATCAGCATGGAAAAAACAATCATGGgtgt                                                  >  1:357044/1‑139 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGgtgt                                                  >  7:311719/1‑139 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGgtgt                                                  >  8:217968/1‑139 (MQ=255)
                                                                                    tATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAAAAATCTTGGttgt                                                  >  6:137194/1‑139 (MQ=255)
                                                                                    tATGAAATATGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATacca                                                                                                                     >  5:24494/1‑72 (MQ=255)
                                                                                    tATGAAATATGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATacca                                                                                                                     <  6:24494/72‑1 (MQ=255)
                                                                                      tGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTtc                                                >  7:154086/1‑139 (MQ=255)
                                                                                            cTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATac                                                                                                                       >  2:349773/1‑62 (MQ=255)
                                                                                            cTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATac                                                                                                                       <  1:349773/62‑1 (MQ=255)
                                                                                            cTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTtg                                          >  1:379389/1‑139 (MQ=255)
                                                                                             tGCTCTTTCACTATAAAATCCCGAACATCACCGGATGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTACCAATGTACACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTAGGTGTTCTCGTtga                                         <  5:202519/139‑1 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTg                                                                                                                                   <  8:54411/48‑1 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTg                                                                                                                                   >  7:54411/1‑48 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCa                                                                                                   <  6:242510/80‑1 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCa                                                                                                   >  5:242510/1‑80 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTtgat                                        >  3:71860/1‑139 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTtgat                                        >  7:359212/1‑139 (MQ=255)
                                                                                              gCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTtgat                                        >  3:290623/1‑139 (MQ=255)
                                                                                                         ctAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCACCTGTGGATCAGCATGGAAAAAACAATCTTGTTTGTTCTCGTTGATGATTTTATTc                              >  2:92167/1‑138 (MQ=255)
                                                                                                                  cGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGgtgt                                                  <  1:46690/109‑1 (MQ=255)
                                                                                                                  cGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGgtgt                                                  >  2:46690/1‑109 (MQ=255)
                                                                                                                  cGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCaattttaatt                    >  5:189243/1‑139 (MQ=255)
                                                                                                                              gAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATTATTTTATTCAAGTTTAATTTTTTCATTGGt         >  3:221804/1‑138 (MQ=255)
                                                                                                                              gAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGtt        <  8:154086/139‑1 (MQ=255)
                                                                                                                               actgATCATAGATTTGAATTGGTGATCCCACCTTACAACCGTTTACCAATGTCCCCACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGtt        <  4:30977/136‑1 (MQ=255)
                                                                                                                               acggATCATAGATTTTTAATGGTTAAACCACATTAAACCCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGtt        <  4:36464/137‑1 (MQ=255)
                                                                                                                                  cataaTAGATTGTAATGGTGATACCACCTTAAAACCCTTTTCCAATGTACCCACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGTTATg     <  5:32470/134‑1 (MQ=255)
                                                                                                                                   aTCATAGATTGTAATGGTGATACCACCTTTCACCCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGTTATGCt   >  1:42496/1‑139 (MQ=255)
                                                                                                                                   aTCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCtt                                                        <  1:191003/86‑1 (MQ=255)
                                                                                                                                   aTCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCtt                                                        >  2:191003/1‑86 (MQ=255)
                                                                                                                                    taaTAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACATGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGTGTGTTCTCGTTGATTATTTTTTTCAATTTTAATTTTTTCATTTGTTTTGGTc  >  4:38608/3‑139 (MQ=255)
                                                                                                                                    tCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGTTATGCTc  >  1:318300/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                    
AAAACCTTTACCTTCGGCTCGTACGATTCGTAATAAACGAAGACATAGCCATCTGCACATTGTTGAAAGAGCAGAATATCACTGTATGAAATCTGCTCTTTCACTCTAAAATCCCGAACATCACCGGAAGGATCATAGATTGTAATGGTGATACCACCTTTCAACCCTTTTCCAATGTAGACACGCCCCTGTGGATCAGCATGGAAAAAACAATCTTGGGTGTTCTCGTTGATGATTTTATTCAATTTTAATTTTTTCATTGGTTATGCTC  >  BS168_reference_genome/4021047‑4021317

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: