Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,087,669 C→T 6.3% H72Y (CAC→TAC)  yxcD → conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,087,6690CT 6.3% 192.5 / 2.1 64H72Y (CAC→TAC) yxcDconserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base C (36/22);  new base T (2/2);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CTGAAGACGTTGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTT  >  BS168_reference_genome/4087532‑4087804
                                                                                                                                         |                                                                                                                                        
cTGAAGACGTTGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACca                                                                                                                                         >  4:15848/1‑139 (MQ=255)
cTGAAGACGTTGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACca                                                                                                                                         <  4:201350/139‑1 (MQ=255)
   aaGACGTTGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACaa                                                                                                                                      >  4:199059/1‑139 (MQ=255)
   aaGACGTTGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACaa                                                                                                                                      >  4:204778/1‑139 (MQ=255)
        gTTGAAGTAAAAGAACTATCGTACAACAAGAGTACAGGGTTATTCGATGAAGCTACATTTTGCCTGCACCACGAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTATAGGAATACAACAAATTc                                                                                                                                  >  7:69166/1‑138 (MQ=255)
         ttGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACca                                                                                                                                         >  2:9047/1‑130 (MQ=255)
         ttGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACca                                                                                                                                         <  1:9047/130‑1 (MQ=255)
                       cTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACcg                                                                                                                  >  1:12878/1‑139 (MQ=255)
                       cTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACcg                                                                                                                  >  6:204519/1‑139 (MQ=255)
                                    aGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGCCATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACTACAACTTCAATCCGGACGTTACCGTCGTGGAAATgg                                                                                                      >  6:31909/1‑137 (MQ=255)
                                    aGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACTACAACTTCAATCAGGACGTTACCGTCGTTGAAtttc                                                                                                      >  6:35904/1‑136 (MQ=255)
                                    aGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGc                                                                                                      >  3:262517/1‑138 (MQ=255)
                                          cAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCCCCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGACTTGCAATTTGa                                                                                               >  4:262402/1‑139 (MQ=255)
                                            gggTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGAt                                                                                              >  2:165874/1‑138 (MQ=255)
                                            gggTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATa                                                                                             <  2:251431/139‑1 (MQ=255)
                                               ttcttcGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTTAATAGCAATATTATaat                                                                                           >  2:40402/1‑137 (MQ=255)
                                                  ttcGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGaaa                                                                                        >  7:30887/1‑138 (MQ=255)
                                                  ttcGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGaaa                                                                                        >  1:107039/1‑138 (MQ=255)
                                                       tGAAGCTACATTTGGCCTGCCCCACAAGCAGCTTATGTCAGATGACCTTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTATCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGc                                                                                  >  4:134318/1‑139 (MQ=255)
                                                       tGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGc                                                                                  >  8:48849/1‑139 (MQ=255)
                                                       tGAAGCTACATTTGGCCTGAACCACAAGCAGCTTATGTCAGATGAAATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGTACGTAACCGGCGTGGAATTTCAATTTTATAAGAAAAAAgcc                                                                                  >  2:54388/1‑137 (MQ=255)
                                                         aaGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCtt                                                                                >  6:61734/1‑139 (MQ=255)
                                                         aaGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCtt                                                                                >  1:160786/1‑139 (MQ=255)
                                                                  ttGGCCTGCACCACAAGCAGCTTATGTCATATGACATTTCAGAAGGCATTATTCAGTTCTTAGTGTAATACCACAAGTTCATTCCGGAGTTTACCGTCTTTTATTTGCATTTTGATAAGAAAAAAGGCTTCTCCGCCt                                                                       <  1:40402/138‑1 (MQ=255)
                                                                     gCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGTATTGCAATTTTATTAGAAAAAAGGCATCTCccccggggg                                                                    >  5:43538/1‑134 (MQ=255)
                                                                         gCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCATAGAGCAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTTATAAGAAAAAAGGCATATCCGCCTTAGtattt                                                               >  8:13414/1‑139 (MQ=255)
                                                                          caccacAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGGAATTTGATAATAAAAAAGGCTTCTCCGCCTTGGTTTTTg                                                              >  3:35003/1‑139 (MQ=255)
                                                                              acaAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCTTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCaa                                                          <  6:207750/139‑1 (MQ=255)
                                                                              acaAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTTAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCaa                                                          <  1:165874/139‑1 (MQ=255)
                                                                              acaAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCaa                                                          <  3:199059/139‑1 (MQ=255)
                                                                                    agcTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAAGCTTCTCCGCCTTGGTTTTTGTCAACGAAg                                                     >  6:129506/1‑138 (MQ=255)
                                                                                      cTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCag                                                  <  5:204519/139‑1 (MQ=255)
                                                                                          tGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTCCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCagaaga                                              <  8:30887/139‑1 (MQ=255)
                                                                                                attaCATTTCAGAAGTCATTATTCAGTTCTTACAGTAATACTACAAATTCAATCCGGACGTTACCGTCGTTTAATTGCAATTTTAAACCAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATaaaa                                        <  5:31909/136‑1 (MQ=255)
                                                                                                aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGTATTTGAATTTTATTATAAAAAAGGCTTGTCCGCCGTTGGGTTTTGGAAAGAAGGCGAAGAATaaa                                         >  1:241993/1‑138 (MQ=255)
                                                                                                aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGCATTTCAAATTGATAAGAAAAAAGGCTTCTCCGCCTTGGGTTTTTTCAACCAAGCAGCAGAATaaa                                         >  1:168475/1‑138 (MQ=255)
                                                                                                aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTGCCCCTTGGTTTTTGTCAACGAAGCAGAAGAATaaa                                         >  2:257268/1‑138 (MQ=255)
                                                                                                aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTGGTATTGCAAATTGTTAAGAAAAAAGGCTTCTCCGCCGTGGTTTTTGTCAACGCAGCAGAAGAATaaa                                         >  3:206865/1‑138 (MQ=255)
                                                                                                aTGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCCTCGTGGAATTGCAATTTTCAAAGTAAAAAGGGTTCGCCGCCATGGGTTTTGTCAAAGAAGGAGCAGAATtaa                                         >  1:65734/1‑138 (MQ=255)
                                                                                                 ttaCATTTAATATGGCATTATTCAGTTCTTCGAGGAATAATACAACTTCAATCCGGACGTTACCGTCGTTGAATTCCATTTTGATAACAAAAACGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATaaaa                                        <  5:35904/136‑1 (MQ=255)
                                                                                                    cATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTTCAAATTGATAAGAAAAAAGGCGTGTCCGCCATGGTTTTTTTCTAAGAAACAGATAAATAAAAAtt                                     >  2:251118/1‑137 (MQ=255)
                                                                                                           taaGCCATTAATCAGTTCTATCCGGAATACCACAACATCACACCGCACGTTACCGCGGTTTAATTCCAATTTCTAAATAAAAAAGGCTTCTCCGCCTTGGATTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  4:262517/138‑1 (MQ=255)
                                                                                                           gAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAACGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  7:11301/139‑1 (MQ=255)
                                                                                                           gAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAACGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  2:160786/139‑1 (MQ=255)
                                                                                                           gAAGGCATTATTAAGTTCTTAGAGGAATCCCACAACTTCAATCCGGCCGTAACCGTCGTTGAATTGCAATTTGATAACAAAAAACGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  2:107039/139‑1 (MQ=255)
                                                                                                           gAAGCCATTAATCAGTTATTACACTAAAACCACAACTTCAATCCGCACGTTACCGGCGTTTAATTGAAATTTAATAAAAAAAACCGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  3:15848/139‑1 (MQ=255)
                                                                                                           catgtcATTATTCTTTTATTAGTGAAACACCACAACTTAAACCGGGCCGTTCCGGTTGTTTATTTGAATTTTAAAAACAAAAACCGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  6:43538/134‑1 (MQ=255)
                                                                                                            aaGGCATTATTCAGTTCTAAGAGCAATTCAACAACTTCAACACGGACGTAACCGTCGTTTAATTCAAAATTCATAAGAAAACCCGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCa                            <  7:104742/139‑1 (MQ=255)
                                                                                                                 attattCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCGCCGCCGTGGTTTTTGTAAACGAAGCAGAACAATAAAAATGTCTTAATCAAACCt                       >  2:186047/1‑139 (MQ=255)
                                                                                                                 attattCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAAATTGATAAAAAAAAAGGCTTCGCCGCCATGGTTTTTGTCAACGAAGCACAATAAAAAAAAAGTATTAATCaaaac                        >  4:79033/1‑136 (MQ=255)
                                                                                                                    attCAGTTCTTTGAGGAATACCACAAGTACAAGGCGGCCTTGACCGTCTTTTAATGGATATTTGATAAGAATAAAGGCTTCTCCGCCTTGGTTTTTGTCCACGAAGCAGAAGAATAAAAATGTCTTAATCAAACATTa                     <  8:69166/138‑1 (MQ=255)
                                                                                                                    attCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTa                     <  2:12878/138‑1 (MQ=255)
                                                                                                                        aGTTCTTAGAGTAATACCACAACTTCAATCCGCACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAACGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcg                <  3:109245/139‑1 (MQ=255)
                                                                                                                        aGTTATTAGAGGAATACCACAACTTCAATCCGGACGTTACCGGCTTTTAATTGCAATTTGAAAAGAAAAACGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcg                <  4:35003/139‑1 (MQ=255)
                                                                                                                         gTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAATGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             >  4:41629/1‑125 (MQ=255)
                                                                                                                         gTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATc                             <  3:41629/125‑1 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGGATTTCAATTTGATAAGAAAAAACGCTTATCCGCCCTGGGTTTTGTCAACGAAGCACAAGAATAAAAATTTCTTAAAAAAACCTTACTCcccg              >  1:128947/1‑139 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTTCAATTTGATAAGAAAAAAGGCTTCTCCCCCTTGGTTTTTTTGAACGAAGCAGAAGAATAAAAATGTATTAATAAAACCTTACTAcgcg              >  8:124283/1‑139 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGaaaa                                                                                       >  4:279108/1‑67 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGaaaa                                                                                       <  3:279108/67‑1 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcgcg              >  5:109001/1‑139 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcgcg              >  2:41367/1‑139 (MQ=255)
                                                                                                                          ttCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCcgcg              >  8:198790/1‑139 (MQ=255)
                                                                                                                                 agGATTACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATACGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTaa       <  3:79033/139‑1 (MQ=255)
                                                                                                                                    agtACCACAACTTCAATCCGGAGGTTACGGTCGTTCAATTCCAATTTGATAAGAAAAACGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGt    <  3:229813/137‑1 (MQ=255)
                                                                                                                                    aaTACCACACCTACAATCCCGACGTTACCGTCGTTTAATTGCAATTTGATAACAAAAACGGCTTCTCCGCCTTGGTTTTTGTCCACGACGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGt    <  5:129506/139‑1 (MQ=255)
                                                                                                                                    aaTACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTTATAAGAAAAAACGCATCTCCGCCTTGGTTTTTGTCAACGAAGCATAATAATAAAAATGTCTTAATAAAACCAAACACCGCGCGGGTAATGt    >  2:27994/1‑139 (MQ=255)
                                                                                                                                      tACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGttt  >  1:91558/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CTGAAGACGTTGAAGTAAAAGAACTATCGTACAACAAGCGTACAGGGTTCTTCGCTGAAGCTACATTTGGCCTGCACCACAAGCAGCTTATGTCAGATGACATTTCAGAAGGCATTATTCAGTTCTTAGAGGAATACCACAACTTCAATCCGGACGTTACCGTCGTTGAATTGCAATTTGATAAGAAAAAAGGCTTCTCCGCCTTGGTTTTTGTCAACGAAGCAGAAGAATAAAAATGTCTTAATCAAACCTTACTCCGCGCGGGTAAGGTTT  >  BS168_reference_genome/4087532‑4087804

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: