Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,785,970 C→G 6.8% P325R (CCG→CGG)  PP_5068 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,785,9700CG 6.8% 177.8 / 3.9 59P325R (CCG→CGG) PP_5068hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

ATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGT  >  NC_002947/5785825‑5786118
                                                                                                                                                 |                                                                                                                                                       
aTCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGc                                                                                                                                                       <  6:34762/148‑1 (MQ=255)
    gCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTAc                                                                                                                                                  >  2:4891/1‑149 (MQ=255)
     cAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTg                                                                                                                                                     <  8:269541/145‑1 (MQ=255)
     cAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTa                                                                                                                                                   <  2:30499/147‑1 (MQ=255)
         cACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGc                                                                                                                                             <  6:120726/149‑1 (MQ=255)
              cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaac                                                                                                                                        <  4:195591/149‑3 (MQ=255)
                                cAGCTACGCCAACCCTGTGCCCATCGCCCACTGTCACCTCAGCGCGCTGGCCCAGGAACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgcgc                                                                                                                       <  5:293472/148‑1 (MQ=255)
                                 aGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCg                                                                                                                     >  3:134368/1‑149 (MQ=255)
                                    tACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                                  >  8:108099/1‑149 (MQ=255)
                                    tACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                                  >  1:159101/1‑149 (MQ=255)
                                      cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCg                                                                                                                                                        <  4:156693/109‑1 (MQ=255)
                                      cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCg                                                                                                                                                        >  3:156693/1‑109 (MQ=255)
                                      cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcgc                                                                                                                 >  4:9720/1‑148 (MQ=255)
                                       gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTTGGCCTGTTCCTGCCGCTGTAACAGGCCCATTACCTGGACCTGATGCGACGgcgca                                                                                                                >  6:244348/1‑147 (MQ=255)
                                       gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGa                                                                                                                                 >  2:126027/1‑131 (MQ=255)
                                       gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGa                                                                                                                                 <  1:126027/131‑1 (MQ=255)
                                         cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaa                                                                                                                                         <  6:263299/121‑2 (MQ=38)
                                         cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaa                                                                                                                                         >  5:263299/1‑120 (MQ=38)
                                             cAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCt                                                                                                         <  4:134368/149‑1 (MQ=255)
                                              aGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                                    >  2:165251/1‑121 (MQ=255)
                                              aGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                                    <  1:165251/121‑1 (MQ=255)
                                               ggTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGa                                                                                                       <  1:4891/149‑1 (MQ=255)
                                                     cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                        <  1:155112/142‑1 (MQ=255)
                                                     cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                        >  2:155112/1‑142 (MQ=255)
                                                       tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTa                                                                                                                                       >  1:148964/1‑109 (MQ=255)
                                                       tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTa                                                                                                                                       <  2:148964/109‑1 (MQ=255)
                                                       tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                                >  7:115976/1‑148 (MQ=255)
                                                       tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                                <  8:115976/148‑1 (MQ=255)
                                                        cGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCAGCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                           <  4:67698/136‑1 (MQ=255)
                                                        cGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                           >  3:67698/1‑136 (MQ=255)
                                                          cccACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGaa                                                                                            >  1:344833/1‑149 (MQ=255)
                                                               tGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgc                                                                                                                         <  3:98220/115‑1 (MQ=255)
                                                               tGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgc                                                                                                                         >  4:98220/1‑115 (MQ=255)
                                                                    ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTAcc                                                                                                                                     >  1:139433/1‑98 (MQ=255)
                                                                    ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTAcc                                                                                                                                     <  2:139433/98‑1 (MQ=255)
                                                                    ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACa                                                                               >  3:344653/1‑149 (MQ=255)
                                                                     cTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GAtgatga                                                                                 >  2:217648/1‑146 (MQ=255)
                                                                     cTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GAtgatga                                                                                 <  1:217648/146‑1 (MQ=255)
                                                                       cAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATgg                                                                            >  5:280011/1‑149 (MQ=255)
                                                                          cgcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGc                                                                         <  7:108099/149‑1 (MQ=255)
                                                                            cgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTg                                                                       <  6:280011/149‑1 (MQ=255)
                                                                                 gCCCCGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGa                                                                  <  4:346855/149‑1 (MQ=255)
                                                                                   ccAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAgc                                                                <  3:9720/149‑1 (MQ=255)
                                                                                     aGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                                >  1:222101/1‑86 (MQ=255)
                                                                                     aGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                                <  2:222101/86‑1 (MQ=255)
                                                                                        cACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCctg                                                           <  4:344653/149‑1 (MQ=255)
                                                                                               caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                        <  8:201756/100‑1 (MQ=255)
                                                                                               caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                        >  7:201756/1‑100 (MQ=255)
                                                                                               caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                    <  2:159101/149‑1 (MQ=255)
                                                                                               caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGCGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCAGACAACCGCCGAA‑‑GATGATGACATGGTGCTTGTCGTGCGGCTGCTGCAGt                                                    >  6:203901/1‑149 (MQ=255)
                                                                                                 ccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCTAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgc                                                  <  5:244348/149‑1 (MQ=255)
                                                                                                    tGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgcatg                                               >  3:92324/1‑149 (MQ=255)
                                                                                                         gAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAcg                                          >  6:89036/1‑149 (MQ=255)
                                                                                                               gAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAATAAGTGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGAgg                                       >  1:294156/1‑149 (MQ=255)
                                                                                                                     cccatcgaagCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTg                                <  1:305358/141‑1 (MQ=255)
                                                                                                                     cTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGa                              >  7:219081/1‑149 (MQ=255)
                                                                                                                            aggtaatGGGCCTGGTACAGCGGCAGGAACAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGc                       >  7:145827/7‑149 (MQ=38)
                                                                                                                               cGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                    <  5:89036/149‑1 (MQ=255)
                                                                                                                                            ccTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAAc       >  5:107847/1‑149 (MQ=255)
                                                                                                                                                 cGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGt  <  8:219081/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
ATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGT  >  NC_002947/5785825‑5786118

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