Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,023,150 C→A 34.0% intergenic (+50/‑200) yhdF → / → bcaP putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,023,1500CA34.0% 98.8 / 14.2 63intergenic (+50/‑200)yhdF/bcaPputative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter
Reads supporting (aligned to +/- strand):  ref base C (11/29);  new base A (19/3);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.13e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTG  >  BS168_reference_genome/1023024‑1023284
                                                                                                                                |                                                                                                                                        
tgctgGCGTCTGACGCATCTTCCTATATGACAGGGCAGACCATTAATGTGTATGGCGGCCGTGTTTTTTTAACGTTAAAAACGAAAAAACAACTCAAAAGCTTGGTTTTTTTTTTTGGTAAAAATTaaaaaaaaaaaat                                                                                                                               >  8:239888/1‑138 (MQ=255)
tgctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAAAGTAAACAAAGAAAAACCCCCTCAAGACCTGGGGTTTTTTGTGGGGTAAAATATAAAAATaaaaaa                                                                                                                               >  5:147283/1‑138 (MQ=255)
tgctgGCGTCTGACGAATCTTCCGATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTGAACGTAATCAAACAAAAACCCCCTCAACACCGGGGTTTTTTTGTGGTGTAAAAATtaaaaataaaaatt                                                                                                                              >  7:124798/1‑137 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACCGAATCAACCAAAAACCCCCACAAGCGCGGGGTTTTTTTGTGGTGTAAAATTTAAAAAAAAAAAAatt                                                                                                                            >  5:197509/1‑139 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACCAAAAAACAACGCAAGAGCAGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAAatt                                                                                                                            >  5:178835/1‑139 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAACGTAATAAAACAAAAAAAACCACAAAAGCTGGTTTTTTTTGTTGTTGAAAAATTTAAAAAAAAAAAAtt                                                                                                                             >  2:53830/1‑139 (MQ=255)
                 tCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAAAGTAATAAACGAAAAACAAACTCAAGACCTGGTTTTTTTTGTGGTGAAAAAATTAAAAAAAAAAAATTTTTAGATTTTTTTat                                                                                                             >  7:11481/1‑137 (MQ=255)
                      cTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTAGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAAtt                                                                                                       <  3:193511/138‑1 (MQ=255)
                        atatGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAAAGTAATCAACCAAAAAAAACCTCAAGAGCGGGTTTTTTTTGTGGTGTAAAATTTCAAAATAAAAAATTTTTAGATTTTTGTTACAAATTTg                                                                                                    >  5:303930/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTCAACGTAATCAAAGAAAAACCAGCTCAAGAGCTGGTTTTTTTTGTGGTGAAAAAATTAAAAATAAAAAATTTTTAGATTTTTATAACAATTTTTTTTTAAttt                                                                                           >  5:82727/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTCAACCGAAACAACCAAAAAAAAACACAAGAGAAGGTGTTTTGTGTGTTGTAAAATTTAAAAAATAAAAATATTTAGAATTTTTTAACAATTTTTGTTTTAAtt                                                                                            >  8:347488/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAAAGTAATAAAAGAAAAAACACCACAAGAGCTGGTTTTTTTTTTGGTGAAAAATTTAAAAAAAAAAAAATTTTAGATTTTTTTTACAATTTTTTTTTTAttt                                                                                           >  6:280660/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTGAACGTAATAAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTGTCGATTTTTGTATAAATTTTTTTTTGAttt                                                                                           >  7:79321/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTGAACGTAATAAACGAAAAAACAGCTCAAGAGCGGGTTTTTTTTGTGGTGTAAAATTTAAAAATAAAAAATTTTTAGATTTTTATAACAATTTTTTTTTTAttt                                                                                           >  7:369091/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATAAACGAAAAAACACCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTGTAGATTTTTATAACAAATTTTTTTTAAttt                                                                                           >  8:417683/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACCTAATCAACCAAAAACCACCTCAAGAGCTGGTTTTTTTTGTGGTGAAAAAATTAAAAATAAAAAAATTTTAGATTTTTTTAACAAATTTTTTTTAAttt                                                                                           >  6:329003/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTGAACGTTATCAACGAAAAACCACCGCAAGAGCTGGTTTTTTGTTTGGTGTAAAATTTAAAAATAAAAAAATTGTGGATTTTTATTAAAAATTTTTTTTTAttt                                                                                           >  7:181616/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTGAACGGAAATAACCAAAAAACAACAGAACAGCGGGTTTTTTTTTTGGTGGAAAAATTAAAAATAAAAAAATTTTTGATTTTTGTTATAATtttttttttgttt                                                                                           >  4:269567/1‑134 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAACGTAATAAAAGAAAAAACAGCTCAAGGGCGGGGTTTTTGTGTGGGGTAAAATTTAAAAAAAAAAAATTTGTGGATTTTTATAAAAAATTTTTTTTTAttt                                                                                           >  5:373758/1‑139 (MQ=255)
                                   cAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCAttttt                                                                                         >  7:417436/1‑139 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                        <  3:1194/139‑1 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCACCTCAAGAGCTGGGTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTCGATTTTTATAAAAATTTTGTTTTAAttttttg                                                                                       >  1:1346/1‑136 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTGTTATTTAAAAGTAATAAAAGAAAAAACAGCGCAAGAGCTGGGTTTTTTTGTTGTGTAAAATTTAAAAAAAAAAAAATTTTGGATTTTTTTTAAAAATTTTTTTTTAttttttg                                                                                       >  4:416361/1‑136 (MQ=255)
                                      aCCATTCATGAGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                      <  5:276252/139‑1 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                   <  4:82927/139‑1 (MQ=255)
                                           tCATGTGAATGGCGGCCGTTTTATTTAAACCTAATAAACCAAAAACCCCCACAAGCGCGGGTGTTTTTGGTGGGGAAAAATTTCAAAAAAAAAAAATTTTGGCTTTTTATAACAAATTTTTTTTTATTTTTGCAAtttt                                                                                 >  6:362581/1‑139 (MQ=255)
                                           tCATGTGAATGGCGGCCGTGTTTTTTTAAAGTTATAAAAGAAAAAACCACACAAGAGCGGGGTTTTTTGGTGGTTTTAAAATTTAAAATAAAAAAATTGTGGTATTTTTTTAAAATTTTTTTTTTATTTTTTGAAtttt                                                                                 >  7:164463/1‑139 (MQ=255)
                                            cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                 <  2:227635/138‑1 (MQ=255)
                                            cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                 <  5:280660/138‑1 (MQ=255)
                                            cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGGTTTCATTTTTGCAAtttt                                                                                 <  5:380500/138‑1 (MQ=255)
                                             aTGTGACTGGCGGACGTTTTATTTCAACGTAATCAACGAAAAACCCGCTCAAGAGCTGGTTTTTTGTGTTGTGAAAAATTTAAAAATAAAAAATTTGTCGATGTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                               <  7:347488/139‑1 (MQ=255)
                                             aTGTGAATGGCGGCCGTTTTATTTAAACTTAATAAACTAAAAACATTCTCACTAGATGGTTTTTTTTTTGGTTAAAAATTGAAAAATAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                               <  8:124798/138‑1 (MQ=255)
                                                tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  1:354102/139‑1 (MQ=255)
                                                tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  8:181616/139‑1 (MQ=255)
                                                tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  7:239888/139‑1 (MQ=255)
                                                tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                            <  3:416361/139‑1 (MQ=255)
                                                  aaTGGCGGCCGTTTTATTTGAAAGTAATCAACGAAAAAACACCTCAAGAGCGGGGTTTTTGTGTGGGGAAAAAATTAAAAAAAAAAAATTTTTCGATTTTTTTTAAAATTTTTTTTTAATTTTTTTAATTTtatatat                                                                           >  3:82927/1‑138 (MQ=255)
                                                  aaTGGCGGCCGTGTTATTTGAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGGTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTTACAATTTTGTTTTCATTTTTGCAATTTtattttt                                                                           >  7:116236/1‑134 (MQ=255)
                                                    tGGCGGCCGTTTTATTTAAACGTAATAAACCAAAAAACAACGCAATAGCTGGGTTTTTTTGTGGTGTAAAATTTAAAAATAAAAAATTTTTTGATTTTTATAAAAAATTTTTTTTATTTTTTTGAATTTtttaaatta                                                                         >  6:102443/1‑136 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  7:96618/138‑1 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  7:417683/138‑1 (MQ=255)
                                                      gcggcCGTTTTATTTCAAAGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAAATTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATCTataata                                                                       <  6:186630/138‑1 (MQ=255)
                                                       cggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATAtt                                                                     <  5:285473/139‑1 (MQ=255)
                                                          ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                  <  3:272134/139‑1 (MQ=255)
                                                             ttttATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAg                                                               <  8:417436/139‑1 (MQ=255)
                                                                        cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAAt                                                    >  1:155231/1‑139 (MQ=255)
                                                                            aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       >  1:143885/1‑116 (MQ=255)
                                                                            aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                       <  2:143885/116‑1 (MQ=255)
                                                                                accaaaaaCCAGCTCAAGAGCTGTTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCAAACAATTTTGTTTTCATTTTTGCAATTTTAAATATAATATTGAGTAGGAGAAAGTAATGATTcgcg                                            <  6:373758/136‑1 (MQ=255)
                                                                                 cGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGATTGATTCGCGa                                           >  8:396507/1‑139 (MQ=255)
                                                                                  gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTATAACAATTTTTTTTTCATTTTTGCATTTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAc                                          <  4:98765/139‑1 (MQ=255)
                                                                                  gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAc                                          <  5:92757/139‑1 (MQ=255)
                                                                                   aaaaaCCCGCGCAAGCGCTGGTTTTTTGTGTGGTTAAAAATTACAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAACGGAATGATTCGCGAcc                                         <  5:362581/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTTGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  3:269567/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  6:303930/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTAATAAATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  3:81912/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGATAAAGGAATGATTCGCGAcc                                         <  8:369091/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                         <  1:53830/139‑1 (MQ=255)
                                                                                    aaaaCCAGATCAATAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTTGCGTTTTTTCAAACATTTTTTTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAAGTATTCGCGAcc                                         <  6:147283/138‑1 (MQ=255)
                                                                                     aaaCCAGCTAAAGAGGTGGTTTTTTGTGTGTTTAAAAATTTAAAAATAAAAAATTTGTCTATTTTTCTAAAATTTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGATAAAGGAATTATTCGCGACCAg                                       <  8:11481/139‑1 (MQ=255)
                                                                                           gCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACaaa                                 >  1:38734/1‑139 (MQ=255)
                                                                                                       ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGaa                       >  5:111324/1‑138 (MQ=255)
                                                                                                                   tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGCCCAGTACAAATATAGTAGAAAAAAAAGaataa           >  5:422027/1‑139 (MQ=255)
                                                                                                                   tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataa           >  4:301202/1‑139 (MQ=255)
                                                                                                                            aTTTCAAAATAAAAAATTTGTCTATTTTTCTAACAATTTTGTTTTAATTTTTCCAATTTTATATATAATATTGAGTAGGAGAACGTTTTGATCCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTg  <  5:102443/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                        
TGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTG  >  BS168_reference_genome/1023024‑1023284

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