Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 199,191 A→C 13.2% N232T (AAC→ACC)  glmM → phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome199,1910AC13.2% 150.5 / 2.7 68N232T (AAC→ACC) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (2/8);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 9.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTTGGCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGGA  >  BS168_reference_genome/199055‑199326
                                                                                                                                        |                                                                                                                                       
tCCTTGGCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAAcg                                                                                                                                       <  5:166399/139‑1 (MQ=255)
tCCTTGGCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAAcg                                                                                                                                       >  3:63772/1‑139 (MQ=255)
tCCTTGGCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAAcg                                                                                                                                       >  3:183779/1‑139 (MQ=255)
  cTTGGCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAAcgcg                                                                                                                                     <  4:183779/139‑1 (MQ=255)
     ggCGACACACCTTTTTGCTGAGTTAGATGCAGATGTTTCTACAATGGGGAATTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGAATCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGa                                                                                                                                   >  8:7401/1‑138 (MQ=255)
     ggCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGCAAAACGCGGAt                                                                                                                                  >  5:32428/1‑139 (MQ=255)
     ggCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAACCGCGGa                                                                                                                                   >  1:353969/1‑138 (MQ=255)
     ggCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGAt                                                                                                                                  >  7:352177/1‑139 (MQ=255)
       cGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGAtct                                                                                                                                <  7:46856/139‑1 (MQ=255)
       cGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGGCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGAtct                                                                                                                                <  8:105998/139‑1 (MQ=255)
         acacacCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGAtctct                                                                                                                              >  1:351967/1‑138 (MQ=255)
              ccTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGACAAAACCGCGGATCTCGGTCtt                                                                                                                         <  5:291311/139‑1 (MQ=255)
                     gCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTc                                                                                                                   <  8:42344/138‑1 (MQ=255)
                       tGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGCAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGAc                                                                                                                >  1:258404/1‑139 (MQ=255)
                       tGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGAc                                                                                                                >  1:194352/1‑139 (MQ=255)
                       tGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGAc                                                                                                                >  1:210617/1‑139 (MQ=255)
                         aTTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGCTCTCGGTCTTGCGTTCGACgg                                                                                                              >  5:355695/1‑139 (MQ=255)
                             aGATGCAGATGTTTCAACAATGGGGACTTCCCCGAACGGATTAACCATTAATGACGGCGGCGGTGCGACGCATCCCGACGCGCGCAGCGCGTTTGTCAAAGAGAAAACCGCGGATCTCGGTCTTGCGTTCGACGGTGa                                                                                                           <  3:122053/138‑1 (MQ=255)
                                   aGATGTTTCTACATTGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGAc                                                                                                    <  4:352939/139‑1 (MQ=255)
                                   aGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGAc                                                                                                    <  3:182432/139‑1 (MQ=255)
                                          tCTACAATGGGGACTTCCCCGAACCGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGc                                                                                                 >  8:111948/1‑135 (MQ=255)
                                              cattgtggACTTCCCCTTACGTATATAACATTAATGACGGCGGGGGGGCGTCTAACCCCGATGGGCTCAGGGCTTTTGTCATATCTAAAAAGGCTGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTg                                                                                          <  8:173899/132‑1 (MQ=255)
                                                 tGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGt                                                                                      >  8:334639/1‑139 (MQ=255)
                                                  ggggACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTc                                                                                     >  5:41465/1‑139 (MQ=255)
                                                        ttCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGa                                                                                >  2:94971/1‑138 (MQ=255)
                                                        ttCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGCGACGGCGACCGCCTGATTGCTGTCGATGa                                                                                >  8:290552/1‑138 (MQ=255)
                                                         tCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTATGTCAAACACAAAACCTCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCAGAATGCGGTCGATGaa                                                                               >  5:235051/1‑138 (MQ=255)
                                                          ccccTAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGaaaa                                                                             <  4:61703/139‑1 (MQ=255)
                                                           cccGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGaaaaa                                                                            <  7:301635/139‑1 (MQ=255)
                                                             cGAACGGTTTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCGCAGCGCTTTTGTCAAAGAGAAAACCGCGGTTCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAgg                                                                          <  6:180525/139‑1 (MQ=255)
                                                             cGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAgg                                                                          <  8:164734/139‑1 (MQ=255)
                                                             cGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAgg                                                                          <  7:188629/139‑1 (MQ=255)
                                                               aaCGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAGGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGa                                                                         >  5:83620/1‑138 (MQ=255)
                                                               aaCGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGCGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGaa                                                                        <  7:290552/139‑1 (MQ=255)
                                                                  ggATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAata                                                                     <  2:194352/139‑1 (MQ=255)
                                                                    aTTAAACATTAATGACGGCGTCGGTTCGACTCGTCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATAtt                                                                   >  8:298548/1‑139 (MQ=255)
                                                                             tAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCGGATTGCTGTCGATGAATAAGGAAATATTGTAGACGGc                                                          >  2:249410/1‑139 (MQ=255)
                                                                             tAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCACCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGc                                                          >  4:304425/1‑139 (MQ=255)
                                                                                  ccGGCGTCGGTTCGACTCATCCCGAAGCGCCCAGCGCGTTTGTCAAAGACAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGAcc                                                      <  5:10772/137‑1 (MQ=255)
                                                                                    ggggtgggTTCGACTCATCCCGAAGCGCTCAGGGCGTTTGTAAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCaaa                                                   <  8:352177/133‑1 (MQ=255)
                                                                                    ggcggcGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGAGGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCaaa                                                   <  6:60081/139‑1 (MQ=255)
                                                                                    ggCGTGGGTTGGACTCATCCCTTAGCGCTCAGCGCGTTTGTCAAAGAGAAAACGGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATTTTGTAGACGGCGACCaaa                                                   <  5:279571/139‑1 (MQ=255)
                                                                                    ggCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATTAAAAAGGAAAAATTTTAGACGGCGCCCac                                                    >  2:170292/1‑137 (MQ=255)
                                                                                    ggCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGAcc                                                      >  4:350215/1‑136 (MQ=255)
                                                                                    ggCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGAcc                                                      <  3:350215/136‑1 (MQ=255)
                                                                                    ggCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCaaa                                                   <  4:157481/139‑1 (MQ=255)
                                                                                    ggCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCaaa                                                   <  6:181304/139‑1 (MQ=255)
                                                                                    ggCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTATTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCaaa                                                   >  3:214181/1‑139 (MQ=255)
                                                                                      cGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGGCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATc                                                 <  8:114990/139‑1 (MQ=255)
                                                                                         cGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAACCGCGGATCTCGGGCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATg                                              <  4:343320/139‑1 (MQ=255)
                                                                                                aCTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACatat                                       >  1:18659/1‑139 (MQ=255)
                                                                                                      cccGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGACTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGaa                                                                        <  4:148382/100‑1 (MQ=255)
                                                                                                      cccGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGACTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGaa                                                                        >  3:148382/1‑100 (MQ=255)
                                                                                                        cGACGCGCTCAGCGCGTTTGTAAAAGAGAAAACCGCGGATCTCGGTCTTGCGTTCGACGGGGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCaaaa                               <  3:120374/139‑1 (MQ=255)
                                                                                                           aGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACAt                            >  5:366408/1‑139 (MQ=255)
                                                                                                               cTCCGCGCTTTTTGCAACGATAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGa                        <  6:173984/139‑1 (MQ=255)
                                                                                                               cTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCaa                                 >  7:17115/1‑130 (MQ=255)
                                                                                                               cTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCaa                                 <  8:17115/130‑1 (MQ=255)
                                                                                                                tCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGAcc                                                      >  1:81872/1‑108 (MQ=255)
                                                                                                                tCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGAcc                                                      <  2:81872/108‑1 (MQ=255)
                                                                                                                    ggggTTTGACAAATAAAAAACCGGGGAGCTCGGTCTGGCGTGGGAGGGTGAGGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCa                   <  4:8468/136‑1 (MQ=255)
                                                                                                                    cgcgTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAAt                     >  7:347538/1‑137 (MQ=255)
                                                                                                                          tGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGc             >  7:249442/1‑139 (MQ=255)
                                                                                                                              aaaGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGtt         <  8:347538/139‑1 (MQ=255)
                                                                                                                                ggagaAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTa       <  6:83620/138‑1 (MQ=255)
                                                                                                                                agagAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTa       <  3:156641/139‑1 (MQ=255)
                                                                                                                                  cgaAAACCGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTTATTGCTGTCGATGAAAAAGTAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTaaa     <  5:79836/138‑1 (MQ=255)
                                                                                                                                   gaAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTaaa     <  6:366408/138‑1 (MQ=255)
                                                                                                                                     aacaCGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGga  <  5:37357/136‑1 (MQ=255)
                                                                                                                                     aaatCGCGGATCTAGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGga  <  7:111948/139‑1 (MQ=255)
                                                                                                                                     aaacCGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTTATTGCTGTCGATGAAAAAGTAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGga  <  4:214181/139‑1 (MQ=255)
                                                                                                                                     aaaaCGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGga  <  2:18659/139‑1 (MQ=255)
                                                                                                                                     aaaaCGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGga  <  6:41465/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TCCTTGGCGACACACCTGTTTGCTGATTTAGATGCAGATGTTTCTACAATGGGGACTTCCCCGAACGGATTAAACATTAATGACGGCGTCGGTTCGACTCATCCCGAAGCGCTCAGCGCGTTTGTCAAAGAGAAAAACGCGGATCTCGGTCTTGCGTTCGACGGTGACGGCGACCGCCTGATTGCTGTCGATGAAAAAGGAAATATTGTAGACGGCGACCAAATCATGTACATATGCTCAAAACATCTGAAATCAGAGGGCCGTTTAAAGGA  >  BS168_reference_genome/199055‑199326

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: