Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,192,6260TC45.0% 86.4 / 194.5 181intergenic (+85/+335)b3051/b3052putative membrane protein/putative kinase
Reads supporting (aligned to +/- strand):  ref base T (56/43);  new base C (50/31);  total (106/75)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 191 total aligned reads displayed.

CAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTC  >  NC_000913/3192551‑3192701
                                                                           |                                                                           
cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt                                                                             <  2:1994045/76‑1 (MQ=255)
  ggAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt                                                                           >  1:4006153/1‑76 (MQ=255)
      atATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCa                                                                       >  2:2269277‑M1/3‑76 (MQ=255)
       aatGCCTGGCGTAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcag                                                                      <  2:1616738/76‑1 (MQ=255)
        atCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTgcagc                                                                     <  2:4625058/76‑1 (MQ=255)
        agcccTGGCCCAAAACGCCGGGGATTTTTTATAGCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcagc                                                                     <  2:5014452/74‑1 (MQ=255)
            cTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAttt                                                                 >  1:5850372/1‑76 (MQ=255)
             tGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGTTCAGGTTGTTGCAGCAtttt                                                                >  2:5086806/1‑76 (MQ=255)
               gCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGc                                                              >  1:5703546/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)
                   aaaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGa                                                          <  1:4329506/76‑1 (MQ=255)
                   aaaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGa                                                          >  1:5404651/1‑76 (MQ=255)
                    aaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGAc                                                         <  2:4460866/76‑1 (MQ=255)
                     aaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAcc                                                        >  1:3802625/1‑75 (MQ=255)
                     aaaGCCGGGGATTTTTTATATCTGCGTTCCGCGAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACg                                                        <  2:232189/76‑1 (MQ=255)
                     aaaGCCGGGGATTTTTTATATCTGCGCTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAcc                                                        >  1:3226906/1‑75 (MQ=255)
                      aaGCCGGGGATTGTTTATATCCGCGTTCCGATAAAAGGTGCAAATGCGCAGGACGCTGCCGCGTTTTGCGCGAccc                                                       >  1:1570815/1‑74 (MQ=255)
                       aGCCGGGGATTTTTTCTATCTGCGTTCCGCTAAAAAGTGCCAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                                                      <  2:5721212/76‑1 (MQ=255)
                         ccGTGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTc                                                    <  2:2659483/76‑1 (MQ=255)
                         ccGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTc                                                    <  1:1013703/76‑1 (MQ=255)
                          cGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCg                                                   <  1:5263050/76‑1 (MQ=255)
                           ggggATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCgg                                                  <  2:5797169/76‑1 (MQ=255)
                           ggggATGTTTTATATCTGCGTTCAGCTAAAAGGTGCAAACGCTCAGGACGGTGCACCGTTTTGCGCGACCGCTCgg                                                  >  2:674375/1‑76 (MQ=255)
                             gtaTTTTGTATATCTGCGGTACGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCgggn                                                <  1:3136036/74‑2 (MQ=255)
                             ggATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCgggg                                                >  1:3023611/1‑76 (MQ=255)
                              gATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGa                                               <  1:4034011/76‑1 (MQ=255)
                                ttttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAg                                             >  2:4201247/1‑76 (MQ=255)
                                ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                                             <  1:4691651/76‑1 (MQ=255)
                                 tttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAgg                                            <  2:2137643/76‑1 (MQ=255)
                                 tttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAgg                                            <  2:5404650/76‑1 (MQ=255)
                                  ttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGc                                           <  2:3023611/76‑1 (MQ=255)
                                  ttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGc                                           <  1:2029542/76‑1 (MQ=255)
                                     tataTCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaa                                        >  2:1913209/1‑76 (MQ=255)
                                     tataTCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa                                        >  1:4151355/1‑76 (MQ=255)
                                      ttatCTGCGTGCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaaa                                       <  2:290930/75‑1 (MQ=255)
                                      atatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaaa                                       >  1:382203/1‑76 (MQ=255)
                                       tatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAAc                                      >  1:4841754/1‑76 (MQ=255)
                                       tatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAAc                                      >  1:5490258/1‑76 (MQ=255)
                                       tatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAAc                                      >  1:3234336/1‑76 (MQ=255)
                                        atCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaaatt                                     >  2:389483/1‑74 (MQ=255)
                                        atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACt                                     >  1:1620923/1‑76 (MQ=255)
                                         tCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaaattg                                    <  1:5784742/76‑4 (MQ=255)
                                         tCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTg                                    >  1:443330/1‑76 (MQ=255)
                                         tCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTg                                    >  2:2801182/1‑76 (MQ=255)
                                         tCTGCGGTCCTCTAACAGCTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCaaaattg                                    <  1:5092275/76‑4 (MQ=255)
                                            gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGcct                                 <  1:1800543/76‑3 (MQ=255)
                                            gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGcct                                 >  1:4416128/1‑74 (MQ=255)
                                            gCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGTGTGACGGCTCGGGGaa                                              <  1:2230451/63‑1 (MQ=255)
                                             cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGcctc                                <  2:2619494/76‑4 (MQ=255)
                                             cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTc                                >  1:3071072/1‑76 (MQ=255)
                                             cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTc                                >  2:627710/1‑76 (MQ=255)
                                             cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTc                                >  1:2194669/1‑76 (MQ=255)
                                               ttCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTctg                              >  1:4268078/1‑74 (MQ=255)
                                               ttCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCag                              <  2:3217313/76‑1 (MQ=255)
                                                tCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTctgg                             >  2:1822569/1‑73 (MQ=255)
                                                tCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga                             <  1:1877734/76‑1 (MQ=255)
                                                  cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAgaga                           >  2:5439902/1‑76 (MQ=255)
                                                  cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAAGCAAAACTGCATCAgaga                           >  1:813535/1‑76 (MQ=255)
                                                    cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCa                               >  1:2195661/1‑70 (MQ=255)
                                                    cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCa                               <  2:2195661/70‑1 (MQ=255)
                                                    cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaaa                         <  2:3308879/76‑1 (MQ=255)
                                                     tAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAg                        <  1:4117376/76‑1 (MQ=255)
                                                      gaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGc                       <  2:5198823/75‑1 (MQ=255)
                                                      aaaagGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGc                       >  2:2243516/1‑76 (MQ=255)
                                                      aaaagGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGc                       <  1:5180345/76‑1 (MQ=255)
                                                      aaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGCGAAAGc                       >  2:1293766/1‑76 (MQ=255)
                                                      aaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGc                       <  1:2810434/76‑1 (MQ=255)
                                                        aaggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCAt                     <  1:1913209/73‑1 (MQ=255)
                                                        aaggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCAt                     >  2:4502007/4‑76 (MQ=255)
                                                         aggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCAtt                    <  2:5207323/74‑1 (MQ=255)
                                                         aaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaag                         >  2:5611579/1‑70 (MQ=255)
                                                         aaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaag                         <  1:5611581/71‑2 (MQ=255)
                                                          ggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTg                   >  1:593975/2‑76 (MQ=255)
                                                          aGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTg                   <  1:5037075/76‑1 (MQ=255)
                                                            tGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCg                 >  2:803533/1‑76 (MQ=255)
                                                             gCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCgc                <  2:1479254/76‑2 (MQ=255)
                                                               aaaTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCgcgg              <  2:4325397/76‑4 (MQ=255)
                                                               aaaTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCgg              >  1:921367/1‑76 (MQ=255)
                                                               aaaTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgg              >  1:2098631/1‑76 (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)
                                                                  tgctCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGt           <  2:5069328/73‑1 (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)
                                                                    ctCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTcc         <  2:1692173/75‑2 (MQ=255)
                                                                    aTCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCt         <  1:4600117/76‑1 (MQ=255)
                                                                     tCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTg        >  1:277072/1‑76 (MQ=255)
                                                                       aGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGCAGGCCACACTGCCTCAGAGAAAGCATTTCGGGGGGTCTGgc      >  1:4061383/1‑76 (MQ=255)
                                                                       aGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTggg      >  2:134959/1‑75 (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)
                                                                           tgttgCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGGCTGGCGCtc  >  1:1225225/1‑76 (MQ=255)
                                                                           |                                                                           
CAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTC  >  NC_000913/3192551‑3192701

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