Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 104,346 G→T 5.8% D259Y (GAT→TAT)  clpC → class III stress response‑related ATPase, AAA+ superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome104,3460GT 5.8% 228.6 / 2.9 69D259Y (GAT→TAT) clpCclass III stress response‑related ATPase, AAA+ superfamily
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base T (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GGTAAAACGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGATGCATCTAATATTTT  >  BS168_reference_genome/104208‑104470
                                                                                                                                          |                                                                                                                            
ggTAAAACGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAgg                                                                                                                              <  1:250051/139‑1 (MQ=255)
     aaCGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGc                                                                                                                         <  5:256288/139‑1 (MQ=255)
     aaCGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGc                                                                                                                         >  6:256288/1‑139 (MQ=255)
     aaCGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGc                                                                                                                         >  4:240223/1‑139 (MQ=255)
     aaCGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGc                                                                                                                         >  2:7837/1‑139 (MQ=255)
     aaCGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGAGAACAGTTGTTGCAGGCACAAAATACCGCGGAGAATTTGAGGAACGc                                                                                                                         >  8:158015/1‑139 (MQ=255)
      aCGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGcc                                                                                                                        >  8:160131/1‑139 (MQ=255)
          cTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTg                                                                                                                      >  8:23561/1‑137 (MQ=255)
            aTCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGCCATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGTATCGCCTgaag                                                                                                                   >  3:290922/1‑138 (MQ=255)
            aTCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGTATCGCCTgaag                                                                                                                   >  1:297340/1‑138 (MQ=255)
             tCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTgaaga                                                                                                                  <  5:222418/138‑1 (MQ=255)
                cAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGTATCGCCTGAAGAAGGt                                                                                                              <  4:290922/139‑1 (MQ=255)
                cAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGTATCGCCTGAAGAAGGt                                                                                                              <  2:297340/139‑1 (MQ=255)
                   aaGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCAt                                                                                                           >  8:82960/1‑139 (MQ=255)
                          ggCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGa                                                                                                     >  8:234838/1‑138 (MQ=255)
                               aGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCg                                                                                               >  5:145017/1‑139 (MQ=255)
                               aGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCg                                                                                               >  5:143542/1‑139 (MQ=255)
                                 cAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGcc                                                                                             >  4:66210/1‑139 (MQ=255)
                                 cAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGc                                                                                              >  7:67047/1‑138 (MQ=255)
                                   aaTTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCcag                                                                                           <  1:95999/139‑1 (MQ=255)
                                       aTCAATAATGTAGTACCGGAAATTTTGCGTGATAAACGTGTGATGACATTAGACTTTGGAAGTGATGTTGGCGGAACAAAATACCGCTGATAATTTGATGTTAGCCTGAAGAAGGTCATGGATGAAATTCGCcaggca                                                                                        <  8:288904/138‑1 (MQ=255)
                                        tCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCcaggcagg                                                                                      <  7:234838/139‑1 (MQ=255)
                                         cAATAATGAAGTACCCGCAATTTTGCGTGATAAACGTGGGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGa                                                                                     <  2:297742/139‑1 (MQ=255)
                                         cAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCcaggcagg                                                                                      >  5:64657/1‑138 (MQ=255)
                                           ataatGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGaaa                                                                                   >  1:114590/1‑139 (MQ=255)
                                           ataatGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGaaa                                                                                   <  4:97817/139‑1 (MQ=255)
                                             aatGAAGTACACGAAATTTTGAGTGTTAAACGTGTTATGACATTAGACATGGGAACAGTTGTTGCCGGAACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAata                                                                                 <  2:84376/139‑1 (MQ=255)
                                              atGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAata                                                                                 <  4:215415/138‑1 (MQ=255)
                                                 aaGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCa                                                                              <  1:221533/138‑1 (MQ=255)
                                                  aGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCAtt                                                                            <  8:67047/139‑1 (MQ=255)
                                                    tACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGATAATTTGAGGATCGCCTGAAGAAGTTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCt                                                                          <  8:91332/139‑1 (MQ=255)
                                                      cccGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTAt                                                                        >  1:40509/1‑139 (MQ=255)
                                                      cccGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTAt                                                                        >  3:275156/1‑139 (MQ=255)
                                                         gAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCa                                                                     <  7:160131/139‑1 (MQ=255)
                                                               ttGCTTGATAAACGTGTGATTAAATTAGACATGGGAACAGTTGTTGCCGGCACAAACTACCGCTGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTAATCATCGAt                                                                <  8:245707/138‑1 (MQ=255)
                                                               ttGCGTGATATACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGAt                                                                <  4:58935/138‑1 (MQ=255)
                                                               ttGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGAt                                                                <  8:122084/138‑1 (MQ=255)
                                                                  cGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAATGTCATGGATGAAATTCGCCACGCAGGAAATATCATTCTATTCATCGATGAg                                                             <  3:157692/138‑1 (MQ=255)
                                                                  cGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAg                                                             <  1:28485/138‑1 (MQ=255)
                                                                  cGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGc                                                            <  7:82960/139‑1 (MQ=255)
                                                                       tAAACGTGTGATTACATTATACATTGGAACAGTTTTTGCCTTAACAAAATAACGCTGAGAATTTGAGGATCGCCTGAAGAAGGTAATGTATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCa                                                        <  4:35449/138‑1 (MQ=255)
                                                                               tgatgaCATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAAt                                               >  7:161016/1‑139 (MQ=255)
                                                                               tgatgaCATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAAt                                               <  3:201395/139‑1 (MQ=255)
                                                                                     cATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGc                                         <  7:288128/139‑1 (MQ=255)
                                                                                        tAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGTCAGGCAGGaaa                                                                                   >  6:191953/1‑94 (MQ=255)
                                                                                        tAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGTCAGGCAGGaaa                                                                                   <  5:191953/94‑1 (MQ=255)
                                                                                         aGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAgg                                      <  5:171340/138‑1 (MQ=255)
                                                                                         aGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAgg                                      <  7:23561/138‑1 (MQ=255)
                                                                                          gACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGc                                     <  1:7837/138‑1 (MQ=255)
                                                                                             aTGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAg                                 <  6:145017/139‑1 (MQ=255)
                                                                                               gggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAg                                 >  3:97810/1‑137 (MQ=255)
                                                                                               gggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGct                               >  4:181360/1‑138 (MQ=255)
                                                                                                ggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCcaggc                                                                                         >  6:157435/1‑80 (MQ=255)
                                                                                                ggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCcaggc                                                                                         <  5:157435/80‑1 (MQ=255)
                                                                                                ggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAg                                 <  4:14883/136‑1 (MQ=255)
                                                                                                ggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAg                                 >  3:14883/1‑136 (MQ=255)
                                                                                                ggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAg                                 >  7:278829/1‑136 (MQ=255)
                                                                                                ggAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAg                                 >  2:85787/1‑136 (MQ=255)
                                                                                                  aaCAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAg                              >  4:10454/1‑137 (MQ=255)
                                                                                                    gaGTTGTTGCCGGCACAAAATCCCGCGGATAATTTGAGGATCGCCTGAAGAAGGTCATGGATTAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAgg                          <  2:40509/138‑1 (MQ=255)
                                                                                                    cAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAgg                          <  4:275156/139‑1 (MQ=255)
                                                                                                    cAGATGTTGCCGGCACAAAATACCGCGTAGAATTTGAGGATCGCCTGTAGAATGTAATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAgg                          <  2:205519/139‑1 (MQ=255)
                                                                                                            gCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTAtt                   >  3:230441/1‑138 (MQ=255)
                                                                                                              cGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGAt                <  3:120024/139‑1 (MQ=255)
                                                                                                              cGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGAt                <  7:228390/139‑1 (MQ=255)
                                                                                                                 cacaAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTc                                                                                                             >  5:116161/1‑43 (MQ=255)
                                                                                                                 cacaAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTc                                                                                                             <  6:116161/43‑1 (MQ=255)
                                                                                                                      aaTACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGATGCATCTaa        <  2:114590/139‑1 (MQ=255)
                                                                                                                           cgcgGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGATGCATCTAATAttt   >  3:154518/1‑139 (MQ=255)
                                                                                                                             cgGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGATGCATCTAATAtttt  >  5:148876/1‑138 (MQ=255)
                                                                                                                                          |                                                                                                                            
GGTAAAACGGCTATCGCAGAAGGTTTGGCACAGCAAATTATCAATAATGAAGTACCCGAAATTTTGCGTGATAAACGTGTGATGACATTAGACATGGGAACAGTTGTTGCCGGCACAAAATACCGCGGAGAATTTGAGGATCGCCTGAAGAAGGTCATGGATGAAATTCGCCAGGCAGGAAATATCATTCTATTCATCGATGAGCTCCATACATTAATCGGGGCAGGCGGAGCAGAAGGTGCTATTGATGCATCTAATATTTT  >  BS168_reference_genome/104208‑104470

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: