Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 209,957 G→T 5.8% G109C (GGT→TGT)  ybcF → putative enzyme; prophage 1 region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome209,9570GT 5.8% 233.3 / 4.3 69G109C (GGT→TGT) ybcFputative enzyme; prophage 1 region
Reads supporting (aligned to +/- strand):  ref base G (27/38);  new base T (2/2);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

ATGTTGAGGAGATTTTTGTAGTAGGAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGTTAATCATAAGGGTGGAAAA  >  BS168_reference_genome/209820‑210094
                                                                                                                                         |                                                                                                                                         
aTGTTGAGGAGATTTTTGTAGTAGGAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAgg                                                                                                                                          <  6:138895/139‑1 (MQ=255)
      aggagATTTTTGTAGTAGGAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTa                                                                                                                                     <  6:142630/138‑1 (MQ=255)
      aggagATTTTTGTAGTAGGAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTAc                                                                                                                                    >  2:147664/1‑139 (MQ=255)
          taTTTTTTTAGTATGAAAAAAAGTTACAAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGt                                                                                                                                 <  2:162365/137‑1 (MQ=255)
           aTTTTTGTAGTAGGAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTc                                                                                                                                >  1:162365/1‑138 (MQ=255)
                       gtaaaaaaaTATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTa                                                                                                                    <  2:228564/133‑1 (MQ=255)
                       ggAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTa                                                                                                                    >  3:213246/1‑138 (MQ=255)
                       ggAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTa                                                                                                                    >  3:167290/1‑138 (MQ=255)
                       ggAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTa                                                                                                                    >  1:356896/1‑138 (MQ=255)
                       ggAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTa                                                                                                                    >  5:411973/1‑138 (MQ=255)
                         aaCAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAAt                                                                                                                 <  1:201120/139‑1 (MQ=255)
                                   aagaagACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAAt                                                                                                           >  3:310266/1‑135 (MQ=255)
                                   aagaagACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAAt                                                                                                           <  4:310266/135‑1 (MQ=255)
                                   aagaagACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAg                                                                                                         >  7:254257/1‑137 (MQ=255)
                                       agaCTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGt                                                                                                     <  7:349806/137‑1 (MQ=255)
                                       agaCTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGa                                                                                                   <  8:12194/139‑1 (MQ=255)
                                       agaCTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGa                                                                                                   <  8:316628/139‑1 (MQ=255)
                                         aCTTCTGCAGGTAATGGACTAACTCAACTTGACACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGaca                                                                                                 <  4:356501/139‑1 (MQ=255)
                                           ttCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGacac                                                                                                >  2:307035/1‑138 (MQ=255)
                                                   gTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGaaaa                                                                                       <  4:283172/139‑1 (MQ=255)
                                                   gTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGaaaa                                                                                       <  2:356896/139‑1 (MQ=255)
                                                   gTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGaaaa                                                                                       <  5:8062/139‑1 (MQ=255)
                                                     aaTGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGaaaa                                                                                       >  4:338452/1‑137 (MQ=255)
                                                         gACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAgtgt                                                                                 <  6:410684/139‑1 (MQ=255)
                                                             aaCTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGa                                                                              >  8:44147/1‑138 (MQ=255)
                                                             aaCTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGGTGat                                                                             >  1:136717/1‑139 (MQ=255)
                                                              aCTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATTGCAGTGACACCATTGAAAAAAGTGTTGata                                                                            <  7:380168/139‑1 (MQ=255)
                                                              aCTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGata                                                                            >  7:175492/1‑139 (MQ=255)
                                                                   aCTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGata                                                                            <  3:402754/134‑1 (MQ=255)
                                                                   aCTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGata                                                                            >  4:402754/1‑134 (MQ=255)
                                                                    cTTGAAACAATGAAAGATAAACTACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTc                                                                      <  2:157099/139‑1 (MQ=255)
                                                                    cTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTc                                                                      >  6:312448/1‑139 (MQ=255)
                                                                    cTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTc                                                                      >  7:154166/1‑139 (MQ=255)
                                                                    cTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTc                                                                      >  4:230310/1‑139 (MQ=255)
                                                                    cTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTc                                                                      <  4:26440/139‑1 (MQ=255)
                                                                    cTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGAGGAATGGCGAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATTTGCTTc                                                                      >  2:116611/1‑139 (MQ=255)
                                                                         aaCAATGAAAGATAAAATACAAACATTTGATTATCTGTTTCAAAATTTCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAca                                                                  <  5:204031/138‑1 (MQ=255)
                                                                         aaCAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAca                                                                  <  4:292145/138‑1 (MQ=255)
                                                                         aaCAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAca                                                                  <  2:225143/138‑1 (MQ=255)
                                                                           cAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatca                                                               <  7:44147/139‑1 (MQ=255)
                                                                           cAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatca                                                               <  1:147664/139‑1 (MQ=255)
                                                                           cAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatca                                                               <  4:167290/139‑1 (MQ=255)
                                                                           cAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatca                                                               <  4:213246/139‑1 (MQ=255)
                                                                            aaTGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatcat                                                              <  3:118101/139‑1 (MQ=255)
                                                                              tGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTATGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGat                                                                             <  7:390990/122‑1 (MQ=255)
                                                                              tGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTATGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGat                                                                             <  3:397575/122‑1 (MQ=255)
                                                                              tGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTATGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGat                                                                             >  4:397575/1‑122 (MQ=255)
                                                                              tGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTATGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGat                                                                             >  8:390990/1‑122 (MQ=255)
                                                                                 aaGATAAAATACATAGATTGGTTTATCTGTTTGAAAATTGCAAGCATTAATTTTTAGGTGGTTCGGTCGATGATTGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGTTTCGACATCATcccc                                                          <  2:179878/138‑1 (MQ=255)
                                                                                    aTAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGt                                                      <  1:15375/139‑1 (MQ=255)
                                                                                     tAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTg                                                     <  5:82906/139‑1 (MQ=255)
                                                                                        aaTACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCg                                                  >  1:264197/1‑139 (MQ=255)
                                                                                          tACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTc                                                <  3:230310/139‑1 (MQ=255)
                                                                                            cAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCat                                              <  3:103271/139‑1 (MQ=255)
                                                                                                 aTTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGtt                                         <  8:175492/139‑1 (MQ=255)
                                                                                                    ggATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTaaa                                      >  4:330272/1‑139 (MQ=255)
                                                                                                         aTCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTa                                  <  3:98305/138‑1 (MQ=255)
                                                                                                          tCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAca                                                                  >  6:131420/1‑105 (MQ=255)
                                                                                                          tCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAca                                                                  <  5:131420/105‑1 (MQ=255)
                                                                                                          tCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTa                                  >  8:61916/1‑137 (MQ=255)
                                                                                                          tCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTa                                  >  8:176437/1‑137 (MQ=255)
                                                                                                           cTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatca                                                               >  8:7629/1‑107 (MQ=255)
                                                                                                           cTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGAcatca                                                               <  7:7629/107‑1 (MQ=255)
                                                                                                               ttgaaaaTTGCAAGCCTGAATTTTTAGGTGGTCCGGTCGCGGAAACGCTAAATGAAAATAGCAGTGACTCCAGTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGGCATATGTTTAAGTTATAGGttt                           >  3:193042/4‑139 (MQ=255)
                                                                                                                      ttGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGtt                    <  1:54444/139‑1 (MQ=255)
                                                                                                                      ttGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGtt                    <  1:112059/139‑1 (MQ=255)
                                                                                                                       tGCAAGCCTTAATTTTTAGGTTGTACGGTAGATGTATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGtt                    <  8:254257/138‑1 (MQ=255)
                                                                                                                        gCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTGATTTGTTTAATTTATAGGGTTATTAGttta                  >  1:256772/1‑137 (MQ=255)
                                                                                                                            gCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGTTAATCAt              <  3:330272/139‑1 (MQ=255)
                                                                                                                                        aGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGTTAATCATAAGGGTGGaaaa  >  5:349558/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                         
ATGTTGAGGAGATTTTTGTAGTAGGAACAAAAGATAAGAAGACTTCTGCAGGTAATGGACTAACTCAACTTGAAACAATGAAAGATAAAATACAAACATTGGATTATCTGTTTCAAAATTGCAAGCCTGAATTTTTAGGTGGTACGGTCGATGAATGGCTAAATGAAAATAGCAGTGACACCATTGAAAAAAGTGTTGATATGATTCGACATCATCCCCTAGTGCCGTCATATGTTAAAGTTAGAGGTTTATTCGTTAATCATAAGGGTGGAAAA  >  BS168_reference_genome/209820‑210094

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