Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,023,150 C→A 28.1% 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,1500CA28.1% 117.9 / 9.8 68intergenic (+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 (14/32);  new base A (15/4);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 7.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTATGTTCTGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAC  >  BS168_reference_genome/1023015‑1023279
                                                                                                                                         |                                                                                                                                  
ccTATGTTATGCTGGCGTCTGACGAATCTTCCTATATGACACGGCCGACCATTAATGTGAATGGCGGCCGTTTATATTAAAAATTAAAAAACAAAAAACAACACAAAAAAATGTTTTTTTTGTTGTGTAAAATTTaaa                                                                                                                                  >  6:154636/1‑138 (MQ=255)
         tgctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTGAAACTAATAAACGAAAAAACAGCAACAGACCGGGGTTTTTTTGTGGTGGAAAAATTTAAAATAAAAAAt                                                                                                                       >  7:5692/1‑139 (MQ=255)
         tgctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAACGTTATAAACCAAAAAACACCACAAGAGCTGGTTTTTTTTTTGGTGTAAAATTTAAAAATaaaaaaa                                                                                                                        >  3:93031/1‑138 (MQ=255)
           ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTTAAAGTTATAAAAGAAAAAACAGCGCAACAGCTGGGTTTTTGTGTGGTGTAAAATTTAAAAATAAAAAAttt                                                                                                                     >  7:194475/1‑139 (MQ=255)
           ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAAAGTAATAAAAGAAAAACAAACACAATAGCTGGTTTTTTTTTTGGTTGAAAAATTAAAAAAAAAAAAAtt                                                                                                                      >  7:124859/1‑139 (MQ=255)
           ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTGAACGTAATCAACGAAAAAACAGCGCAAGAGCGGGTTTGTTTTGTGGTGTAAAAATTAAAAATTAAAAAttt                                                                                                                     >  4:113815/1‑139 (MQ=255)
                gTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACCTAATAAACCAAAAAACACCACAAGAGCGGGGTTTTTTTTTGGTGAAAAATTTAAAAAAAAAAAAattttgga                                                                                                                >  7:94288/1‑134 (MQ=255)
                       gAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCt                                                                                                         <  1:201669/139‑1 (MQ=255)
                        aaTCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAAAAAAAAATTTTTCGAtttttgt                                                                                                         >  8:213827/1‑136 (MQ=255)
                         aTCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTaa                                                                                                       <  7:213827/139‑1 (MQ=255)
                                 atatGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAAAGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTTTCGATTTTTGTAACAAttttt                                                                                               >  7:272714/1‑138 (MQ=255)
                                          ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATAAACGAAAAAACAGCTCAAGAGCTGGGTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTGAttt                                                                                      >  7:115968/1‑139 (MQ=255)
                                          ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTAGATTTTTGTAACAATTTTGTTTTCAttt                                                                                      >  6:308441/1‑139 (MQ=255)
                                          ggCAGACCATTCATGTGAATGGCGGCCGTTGTTTTTTTTAGGGTTAAACGAAAAAAAAACGCAAGATCTGGGTGTTTGTGTGGGGTAATATTTTAAAATAAAAAATTTTTTCGTTTTTTTTTGAAATtttttttttttt                                                                                       >  2:338210/1‑135 (MQ=255)
                                          ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAAACTAAAAAACGAAAAAACAACTCAAAATCTGGGTTTTTTTGTGGTGAAAAATTTAAAAATAAAAAAATTTTAGATTTTTATTAAAAATTTTTTTTAAttt                                                                                      >  4:80958/1‑139 (MQ=255)
                                            cAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAAAGTAATAAACCAAAAAACCGCTCAAAAGCTGGTTTTTTTTGTTTTGAAAAATTTAAAAAAAAAAAAATTTTCGATTTTTTGAAAAAATTTTTTTTAAatttt                                                                                    >  2:280528/1‑139 (MQ=255)
                                             aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                   <  5:321186/139‑1 (MQ=255)
                                             aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                   <  8:94288/139‑1 (MQ=255)
                                               aCCATTCATGTGAATGGCGGCCGTTTTTTTTGAAAGTAATAAACGAAAAACCAGCTCAAGAGCTGGGTTTTTTTGTGGTGTAAAATTTAAAAATAAAAAATTTTTCGATTTTTATAACAATTTTTTTTTGAttttttc                                                                                  >  5:299224/1‑136 (MQ=255)
                                               aCCATTCATGTGAATGGCGGCCGTTTTATTTCAACGGAATAAACGAAAAAACAGCGCAAGGGCTGGTTTTTTGTGTGGGGTTAAATTTGAAAATAAAAAATTTGTCGATTTTTGTTACAATTTTGTTTTGATTTTTgg                                                                                  >  1:56561/1‑137 (MQ=255)
                                               aCCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACGAAAAAACAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAAATTAAAAATAAAAAATTTTTCGATTTTTCTAACAATTTTGTTTTAAttttttg                                                                                  >  3:351313/1‑136 (MQ=255)
                                               aCCATTCATGTGAATGGCGGCCGTTTTATTTAAACGAAAAAAACGAAAAACCACCTAAAGAGCAGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAAATTGTCGCTATTTCTAACAATTTTGTTTTCATTTTTga                                                                                  >  1:4553/1‑137 (MQ=255)
                                                  aTTCCTGTGAATGGCGGCAGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGTTGAAAAATTTCAAAATAAAAAATTTGTCTATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  3:359142/139‑1 (MQ=255)
                                                  aTTCATGTGAATGGCAGCCGTTTTATTTCAAAGTAATCAACGAAAAACCATCTCAAGATCTTGTTTTTTGTGTTGTGAAAAATTTAAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  4:210455/139‑1 (MQ=255)
                                                    tCATGTGAATGGCGGCCGTTTTTTTTGAACGTAATAAACGAAAAAACAGCAAAAGAGCTGGTTTTTTTTTTTGTGAAAAAATTAAAAAAAAAAAATTTTTCGATTTTTATAAAAATTTTTTTTTGATTTTTTGAAtttt                                                                            >  8:75541/1‑139 (MQ=255)
                                                     cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                            <  6:201181/138‑1 (MQ=255)
                                                     cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                            <  4:144566/138‑1 (MQ=255)
                                                     cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                            <  3:296336/138‑1 (MQ=255)
                                                     cATGTGAATGGCGGCAGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                            <  4:351313/138‑1 (MQ=255)
                                                      aTGTGACTGGCGGCCGGTTTAATTCAACCTAATCAACGAAAACCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  6:329514/139‑1 (MQ=255)
                                                      aTGTGAATGGCGGCCGTTTTTTTTGAAACTAATCAACGAAAAACCCGCTAAAGAGCAGGTTTTTTTTGTGGTGAAAAATTTAAAAAAAAAAAAATTTTAGATTTTTATAAAAATTTTTTTTTAATTTTTTGAATTTtat                                                                          >  6:264083/1‑139 (MQ=255)
                                                      aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  3:273656/139‑1 (MQ=255)
                                                      aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  3:46079/139‑1 (MQ=255)
                                                         tgAATGGCGGCCGTTTTATTTCAACGTAATCAAAGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAATATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                       <  1:60158/139‑1 (MQ=255)
                                                         tgAATGGCGGCAGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatat                                                                        <  1:338210/138‑1 (MQ=255)
                                                         tgAATGGCGGACGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                       <  8:124859/139‑1 (MQ=255)
                                                           aaTGGCGGCCGTGTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatatat                                                                      >  5:134940/1‑138 (MQ=255)
                                                             tGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAAACAGCTCAAGGGCTGGTTTTTTTTGTGGTGAAAAATTTCAAAATAAAAAAATTTTCGATTTTTATAACAATTTTTTTTTCATTTTTTAAATTTTATATataa                                                                    >  5:266222/1‑138 (MQ=255)
                                                              ggcggcCGTGTTATTTGAACGTAATCAAAGAAAAAAACGCTCAAGAGCTGGGTTTTTGTGTGGTGATAAATTTAAAAATAAAAAATTTGTCGATTTTTGTAAAAATTTTTTTTTGAATTTTTGAATTTtatattataaa                                                                  >  4:147225/1‑133 (MQ=255)
                                                               gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                  <  2:27023/138‑1 (MQ=255)
                                                                   ccGTTTTATTTCTACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAATAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                             <  3:327629/139‑1 (MQ=255)
                                                                   ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                             <  7:128736/139‑1 (MQ=255)
                                                                   ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                             <  3:26936/139‑1 (MQ=255)
                                                                   ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                             <  7:75541/139‑1 (MQ=255)
                                                                          tATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGga                                                        <  6:134940/138‑1 (MQ=255)
                                                                                  cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGTAAt                                               >  5:52993/1‑139 (MQ=255)
                                                                                      aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAtt                                           >  3:222892/1‑139 (MQ=255)
                                                                                          aCGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTcgcg                                       >  4:312075/1‑139 (MQ=255)
                                                                                            gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTAGCGAc                                     >  8:34883/1‑139 (MQ=255)
                                                                                            gAAAAACCACCTCAAGAGCTGGTTTTTTGAGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATATTGTTTTCATTTTTGCAATTTTATATATAAAATTGAGTAGGAGAAAGGAATGATTCGCGa                                      >  8:34876/1‑138 (MQ=255)
                                                                                             aaaaaCCTGCTCAACAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAATTTTGTCGATTTTTCTAACATTTTTGTTTTCTTTTTTGAAATTTTTTATATAATATTGAGTAGGTCAAAGGAATGATACGCGAcc                                    <  3:257553/139‑1 (MQ=255)
                                                                                             aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                    <  5:308441/139‑1 (MQ=255)
                                                                                              aaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                    <  4:93031/138‑1 (MQ=255)
                                                                                              aaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCa                                   <  5:264083/139‑1 (MQ=255)
                                                                                              aaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                    <  3:80958/138‑1 (MQ=255)
                                                                                                     gCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGAGTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGTATGAATCGCTACCAGTACaaa                            >  4:281262/1‑139 (MQ=255)
                                                                                                         aaGTGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAAtata                        <  8:115968/139‑1 (MQ=255)
                                                                                                                ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCTACCAGTACAAATATAGTAGaa                  >  7:217287/1‑138 (MQ=255)
                                                                                                                 gTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          >  5:155471/1‑81 (MQ=255)
                                                                                                                 gTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  6:155471/81‑1 (MQ=255)
                                                                                                                           tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataa      >  3:64271/1‑139 (MQ=255)
                                                                                                                           tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataa      >  1:120310/1‑139 (MQ=255)
                                                                                                                           tggtgAAAAATTTCAAAAAAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGCATAAAGTAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata       <  8:71543/138‑1 (MQ=255)
                                                                                                                              ttttaaaaTTAAAAAATAAAATATTTTTTTATTTTTATAATATTTTTTTTTTATTTTTTGAAATTTTATTTATAATATTTAGTAGAATAACGTATTTTTTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataat     <  4:64271/134‑1 (MQ=255)
                                                                                                                              gttaaaaaTTTCAAAATAAAAATTTTGTCTATTTTTCTAACAATTTTTTTTTTATTTTTGAAATTTTTAATATAATTTTGAGTTGGTGAAAGGAATGTTTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATg    <  2:4553/136‑1 (MQ=255)
                                                                                                                                 aaaaaTTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc  <  3:147225/138‑1 (MQ=255)
                                                                                                                                 aaaaaTTTCAAAATAAAAAATTTGTCGATTTTTCTAAAAATTTTGTTTTCATTTTTCCAATTTAATATATAATATTCAGTAGCAGAAAGTAATGATCCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc  <  1:218532/138‑1 (MQ=255)
                                                                                                                                 aaaaaTTTAAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTTAGTAGGAGAAAGGAATTATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc  <  6:266222/138‑1 (MQ=255)
                                                                                                                                 aaaaaTTTAAAAATAAAAAATTTGTCCATTTTAATAACAATTTTGTTTTAATTTTTGAAATTTTATATAAAATATTGAGTAGGACAAAGGAATGATCCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc  <  3:201258/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
CCTATGTTCTGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAC  >  BS168_reference_genome/1023015‑1023279

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: