Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,563,817 A→G 7.1% intergenic (+83/+39) PP_1372 → / ← pitB hypothetical protein/phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,563,8170AG 7.1% 156.5 / 2.6 56intergenic (+83/+39)PP_1372/pitBhypothetical protein/phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/2);  total (30/26)
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.06e-01

GTTGGGGCGGCAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGC  >  NC_002947/1563674‑1563964
                                                                                                                                                  |                                                                                                                                                     
gTTGGGGCGGCAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGgc                                                                                                                                                 >  4:303676/1‑149 (MQ=255)
   ggggCGGCAGTTGGTGCGCGGGTTGATGGGGAGCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGAACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCAATGAAAAAGgcgc                                                                                                                                               >  1:204653/1‑149 (MQ=255)
   ggggCGGCAGTTGGTGCGCGGGTTGATGGGGTACGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAATCCCGCTCCGGCAAAGAACAAGGTGGCGAAGGCATAACGGCTATGAAAAAGgcgc                                                                                                                                               >  1:200043/1‑149 (MQ=255)
      gcggcAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGccc                                                                                                                                             >  7:184931/1‑147 (MQ=255)
            gTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCATCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGAGACCGAGGgcgc                                                                                                                                     >  1:243783/1‑149 (MQ=255)
             ttGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGcc                                                                                                                                    >  8:5377/1‑149 (MQ=255)
              ttgTGCGCGGGTTGATGGGGTCGTTGTTAGGTGGCAAGAAAAGGTATTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGcc                                                                                                                                    <  2:204653/146‑1 (MQ=255)
              tGGTGCGCGGGTTGATGGGGTCGTTTTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTTGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGcc                                                                                                                                    <  2:200043/148‑1 (MQ=255)
                gTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAGAAGGCGCCCTCGGGCGCCttt                                                                                                                                 >  6:7013/1‑148 (MQ=255)
                   cgcgGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCGCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAGAAGGCGCCCTCGGGCGCCtttttc                                                                                                                              <  5:7013/149‑5 (MQ=255)
                            aTGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGgt                                                                                                                     <  4:282954/149‑1 (MQ=255)
                               gggTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGgtgt                                                                                                                   >  1:29338/1‑148 (MQ=255)
                                       ttgttAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCttca                                                                                                                                >  4:220664/1‑129 (MQ=255)
                                       ttgttAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCCCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCttca                                                                                                                                <  3:220664/129‑1 (MQ=255)
                                            aGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAATGAAGGCGccctc                                                                                                                                           >  8:204847/1‑112 (MQ=255)
                                            aGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAATGAAGGCGccctc                                                                                                                                           <  7:204847/114‑3 (MQ=255)
                                            aGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGAGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAAt                                                                                                       >  2:206118/1‑149 (MQ=255)
                                                  cAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCttt                                                                                                 <  4:220915/149‑1 (MQ=255)
                                                    aGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTgg                                                                                               <  8:184931/149‑1 (MQ=255)
                                                          ggTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgcca                                                                                         <  7:5377/149‑1 (MQ=255)
                                                           gTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgccag                                                                                        >  5:323618/1‑149 (MQ=255)
                                                                    cTGCGCCGCTCCATACCCATGGTAGGACAAGCCCGCTCCGTCAAAGAACAAGTTGGCGTAGTCATTACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGATTTGGACGCCAGCCAGAAc                                                                                 <  4:244263/147‑1 (MQ=255)
                                                                     tGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAgg                                                                              >  4:125246/1‑149 (MQ=255)
                                                                      gcgccgcgccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAggc                                                                             >  5:53119/1‑149 (MQ=255)
                                                                        gccgcgccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAggccg                                                                           <  2:29338/149‑1 (MQ=255)
                                                                         ccgcgccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAggccgg                                                                          >  5:294792/1‑149 (MQ=255)
                                                                             gccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGgccgc                                                                      >  5:176515/1‑149 (MQ=255)
                                                                                aTACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAg                                                                   >  4:181406/1‑149 (MQ=255)
                                                                                 tACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAgg                                                                  <  3:303676/149‑1 (MQ=255)
                                                                                    ccTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGccagccag                                                                                    <  2:90060/128‑1 (MQ=255)
                                                                                    ccTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGccagccag                                                                                    >  1:90060/1‑128 (MQ=255)
                                                                                    ccTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGcc                                                                >  5:290263/1‑148 (MQ=255)
                                                                                      ttGATAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCAt                                                             >  5:213590/1‑149 (MQ=255)
                                                                                         gTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCg                                                           <  3:156615/148‑1 (MQ=255)
                                                                                         gTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCg                                                           <  6:294792/148‑1 (MQ=255)
                                                                                          tAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGag                                                         >  8:104344/1‑149 (MQ=255)
                                                                                            gCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggc                                                       <  6:323618/149‑1 (MQ=255)
                                                                                                       cTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCa                                           <  6:176515/149‑1 (MQ=255)
                                                                                                          cGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAg                                          >  1:226066/1‑148 (MQ=255)
                                                                                                          cGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGa                                        >  8:185031/1‑149 (MQ=255)
                                                                                                          cGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGa                                        <  1:206118/149‑1 (MQ=255)
                                                                                                             cAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGAcc                                      >  4:195895/1‑148 (MQ=255)
                                                                                                             cAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGAcc                                      <  3:195895/148‑1 (MQ=255)
                                                                                                              aaaGAACACGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAggccg                                                                           <  4:6961/111‑1 (MQ=255)
                                                                                                              aaaGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCa                                    >  8:150701/1‑149 (MQ=255)
                                                                                                               aacaaCAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCAt                                   <  7:104344/146‑1 (MQ=255)
                                                                                                                  aaCAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGc                                 >  7:267277/1‑148 (MQ=255)
                                                                                                                   aCAAGGTGGCGTAGGCATAACGGCTATGAAACAGGCGCCCGAGGGCTCCTTCATTCATTGCGGGGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCATGGGTCCTGAccctgcct                                <  6:182061/149‑6 (MQ=255)
                                                                                                                   aCAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGc                                 >  6:210898/1‑147 (MQ=255)
                                                                                                                           gCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGAt                        <  2:226066/148‑1 (MQ=255)
                                                                                                                             gTAGGAATAACGGATTTGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGt                     <  6:213590/149‑1 (MQ=255)
                                                                                                                                 gCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGtttt                  >  1:205294/1‑148 (MQ=255)
                                                                                                                                 gCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTc                 <  7:185031/149‑1 (MQ=255)
                                                                                                                                    tAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                >  8:85193/1‑147 (MQ=255)
                                                                                                                                       cGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGAcc                                      <  4:89775/122‑1 (MQ=255)
                                                                                                                                       cGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGAcc                                      >  3:89775/1‑122 (MQ=255)
                                                                                                                                                aaaaGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGc  <  7:150701/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTTGGGGCGGCAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGC  >  NC_002947/1563674‑1563964

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