Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,023,150 C→A 30.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,1500CA30.0% 118.9 / 8.2 76intergenic (+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 (8/39);  new base A (18/5);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 9.99e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTGT  >  BS168_reference_genome/1023026‑1023285
                                                                                                                              |                                                                                                                                       
ctgGCGTCTGGCGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGGCGGGTTTTTTTAAACGCATAAAAAAAAAAACCCCCGCAGAGCGGGGGTTTTTTGTGGTGGTTAAATTAAATAATAAAAAAtt                                                                                                                            >  8:367663/1‑139 (MQ=255)
ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATTGCGGCCGTGTTTTTTAAAAGAAAAAAAACAAAAAACAACAAAAAAGGGGGGTTTTTTTGTTGTGTAAAATTTAAAAAAAAAAAAatt                                                                                                                           >  6:373642/1‑139 (MQ=255)
ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTTAACGTAATCAAACAAAAAAAACCGCAACAGCGGGTGTTTTTTTTGGTGAAAAATTTaaaaaaaaaaaaaat                                                                                                                           >  3:270988/1‑136 (MQ=255)
ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGGTTTTTTTGAACGTAATCAACGAAAAAACACCTCAACAGCGGGGGTTTTGTGGGGTGGAAAATTTAAAAAAAAAAAAatt                                                                                                                           >  2:245348/1‑139 (MQ=255)
             aaTCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAAACTAATAAACGAAAAAACACCAAAAAAGCGGGGTTTTTTGGTGGTGAAAAATTTAAAAATAAAAAAATTTTTGATTTTTca                                                                                                               >  8:202838/1‑137 (MQ=255)
               tCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAAAGTAATAAAAGAAAAAAAACCTCAAGAGCTGGTTTTTTGTTTGGGTGTAAATTTAAAAAAAAAAAAATTTTAGATTTTTATaa                                                                                                             >  8:79811/1‑138 (MQ=255)
                         tGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTCAACGTAATCAACGAAAAAACAGCTCAAGAGCTGGTTTTTTGTTTGGTGTAAAATTTGAAAAAAAAAAATTTGTGGATTTTTGTAAAAATTTTGttt                                                                                                  >  1:265880/1‑139 (MQ=255)
                            cAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAAAGTAATAAACGAAAAACCCGCTCAAGAGCTGGGTTTTTTTGTGGTGAAAAATTTAAAAATAAAAAATTTTTCTATTTTTATAAAAATTTTGttttaa                                                                                               >  6:310468/1‑137 (MQ=255)
                            cAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTCAACGTAATCAAAGAAAAACCAGCACAAGCGCTGGTTTTTTGTGTGGTGTAAAATTTGAAAATAAAAAATTTTTAGATTTTTGTAACAATTTTGttttaa                                                                                               >  2:80374/1‑137 (MQ=255)
                               ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAACGTAATAAACGAAAAAACAGCGCAAGAGCTGGTTTTTTGTTTGGTGGAAAATTTAAAAAAAAAAAAATTTTGGATTTTTGTTAAAATttttttttgtttt                                                                                            >  1:84421/1‑134 (MQ=255)
                               ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTGAAAGTAATAAACGAAAAACCACCTCAACAGCTGGTTTTTTGTGTGGTGTAAAATTTAAAAATAAAAAATTTTTCGATTTTTTTAAAAATTTTTTTTTAAttt                                                                                            >  3:64149/1‑139 (MQ=255)
                               ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTGAAAGTAATAAACGAAAAAACAGCTCAAGAGCTGGTTTTTTTTGTGGTGTAAAAATTAAAAATAAAAAATTTTTGGATTTTTGTAACAATTTTTTTTTGAttt                                                                                            >  1:61847/1‑139 (MQ=255)
                               ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTCAACGTAATAAACGAAAAAACACCACAACAGCTGGTTTTTTTTGTGGGTAAAAATATCAAAAAAAAAAATTTTTGGATTTTTCTAACAATGTTGTTTTCAttt                                                                                            >  5:429018/1‑139 (MQ=255)
                                 cAGACCATTCATGTGAATGGCGGCCGTTTTATTTGAAACTAATAAACCAAAAAACACCGCAACAGCTGGTTTTTTGTGTTGTGATAAAATTGAAAATAAAAAAATTTTAGATTTTTATAAAAAATTTTTTTTAAttttt                                                                                          >  3:288607/1‑139 (MQ=255)
                                    aCCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATAAACGAAAAAACACCACAAGAGCGGGTTTTTTGTGTGGTGTAAAATTTAAAAATAAAAAAATTTTGGATTTTTTTTAAAATTTTTTTTTAAttttttg                                                                                        >  2:91773/1‑136 (MQ=255)
                                    aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCa                                                                                       <  3:242773/139‑1 (MQ=255)
                                    aCCATTCATGTGAATGGCGGCCGTGTTTTTTTAAAATAATAAACCAAAAAACAACTAAAGAACTTTTTTTTTTTTTGGTGAATAATTTTAAAAAAAAAAAATTTTTTATTTTTTTTACAAATATTTTTTTAttttttt                                                                                        >  4:242773/1‑136 (MQ=255)
                                    aCCATTCATGTGAATGGCGGCCGTGTTTTTTCAACGTAATAAACGAAAAACCCCCTACAGACCTGGTGTTTTTTGTGGGGGAAAATATCAAAAAAAAAAACTTTAGGAATTTTGTATCAAATTTGTTTTGAgtttttc                                                                                        >  6:206865/1‑136 (MQ=255)
                                       gacagtgtgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                    <  8:50395/134‑1 (MQ=255)
                                       aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                    <  3:137983/139‑1 (MQ=255)
                                         tCATGTGAATGGCGGCCGTTTTTTTTGAAAGGAAAAAACCAAAAAACAACAAAAGAGCTGGTTTTTTTTGTTGTTTAAAATTTGAAAATAAAAAAATTGTGGATTTTTATTAAAAATTTTTTTTAATTTTTTGAAtttt                                                                                  >  2:360447/1‑139 (MQ=255)
                                         tCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATAAAACAAAAACCAGCTAAAGAGCGGGTTTTTTGTGTTGTGTAAAATTTAAAAAAATAAAATTTTTTTGTTTTTGTTTCAATTTTTTTTTGATTTTTTCAtttt                                                                                   >  5:332641/1‑138 (MQ=255)
                                         tCATGTGAATGGCGGCCGTGTTATTTGAACGGAATAAACGAAAAACCAGCGCAAGAGCGGGGGTTTTGTGTGGTGTAAAATTTAAAAATAAAAAATTTTTCGATTTTTTTAACAATTTTTTTTTGATTTTTTGAAtttt                                                                                  >  2:246309/1‑139 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  7:62131/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  4:288607/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  7:79811/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  1:246309/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAAAGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  7:179011/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTAATTCAACGTAATCAACGAAAAACCACCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  5:42979/138‑1 (MQ=255)
                                              tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                             <  3:396039/139‑1 (MQ=255)
                                              tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAAATTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                             <  3:398862/139‑1 (MQ=255)
                                                aaTGGCGGCCGTTTTATTTGAACCTAATCAACCAAAAAACACCTCAAGAGCTGGGTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTGGATTTTTGTAACAAATTTTTTTTAATTTTTGGAATTTtttttat                                                                            >  5:417545/1‑138 (MQ=255)
                                                aaTGGCGGCCGTTTTATTTCAACGTAATAAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTTTGGTGAAAAATTTCAAAAAAAAAAATTTTTAGATTTTTATAACAATTTTTTTTTCATTTTTTGAATTTtttatat                                                                            >  5:298106/1‑138 (MQ=255)
                                                aaTGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTGGATTTTTATAACAATTTTGTTTTAATTTTTGCAATTTtatatat                                                                            >  8:326249/1‑138 (MQ=255)
                                                aaTGGCGGCCGTGTTTTTTGAACGTAATAAACCAAAAACCAGCTCAAGAGCGGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTGGGTTTTTTTAAAAATTTTTTTTTTATTTTTTGAAtttttttttt                                                                            >  8:79630/1‑132 (MQ=255)
                                                aaTGGCGGCCGTGTTTTTTAAAAGTAATAAACGAAAAAACACCTCAAGAGCGGGTTTTTTGTGTGGTGAAAAAATTAAAAAAAAAAAAATTTTGGATTTTTATAAAAATTTTTTTTTGATTTTTGAAAATTtttttat                                                                            >  8:225595/1‑138 (MQ=255)
                                                  tGGCGGCCGTTTTATTTAAAAGTAATAAACGAAAAACCAGCACAAAAGCGGGTTTTTTTTTTGGGGAAAAAATTAAAAATAAAAAAATTTTAGATTTTTATAAAAAATTTTTTTTAATTTTTTCAATTTTATATataa                                                                          >  2:343316/1‑138 (MQ=255)
                                                    gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTTTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                        <  7:18207/138‑1 (MQ=255)
                                                    gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                        <  2:153632/138‑1 (MQ=255)
                                                    gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                        <  6:277499/138‑1 (MQ=255)
                                                    gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                        <  5:298747/138‑1 (MQ=255)
                                                    gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                        <  4:64149/138‑1 (MQ=255)
                                                     cggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATAtt                                                                      <  2:13641/139‑1 (MQ=255)
                                                        gcCGTTTTAATTCAACGTAATCAACGAAAAAACAGCAGAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTGTAGATTTTTATAACAATTTTTTTTTCATTTTTGGAATTTTTTATTTAatatgta                                                                    >  7:50395/1‑136 (MQ=255)
                                                         ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                   <  1:91773/139‑1 (MQ=255)
                                                         cagTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                                   <  3:226767/137‑1 (MQ=255)
                                                                       cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                    <  5:140293/108‑1 (MQ=255)
                                                                       cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                                    >  6:140293/1‑108 (MQ=255)
                                                                                cGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGa                                            <  2:61847/139‑1 (MQ=255)
                                                                                 gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAc                                           <  5:373642/139‑1 (MQ=255)
                                                                                 gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGTAATGATTCGCGAc                                           <  7:203741/139‑1 (MQ=255)
                                                                                 gAAAAACCAGCTAAAGAGCGGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAc                                           <  7:112271/139‑1 (MQ=255)
                                                                                 caaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAc                                           <  7:225595/138‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCACGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  1:354438/139‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCAATAGATGGTTTTTTGTTTGGTGAAAAATTTAAAAATAAAAAATTTGTCTATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGATAAAGGAATGATTCGCGAcc                                          <  1:343316/139‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  7:202838/139‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGTTTAAAAATTTCAAAATAAAAAATTTGTCTATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  6:429018/139‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  7:79630/139‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  7:367663/139‑1 (MQ=255)
                                                                                  aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  3:36285/139‑1 (MQ=255)
                                                                                  aaaaaCAAGCTCAAGTGTTGGTTTTTTGTGTGGTTAAAAATTTTAAAATAAAATATTTGTCTATTTTTCTAAAAATTTTGTTTTAATTTTTGCAATTTTATATATAATTTTGATTAGGCCAAAGGAATGATTCGCGAcc                                          <  4:393030/139‑1 (MQ=255)
                                                                                   aaaaCCAGCTCAAGAGCTGTTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTTTTTTCATTTTTGCAATTTTATATAAAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  3:377574/138‑1 (MQ=255)
                                                                                   aaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  4:126450/138‑1 (MQ=255)
                                                                                     acccAGCCCACTAGCTGGTTTTTTGGGGGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTCATTTTTGCAATTTTATATATAAAATTGAGTTGGACAACGGAAGGATTCGCGACCAGt                                       <  4:322955/137‑1 (MQ=255)
                                                                                        cAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  >  6:328472/1‑93 (MQ=255)
                                                                                        cAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                                  <  5:328472/93‑1 (MQ=255)
                                                                                           cTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          >  2:290965/1‑130 (MQ=255)
                                                                                           cTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                          <  1:290965/130‑1 (MQ=255)
                                                                                                       ttttttGTGTGTTTAAAATTTTCAAAAAAAAAATTTTGTCGATTTTTATAAAAATTTTGTTTTAATTTTTGCAATTTTATATATAATTTTGATTAGGAGAAAGGAAGTATTCGCGACCAGTACAAATATAGTAGaaaaa                     <  6:128502/139‑1 (MQ=255)
                                                                                                           ttGTTTGGTTAAAATTTTAAAAATAAAATATTTGTCGATTTTTATAACATTTTTTTTTTTTTTTTTGATATTTTTTATATATTATTTAGTTGGAGAAAGTATTTTTTCGCGACCAGTACAAATATAGTAGaaaaaaaa                  <  7:370740/138‑1 (MQ=255)
                                                                                                           ttGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCTATTTTTCTAAAAATTTTTTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATTATTCGCGACCAGTACAAATATAGTAGaaaaaaaa                  <  1:360447/138‑1 (MQ=255)
                                                                                                           ttGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGaaaaaa                    >  6:95085/1‑136 (MQ=255)
                                                                                                              tgtgGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATTATGCGCGACCAGTACAAATATAGTAGAAAAAAAAGaat              <  2:265880/139‑1 (MQ=255)
                                                                                                                tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATTTATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTTGAAAAAAAAGaataa            >  8:41159/1‑139 (MQ=255)
                                                                                                                   gttaaaaaTTTAAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTCATTTTTTCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATg          <  3:398676/136‑1 (MQ=255)
                                                                                                                        aaaTTTAAAAATAAAAAATTTGTCTATTTTTCTAAAAATTTTGTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCAt     <  5:310468/139‑1 (MQ=255)
                                                                                                                           tttCAAAATAAAAAATTTGTCTATTTTTATAACAATTTTGTTTTCATTTTTTCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATGCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTGt  <  1:245348/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                       
CTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGACCATTGT  >  BS168_reference_genome/1023026‑1023285

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: