Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,796,836 G→T 10.0% intergenic (+20/‑112) PP_t43 → / → PP_2449 tRNA‑Val/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,796,8360GT10.0% 170.2 / 9.4 60intergenic (+20/‑112)PP_t43/PP_2449tRNA‑Val/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/30);  new base T (3/3);  total (27/33)
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.27e-01

AAAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCG  >  NC_002947/2796689‑2796981
                                                                                                                                                   |                                                                                                                                                  
aaaaGAGCGTTACAGCTACGCCACGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGa                                                                                                                                                   <  8:262641/149‑1 (MQ=255)
       cTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                            <  4:23684/149‑1 (MQ=255)
        tttACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCg                                                                                                                                           >  7:251808/1‑149 (MQ=255)
           aCAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGcc                                                                                                                                         <  6:133456/148‑1 (MQ=255)
             aGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCtt                                                                                                                                      >  6:317364/1‑149 (MQ=255)
               gTACGCATCGGGCCACTAATATGCGCCCCACACGACAGGCACATAGCTCAGTTGGATAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGTCGATCAGAAATGATCGCCCtttt                                                                                                                                    >  5:266476/1‑149 (MQ=255)
                tACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCttttt                                                                                                                                   <  1:95531/149‑1 (MQ=255)
                     aTCTGGCCACTAATATGCGCACAACACGACAGGCACATAGCTCAGATGGTTAGAGCACCACCATGAAATGGTGGGGGTCGTTGGTTCGAGTCAAATTGTGCCTACCAAATCGAAAAAAGGCGATCAGAAATTATAGAACTTTTTGAAtt                                                                                                                             >  8:114235/1‑149 (MQ=17)
                        gggCCACTAATATGCTCACCACACGACAGGCACATAGCTCAGTTGGTTAGTGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTgc                                                                                                                          <  8:296918/149‑1 (MQ=255)
                            cACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTggtagg                                                                                                                      <  5:97888/149‑6 (MQ=14)
                                ccgatGCGTACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAag                                                                                                                  <  5:324204/146‑1 (MQ=255)
                                         accacACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCgtgt                                                                                                          >  5:129889/1‑148 (MQ=255)
                                            acacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACg                                                                                                      <  1:161947/149‑1 (MQ=255)
                                                gACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                  <  1:158517/149‑1 (MQ=255)
                                                gACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                  >  2:250509/1‑149 (MQ=255)
                                                 aCAGTCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTgg                                                                                                 <  4:213714/149‑1 (MQ=255)
                                                 aCAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTgg                                                                                                 <  3:72195/149‑1 (MQ=255)
                                                         aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATAAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTc                                                                                         >  8:355219/1‑149 (MQ=255)
                                                         aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGCTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCt                                                                                                                                       <  1:169600/104‑1 (MQ=14)
                                                         aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGCTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCt                                                                                                                                       >  2:169600/1‑104 (MQ=14)
                                                           aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCg                                                                                                                     <  3:372338/119‑1 (MQ=255)
                                                           aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCg                                                                                                                     >  4:372338/1‑119 (MQ=255)
                                                             cTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCggg                                                                                     <  8:2419/149‑1 (MQ=255)
                                                                aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGc                                                                                    <  6:45436/147‑1 (MQ=255)
                                                                 gTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGCTCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGACGCGGCGTGTATCGTGTGGTCGCGGTCCGGGCaga                                                                                 >  4:160413/1‑147 (MQ=21)
                                                                        agagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGc                                                                          >  1:86916/1‑149 (MQ=255)
                                                                                 ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa                                                                 >  7:90731/1‑149 (MQ=255)
                                                                                     gACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATc                                                               <  1:250509/147‑1 (MQ=255)
                                                                                      aCATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCg                                                             <  5:317364/148‑1 (MQ=255)
                                                                                      aCATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGc                                                            <  6:224221/149‑1 (MQ=255)
                                                                                       cATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGcc                                                           <  7:355219/149‑1 (MQ=255)
                                                                                           gtggGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCtt                                                       >  5:371825/1‑149 (MQ=255)
                                                                                              ggggTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                    <  1:359167/133‑1 (MQ=255)
                                                                                              ggggTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                    >  2:359167/1‑133 (MQ=255)
                                                                                                    gTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg                                               <  1:150661/148‑1 (MQ=255)
                                                                                                       ggTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCGTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACTTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCTAATCTGGAg                                            <  5:4126/148‑1 (MQ=255)
                                                                                                        gTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGccc                                                                                                                                        <  5:73206/56‑1 (MQ=255)
                                                                                                        gTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGccc                                                                                                                                        >  6:73206/1‑56 (MQ=255)
                                                                                                               tCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACAcgcg                                   >  4:62/1‑149 (MQ=255)
                                                                                                                 cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTg                                 >  1:63843/1‑149 (MQ=255)
                                                                                                                    ttGTGCCTACCAAATCGAAAAAGGTCGATCAGAAATGCTCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTgcgc                              <  6:266476/149‑1 (MQ=255)
                                                                                                                          cTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAg                                            >  5:348857/1‑129 (MQ=255)
                                                                                                                          cTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAg                                            <  6:348857/129‑1 (MQ=255)
                                                                                                                           tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCg                                                                                       <  8:137507/85‑1 (MQ=255)
                                                                                                                           tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCg                                                                                       >  7:137507/1‑85 (MQ=255)
                                                                                                                           tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                       >  7:117034/1‑149 (MQ=255)
                                                                                                                            aCCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTcgcg                                                                                            >  6:51323/1‑79 (MQ=255)
                                                                                                                            aCCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTcgcg                                                                                            <  5:51323/79‑1 (MQ=255)
                                                                                                                             ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                       <  4:169906/147‑1 (MQ=255)
                                                                                                                             ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGc                      <  2:86916/148‑1 (MQ=255)
                                                                                                                                     aaaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt             <  6:371825/149‑1 (MQ=255)
                                                                                                                                      aaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCTGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATg            <  8:251808/149‑1 (MQ=255)
                                                                                                                                         aGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa         >  4:206070/1‑149 (MQ=255)
                                                                                                                                         aGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa         <  6:129889/149‑1 (MQ=255)
                                                                                                                                          gggCGATCATTTCTGATCGCCCTTTTTGCGATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg           >  1:28537/1‑147 (MQ=255)
                                                                                                                                          gggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAg        >  6:245362/1‑149 (MQ=255)
                                                                                                                                             cGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGc     >  6:256182/1‑149 (MQ=255)
                                                                                                                                               aTCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa         <  3:344937/143‑1 (MQ=255)
                                                                                                                                               aTCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa         >  4:344937/1‑143 (MQ=255)
                                                                                                                                                tCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCg  >  5:223883/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCG  >  NC_002947/2796689‑2796981

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: