Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,192,6250TA44.1% 77.6 / 203.5 180intergenic (+84/+336)b3051/b3052putative membrane protein/putative kinase
Reads supporting (aligned to +/- strand):  ref base T (56/44);  new base A (50/29);  total (106/74)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 193 total aligned reads displayed.

CCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCT  >  NC_000913/3192550‑3192700
                                                                           |                                                                           
ccAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt                                                                             >  1:2655440/1‑76 (MQ=255)
 cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt                                                                            <  2:1994045/76‑1 (MQ=255)
 cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt                                                                            <  2:5030845/76‑1 (MQ=255)
    gAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgtt                                                                         <  1:2094773/76‑1 (MQ=255)
      ataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgc                                                                       >  1:4897260/1‑76 (MQ=255)
       taatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgca                                                                      >  1:1646894/1‑76 (MQ=255)
       taatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgca                                                                      >  1:5835717/1‑76 (MQ=255)
       atATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCa                                                                      >  2:2269277‑M1/3‑76 (MQ=255)
       atATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCa                                                                      <  1:607175‑M1/74‑1 (MQ=255)
        aatGCCTGGCGTAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcag                                                                     <  2:1616738/76‑1 (MQ=255)
         atCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTgcagc                                                                    <  2:4625058/76‑1 (MQ=255)
         agcccTGGCCCAAAACGCCGGGGATTTTTTATAGCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcagc                                                                    <  2:5014452/74‑1 (MQ=255)
           cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTgcagcgt                                                                  >  1:5207324/1‑74 (MQ=255)
             cTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGttt                                                                >  1:5041498/1‑76 (MQ=255)
                 cccAAAACGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCg                                                            <  1:2306890/76‑1 (MQ=255)
                 cccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCg                                                            >  2:3446698/1‑76 (MQ=255)
                 cccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCg                                                            <  2:942101/76‑1 (MQ=255)
                    aaaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGa                                                         <  1:4101497/76‑1 (MQ=255)
                    aaaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGa                                                         <  1:4329506/76‑1 (MQ=255)
                     aaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAc                                                        >  2:2632219/1‑76 (MQ=255)
                     aaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGAc                                                        <  2:4460866/76‑1 (MQ=255)
                      aaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAcc                                                       <  2:5086161/76‑2 (MQ=255)
                      aaaGCCGGGGATTTTTTATATCTGCGTTCCGCGAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACg                                                       <  2:232189/76‑1 (MQ=255)
                      aaaGCCGGGGATTTTTTATATCTGCGCTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAcc                                                       >  1:3226906/1‑75 (MQ=255)
                        aGCCGGGGATTTTTTCTATCTGCGTTCCGCTAAAAAGTGCCAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                                                     <  2:5721212/76‑1 (MQ=255)
                          ccGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTc                                                   <  1:1013703/76‑1 (MQ=255)
                          ccGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTc                                                   >  1:819524/1‑76 (MQ=255)
                           cGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGGTGCAGCGTTTTGCGTGACCGCTCg                                                  >  2:2688706/1‑76 (MQ=255)
                            ggggATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCgg                                                 <  2:5797169/76‑1 (MQ=255)
                            ggggATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCgg                                                 >  2:3987077/1‑76 (MQ=255)
                            ggggATGTTTTATATCTGCGTTCAGCTAAAAGGTGCAAACGCTCAGGACGGTGCACCGTTTTGCGCGACCGCTCgg                                                 >  2:674375/1‑76 (MQ=255)
                              gtaTTTTGTATATCTGCGGTACGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCgggn                                               <  1:3136036/74‑2 (MQ=255)
                              ggATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCgggg                                               >  1:3023611/1‑76 (MQ=255)
                               gATTTTTTATATCTGCGCTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTGCGTGACCGCTCGGGGaa                                             <  1:1880510/76‑1 (MQ=255)
                                 ttttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAg                                            >  2:4201247/1‑76 (MQ=255)
                                 ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                                            >  2:5156644/1‑76 (MQ=255)
                                  tttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAgg                                           <  2:2137643/76‑1 (MQ=255)
                                  tttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAgg                                           >  1:158347/1‑76 (MQ=255)
                                  tttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAgg                                           <  2:5404650/76‑1 (MQ=255)
                                   ttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGc                                          >  1:213810/1‑76 (MQ=255)
                                     ttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaa                                        >  1:1275581/1‑76 (MQ=255)
                                     ttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaa                                        >  1:4741314/1‑76 (MQ=255)
                                      tataTCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaa                                       >  2:1913209/1‑76 (MQ=255)
                                      tataTCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa                                       >  1:4151355/1‑76 (MQ=255)
                                        tatCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaaat                                     <  1:1972907/76‑2 (MQ=255)
                                        tatCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaaat                                     >  2:1800543/1‑75 (MQ=255)
                                        tatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAAc                                     >  1:5490258/1‑76 (MQ=255)
                                             gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGGAGGCAAAATTGcct                                >  1:4654623/1‑74 (MQ=255)
                                             gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGcct                                >  2:2981242/1‑74 (MQ=255)
                                             gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGcct                                >  1:4416128/1‑74 (MQ=255)
                                             gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGcct                                >  1:5107487/1‑74 (MQ=255)
                                             gCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGTGTGACGGCTCGGGGaa                                             <  1:2230451/63‑1 (MQ=255)
                                              cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTc                               >  2:627710/1‑76 (MQ=255)
                                              cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTc                               >  1:3071072/1‑76 (MQ=255)
                                                ttCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTctg                             >  1:4268078/1‑74 (MQ=255)
                                                ttCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCag                             <  2:3217313/76‑1 (MQ=255)
                                                 tCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga                            <  1:3925238/76‑1 (MQ=255)
                                                  ccGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag                           >  1:5851352/1‑76 (MQ=255)
                                                   cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAgaga                          >  2:5439902/1‑76 (MQ=255)
                                                   cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAgaga                          >  1:989622/1‑76 (MQ=255)
                                                    gCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaa                         <  1:1864376/76‑1 (MQ=255)
                                                     cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCa                              <  2:2195661/70‑1 (MQ=255)
                                                     cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCa                              >  1:2195451/1‑70 (MQ=255)
                                                     cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaaa                        <  2:3308879/76‑1 (MQ=255)
                                                      tAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAg                       <  1:4117376/76‑1 (MQ=255)
                                                       aaaagGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGc                      >  2:2243516/1‑76 (MQ=255)
                                                       aaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGCGAAAGc                      >  2:1293766/1‑76 (MQ=255)
                                                       aaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGc                      <  1:2810434/76‑1 (MQ=255)
                                                         aaggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCAt                    >  2:4502007/4‑76 (MQ=255)
                                                          aggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCAtt                   <  2:2872922/74‑1 (MQ=255)
                                                          aggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCAtt                   <  2:5207323/74‑1 (MQ=255)
                                                          aggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGACTCTGGGAAAGCAtt                   >  1:4036742/3‑76 (MQ=255)
                                                          aaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaag                        >  2:5611579/1‑70 (MQ=255)
                                                          aaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaag                        <  1:5611581/71‑2 (MQ=255)
                                                           aGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGGAAAACTGCATCAGAGAAAGCATTg                  >  1:5710626/1‑76 (MQ=255)
                                                            gTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTgg                 >  1:3798159/1‑75 (MQ=255)
                                                             tGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCg                >  2:803533/1‑76 (MQ=255)
                                                              gCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCgc               <  2:1479254/76‑2 (MQ=255)
                                                               cAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTg              <  2:3026393/76‑1 (MQ=255)
                                                                aaaTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCgcgg             <  2:1828225/76‑4 (MQ=255)
                                                                aaaTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCgcgg             <  2:4325397/76‑4 (MQ=255)
                                                                aaaTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCgg             >  1:4436147/1‑76 (MQ=255)
                                                                aaaTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgg             <  1:5439903/76‑1 (MQ=255)
                                                                aaaTGAGCAGGACGTTGCAGCGTTTTGCGTGCCCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCgcgg             >  1:5310474/1‑73 (MQ=255)
                                                                 aaTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTggg            >  2:3581000/1‑76 (MQ=255)
                                                                 aaTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTggg            <  2:1548552/76‑1 (MQ=255)
                                                                  aTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgggg           >  2:2566356/1‑76 (MQ=255)
                                                                  aTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgggg           >  1:2402860/1‑76 (MQ=255)
                                                                   tGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTggggc          >  2:5539625/1‑75 (MQ=255)
                                                                    gATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTc         >  2:3052090/1‑76 (MQ=255)
                                                                    gATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTc         >  2:4169955/1‑76 (MQ=255)
                                                                     ctCAGGGCGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTcc        <  2:334345/75‑2 (MQ=255)
                                                                      tCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTg       >  1:277072/1‑76 (MQ=255)
                                                                      tCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTg       <  1:3820099/76‑1 (MQ=255)
                                                                        aGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTggg     >  2:134959/1‑75 (MQ=255)
                                                                        aGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGgc     >  1:4021653/1‑76 (MQ=255)
                                                                        aGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGGGGGGTCTGgc     >  2:2310317/1‑76 (MQ=255)
                                                                        aGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAAATGCATCAGAGAAAGCATTGCGTGGGGGCTGgc     >  1:1075896/1‑76 (MQ=255)
                                                                          gttgttGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGgcgc   <  1:2617809/76‑1 (MQ=255)
                                                                           ttgttgCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCt  >  1:3056806/1‑76 (MQ=255)
                                                                           |                                                                           
CCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCT  >  NC_000913/3192550‑3192700

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: