Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 979,852 T→G S256A (TCC→GCC)  rpsA → 30S ribosomal protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606979,8520TG98.6% 1506.5 / 2.1 352S256A (TCC→GCC) rpsA30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (164/182);  total (166/186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 359 total aligned reads displayed.

GAGCGAAATCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTATCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTG  >  REL606/979779‑979926
                                                                         |                                                                          
gAGCGAAATCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGcc                                                                          >  2:6596661/1‑76 (MQ=255)
    gAAATCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTgg                                                                      <  2:96033/76‑1 (MQ=255)
    gAAATCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTgg                                                                      <  1:9537966/76‑1 (MQ=255)
     aaaTCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTggg                                                                     <  1:3114282/76‑1 (MQ=255)
       aTCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGcc                                                                   >  1:2547505/1‑76 (MQ=255)
         cGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTg                                                                 >  1:7347165/1‑76 (MQ=255)
         cGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTg                                                                 <  1:8937520/76‑1 (MQ=255)
         cGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTg                                                                 <  2:5443730/76‑1 (MQ=255)
          gTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGa                                                                >  2:2006409/1‑76 (MQ=255)
          gTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGa                                                                <  1:3288429/76‑1 (MQ=255)
          gTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGa                                                                <  2:4324539/76‑1 (MQ=255)
          gTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGa                                                                <  2:9877073/76‑1 (MQ=255)
           tCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGaa                                                               <  1:2854350/76‑1 (MQ=255)
           tCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGaa                                                               <  2:8743548/76‑1 (MQ=255)
            cAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGaaa                                                              >  2:3881843/1‑76 (MQ=255)
            cAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGaaa                                                              <  2:7031940/76‑1 (MQ=255)
            cAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGaaa                                                              <  1:7631266/76‑1 (MQ=255)
             aaCGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAAc                                                             <  1:10199223/76‑1 (MQ=255)
             aaCGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAAc                                                             <  1:10052144/76‑1 (MQ=255)
               cGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAg                                                           >  2:9837402/1‑76 (MQ=255)
               cGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAg                                                           >  2:8550907/1‑76 (MQ=255)
                 tGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCt                                                         <  2:3493425/76‑1 (MQ=255)
                  gggCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTg                                                        >  2:6143653/1‑76 (MQ=255)
                    gCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTggg                                                      <  1:4498012/76‑1 (MQ=255)
                      gacgaAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTg                                                        >  1:7736804/1‑72 (MQ=255)
                        cgaAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGaa                                                  >  2:1019190/1‑76 (MQ=255)
                          aaaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGa                                                <  1:3509116/76‑1 (MQ=255)
                          aaaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGa                                                >  1:8131355/1‑76 (MQ=255)
                          aaaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGa                                                <  1:591484/76‑1 (MQ=255)
                           taTCAATGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGAt                                               >  2:4687912/2‑76 (MQ=255)
                           aaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGAt                                               <  1:6505812/76‑1 (MQ=255)
                           aaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGAt                                               >  1:6965829/1‑76 (MQ=255)
                           aaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGAt                                               <  1:374451/76‑1 (MQ=255)
                           aaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGAt                                               <  2:1089466/76‑1 (MQ=255)
                           aaTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGAt                                               >  1:8787207/1‑76 (MQ=255)
                            aTCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATc                                              <  2:6452021/76‑1 (MQ=255)
                             tCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATcc                                             >  2:1397508/1‑76 (MQ=255)
                               aCTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGt                                           <  1:1190459/76‑1 (MQ=255)
                               aCTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGt                                           <  2:10064501/76‑1 (MQ=255)
                                  gTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTggg                                        <  2:3789947/76‑1 (MQ=255)
                                  gTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTggg                                        <  1:707360/76‑1 (MQ=255)
                                   ttAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGt                                       <  2:7412145/76‑1 (MQ=255)
                                   ttAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGt                                       >  2:2619426/1‑76 (MQ=255)
                                    tAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTa                                      <  2:5684928/76‑1 (MQ=255)
                                    tAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTa                                      <  1:386261/76‑1 (MQ=255)
                                     aaaGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAg                                     >  1:194778/1‑76 (MQ=255)
                                     aaaGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAg                                     >  2:4174106/1‑76 (MQ=255)
                                       aGTGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCt                                   >  1:8525880/1‑76 (MQ=255)
                                         tGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTAt                                 <  2:5629390/76‑1 (MQ=255)
                                         tGCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTAt                                 <  1:2233481/76‑1 (MQ=255)
                                          gCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATc                                <  2:7824936/76‑1 (MQ=255)
                                          gCTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATc                                >  1:6329073/1‑76 (MQ=255)
                                           cTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCg                               >  1:5465594/1‑76 (MQ=255)
                                           cTGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCg                               <  1:4316381/76‑1 (MQ=255)
                                            tGAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGc                              >  1:2762398/1‑76 (MQ=255)
                                             gAAGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCt                             >  2:5038701/1‑76 (MQ=255)
                                              aaGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTa                            <  2:3964053/76‑1 (MQ=255)
                                              aaGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTa                            >  2:8424048/1‑76 (MQ=255)
                                               aGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaa                           <  2:3172562/76‑1 (MQ=255)
                                               aGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaa                           <  2:9279692/76‑1 (MQ=255)
                                               aGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaa                           >  2:8260652/1‑76 (MQ=255)
                                               aGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaa                           <  2:1529248/76‑1 (MQ=255)
                                               aGTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaa                           <  2:5171547/76‑1 (MQ=255)
                                                gTTCGACCGCGAACGTCCCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaaa                          <  2:9236306/76‑1 (MQ=255)
                                                gTTCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTaaa                          <  1:5649761/76‑1 (MQ=255)
                                                 ttCGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAAc                         >  1:170796/1‑76 (MQ=255)
                                                   cGACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGt                       >  2:2387988/1‑76 (MQ=255)
                                                    gACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGtt                      >  1:9356787/1‑76 (MQ=255)
                                                    gACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGtt                      >  1:1506724/1‑76 (MQ=255)
                                                    gACCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGtt                      >  2:4452262/1‑76 (MQ=255)
                                                     aCCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTa                     <  1:884936/76‑1 (MQ=255)
                                                     aCCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTa                     <  2:480906/76‑1 (MQ=255)
                                                     aCCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTa                     >  1:4039097/1‑76 (MQ=255)
                                                     aCCGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTa                     >  1:6852392/1‑76 (MQ=255)
                                                      ccGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTAt                    >  1:4823968/1‑76 (MQ=255)
                                                      ccGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTAt                    <  1:2390182/76‑1 (MQ=255)
                                                      ccGCGAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTAt                    <  2:2938325/76‑1 (MQ=255)
                                                        gcgAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATcc                  >  1:2121158/1‑76 (MQ=255)
                                                        gcgAACGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATcc                  >  2:816051/1‑76 (MQ=255)
                                                           aaCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGa               <  2:1861907/76‑1 (MQ=255)
                                                           aaCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGa               <  1:8636720/76‑1 (MQ=255)
                                                           aaCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGa               >  2:6374303/1‑76 (MQ=255)
                                                           aaCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGa               >  1:9680573/1‑76 (MQ=255)
                                                            aCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGaa              <  1:5213644/76‑1 (MQ=255)
                                                            aCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGaa              <  1:599514/76‑1 (MQ=255)
                                                            aCGTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGaa              <  2:215239/76‑1 (MQ=255)
                                                              gtccccGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAgg            <  2:1649432/73‑1 (MQ=255)
                                                              gTACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAgg            >  2:7287387/1‑76 (MQ=255)
                                                               tACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGt           >  1:9112780/1‑76 (MQ=255)
                                                               tACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGt           >  2:6740122/1‑76 (MQ=255)
                                                               tACCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGt           <  1:7942829/76‑1 (MQ=255)
                                                                aCCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTa          <  2:5735425/76‑1 (MQ=255)
                                                                aCCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTa          >  2:3701875/1‑76 (MQ=255)
                                                                aCCCGTGTAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTa          >  2:1292529/1‑76 (MQ=255)
                                                                   cGTGTAGCCGTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCa       <  2:8875976/76‑1 (MQ=255)
                                                                    gtgtAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCaa      <  2:6869732/76‑1 (MQ=255)
                                                                     tgtAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCaaa     >  2:7622984/1‑76 (MQ=255)
                                                                     tgtAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCaaa     <  1:6120922/76‑1 (MQ=255)
                                                                      gtAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAac    >  1:1328664/1‑76 (MQ=255)
                                                                      gtAGCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAac    <  1:4183625/76‑1 (MQ=255)
                                                                         |                                                                          
GAGCGAAATCGTCAACGTGGGCGACGAAATCACTGTTAAAGTGCTGAAGTTCGACCGCGAACGTACCCGTGTATCCCTGGGCCTGAAACAGCTGGGCGAAGATCCGTGGGTAGCTATCGCTAAACGTTATCCGGAAGGTACCAAACTG  >  REL606/979779‑979926

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