Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 783,336 T→A L274Q (CTG→CAG)  ybhI → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606783,3360TA97.7% 683.1 / 3.9 174L274Q (CTG→CAG) ybhIpredicted transporter
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base A (95/74);  total (98/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 185 total aligned reads displayed.

AGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCACCGT  >  REL606/783264‑783408
                                                                        |                                                                        
agGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATgcagct                                                                       >  1:4206472/1‑76 (MQ=255)
agGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATgcagct                                                                       >  2:1227116/1‑76 (MQ=255)
agGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATgcagct                                                                       >  2:3758492/1‑76 (MQ=255)
  gTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAGAATCCGCGAAAAAATGCAGCTCg                                                                     >  1:1918290/1‑76 (MQ=255)
   tGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCgg                                                                    >  2:1068498/1‑76 (MQ=255)
   tGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCgg                                                                    >  2:4738516/1‑76 (MQ=255)
     gATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGgtg                                                                  >  2:173115/1‑76 (MQ=255)
     gATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGgtg                                                                  <  1:3432848/76‑1 (MQ=255)
     gATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGgtg                                                                  <  1:2123826/76‑1 (MQ=255)
      aTAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGgtgt                                                                 >  1:6014003/1‑76 (MQ=255)
      aTAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGgtgt                                                                 <  2:4050033/76‑1 (MQ=255)
      aTAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGgtgt                                                                 >  1:1843129/1‑76 (MQ=255)
       tAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTc                                                                >  1:4378881/1‑76 (MQ=255)
       tAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTc                                                                >  2:4637578/1‑76 (MQ=255)
         aCAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCtt                                                              >  1:1121758/1‑76 (MQ=255)
         aCAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCtt                                                              <  2:4373622/76‑1 (MQ=255)
          cAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGAGCGATGAAAATCCGCGAAAAGATGCTGCTCCGTGTCttt                                                             <  2:5064511/76‑1 (MQ=255)
          cAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCttt                                                             >  1:4725659/1‑76 (MQ=255)
          cAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCttt                                                             <  2:2627507/76‑1 (MQ=255)
          cAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCttt                                                             <  2:2910499/76‑1 (MQ=255)
          cAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCttt                                                             <  2:5597302/76‑1 (MQ=255)
           aaaaCCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTtg                                                            >  2:1837884/1‑76 (MQ=255)
           aaaaCCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTtg                                                            >  2:3923940/1‑76 (MQ=255)
           aaaaCCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTtg                                                            <  1:2418713/76‑1 (MQ=255)
            aaaCCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTtgt                                                           >  1:2504561/1‑76 (MQ=255)
             aaCCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTtgtg                                                          <  1:1662043/76‑1 (MQ=255)
             aaCCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTtgtg                                                          <  2:4376414/76‑1 (MQ=255)
               ccATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCt                                                        >  2:5785202/1‑76 (MQ=255)
               ccATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCt                                                        <  2:4618857/76‑1 (MQ=255)
                cATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTg                                                       >  1:4546954/1‑76 (MQ=255)
                  tCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGgc                                                     <  1:1433663/76‑1 (MQ=255)
                    gCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCGGCGAAAAAATGCAGCTCGGTGGCTTTGTGCTGgcgc                                                   >  1:320478/1‑76 (MQ=255)
                    gCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTCGTGCTGgcgc                                                   >  2:5771696/1‑76 (MQ=255)
                    ccTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGgcgc                                                   >  2:868437/2‑76 (MQ=255)
                     cTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgct                                                  >  2:386585/1‑76 (MQ=255)
                     cTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgct                                                  >  1:5195323/1‑76 (MQ=255)
                     cTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgct                                                  >  1:4757525/1‑76 (MQ=255)
                     cTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgct                                                  >  2:3151174/1‑76 (MQ=255)
                      tAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctg                                                 >  1:4965945/1‑76 (MQ=255)
                      tAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctg                                                 >  1:1794077/1‑76 (MQ=255)
                      tAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctg                                                 >  2:3211454/1‑76 (MQ=255)
                      tAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctg                                                 <  2:3033698/76‑1 (MQ=255)
                       aaaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgc                                                >  2:5187459/1‑76 (MQ=255)
                       aaaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgc                                                >  2:3731726/1‑76 (MQ=255)
                       aaaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgc                                                >  1:324568/1‑76 (MQ=255)
                       aaaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgc                                                >  1:5103906/1‑76 (MQ=255)
                       aaaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgc                                                <  1:6273592/76‑1 (MQ=255)
                        aaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgct                                               >  2:3254002/1‑76 (MQ=255)
                        aaGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCgctgct                                               <  2:4206472/76‑1 (MQ=255)
                         aGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGctgctg                                              <  2:4710870/76‑1 (MQ=255)
                           cGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTgctggg                                            <  2:6028127/76‑1 (MQ=255)
                            gggCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTgctgggc                                           >  1:3851156/1‑76 (MQ=255)
                              gCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTgctgggctg                                         <  2:5094476/76‑1 (MQ=255)
                               ccTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGTGCTgctgggttgg                                        >  2:1585186/1‑76 (MQ=255)
                               ccTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTgctgggctgg                                        >  1:4847635/1‑76 (MQ=255)
                                cTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGa                                       <  2:1619584/76‑1 (MQ=255)
                                cTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGa                                       >  1:3710558/1‑76 (MQ=255)
                                cTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGa                                       >  2:2298225/1‑76 (MQ=255)
                                 ttGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGAt                                      <  1:4099930/76‑1 (MQ=255)
                                 ttGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGAt                                      >  1:2686327/1‑76 (MQ=255)
                                  tGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGAtt                                     <  1:2480354/76‑1 (MQ=255)
                                   gCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGAttt                                    <  2:3940626/76‑1 (MQ=255)
                                    ccGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGAtttt                                   >  1:5043808/1‑76 (MQ=255)
                                     cGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTc                                  >  2:1349642/1‑76 (MQ=255)
                                     cGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTc                                  >  2:5960363/1‑76 (MQ=255)
                                     cGAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTc                                  >  1:85019/1‑76 (MQ=255)
                                      gAACTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCa                                 <  1:6201598/76‑1 (MQ=255)
                                       aaCTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAg                                <  1:2103150/76‑1 (MQ=255)
                                       aaCTAGGACCGAGGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCat                                >  2:2728727/1‑75 (MQ=255)
                                        aCTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGt                               <  1:4122056/76‑1 (MQ=255)
                                         cTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTa                              >  2:34509/1‑76 (MQ=255)
                                         cTAGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTa                              >  1:753528/1‑76 (MQ=255)
                                           aGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAg                            <  2:3157029/76‑1 (MQ=255)
                                           aGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAg                            >  1:534391/1‑76 (MQ=255)
                                           aGGACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAg                            >  1:6153343/1‑76 (MQ=255)
                                            ggACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGt                           <  1:2866103/76‑1 (MQ=255)
                                            ggACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGt                           >  1:4658756/1‑76 (MQ=255)
                                            ggACCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGt                           >  2:5490903/1‑76 (MQ=255)
                                              aCCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGtct                         >  2:2815986/1‑76 (MQ=255)
                                              aCCGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGtct                         >  2:4836634/1‑76 (MQ=255)
                                               ccGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGtctc                        <  2:1237440/76‑1 (MQ=255)
                                                cGATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTctct                       <  2:1848458/76‑1 (MQ=255)
                                                 gATGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTg                      <  1:1624283/76‑1 (MQ=255)
                                                   tGAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTggg                    <  2:6200304/76‑1 (MQ=255)
                                                    gAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTgggg                   >  1:3148339/1‑76 (MQ=255)
                                                    gAAAATCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTgggg                   <  2:5305109/76‑1 (MQ=255)
                                                     aaaaTCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTggggg                  >  2:5417720/1‑76 (MQ=255)
                                                     aaaaTCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTggggg                  >  1:3109698/1‑76 (MQ=255)
                                                      aaaTCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGt                 >  2:195238/1‑76 (MQ=255)
                                                      aaaTCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGt                 <  1:2528924/76‑1 (MQ=255)
                                                         tCCGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTtga              <  2:2230886/76‑1 (MQ=255)
                                                          ccGCGAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTtgat             <  1:667125/76‑1 (MQ=255)
                                                             cgAAAAAATGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGaa          <  1:1314838/76‑1 (MQ=255)
                                                               aaaaaaTGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATc        >  1:1170492/1‑76 (MQ=255)
                                                                aaaaaTGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATcc       <  1:2407366/76‑1 (MQ=255)
                                                                aaaaaTGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATcc       <  1:262354/76‑1 (MQ=255)
                                                                aaaaaTGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATcc       <  2:4634910/76‑1 (MQ=255)
                                                                 aaaaTGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCa      >  2:4335075/1‑76 (MQ=255)
                                                                    aTGCAGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCACCg   <  2:863358/76‑1 (MQ=255)
                                                                     tgcagcTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCACCGt  >  1:572827/1‑76 (MQ=255)
                                                                        |                                                                        
AGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCACCGT  >  REL606/783264‑783408

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