Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,371,827 C→T 7.0% C37Y (TGT→TAT)  carB ← carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,371,8270CT 7.0% 169.5 / 2.4 57C37Y (TGT→TAT) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (28/25);  new base T (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GCGCTCTTTCTCGATGATCTTGGCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAAT  >  NC_002947/5371682‑5371974
                                                                                                                                                 |                                                                                                                                                   
gcgcTCTTTCTCGATGATCTTGGCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAgg                                                                                                                                                  >  3:168922/1‑149 (MQ=255)
          tcGATGATCTTGGCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGcc                                                                                                                                        >  8:31479/1‑149 (MQ=255)
               gatCTTGGCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGag                                                                                                                                    >  3:49736/1‑148 (MQ=255)
                 tCTTGGCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGGCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCGTTACCGGCCTGGGCGCCGGagtag                                                                                                                                 >  4:173266/1‑149 (MQ=255)
                     ggCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTc                                                                                                                               <  4:30797/147‑1 (MQ=255)
                         aCCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTATAGGCCTGGGCGCCGGAGTAGTCGAATTc                                                                                                                         >  8:158784/1‑149 (MQ=255)
                         aCCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTATAGGCCTGGGCGCCGGAGTAGTCGAATTc                                                                                                                         >  2:258212/1‑149 (MQ=255)
                         aCCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTc                                                                                                                         <  7:195785/149‑1 (MQ=255)
                            gATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCa                                                                                                                      <  6:57913/149‑1 (MQ=255)
                                  cACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCAGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCGTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTg                                                                                                                >  8:39743/1‑149 (MQ=255)
                                    cTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTATAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGc                                                                                                              <  1:258212/149‑1 (MQ=255)
                                    cTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTATAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGc                                                                                                              <  7:158784/149‑1 (MQ=255)
                                        aTCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGAt                                                                                                          >  3:226005/1‑149 (MQ=255)
                                         tCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATc                                                                                                         <  1:210990/149‑1 (MQ=255)
                                            gCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGa                                                                                                                             <  1:234226/126‑1 (MQ=255)
                                            gCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGa                                                                                                                             >  2:234226/1‑126 (MQ=255)
                                            gCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAAtt                                                                                                                          >  8:58900/1‑129 (MQ=255)
                                            gCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAAtt                                                                                                                          <  7:58900/129‑1 (MQ=255)
                                              tCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGAt                                                                                                    >  8:97645/1‑149 (MQ=255)
                                                 aTGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAATCCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCgg                                                                                                 <  1:171940/149‑1 (MQ=255)
                                                 aTGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGAt                                                                                                          <  3:152812/140‑1 (MQ=255)
                                                 aTGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGAt                                                                                                          >  4:152812/1‑140 (MQ=255)
                                                 aTGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCgg                                                                                                 <  7:31479/149‑1 (MQ=255)
                                                    tAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGc                                                                                               >  6:97903/1‑148 (MQ=255)
                                                     aGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTc                                                                                                                               >  3:29033/1‑115 (MQ=255)
                                                     aGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTc                                                                                                                               <  4:29033/115‑1 (MQ=255)
                                                     aGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCa                                                                                             >  4:154641/1‑149 (MQ=255)
                                                      ggTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAg                                                                                            >  4:107362/1‑149 (MQ=255)
                                                       gTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAgc                                                                                           <  4:49736/149‑1 (MQ=255)
                                                            gTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGcc                                                                                        >  1:11181/1‑147 (MQ=255)
                                                            gTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCga                                                                                      <  7:134159/149‑1 (MQ=255)
                                                             tCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGcc                                                                                        >  8:131361/1‑146 (MQ=255)
                                                               aGCCATGGCCGGGTCGGTCATGATGGTGTCTGGGTTGGAGTTAACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATAACGATAGGGCCAGCGCCGATaa                                                                                   >  7:135119/1‑149 (MQ=255)
                                                                 ccATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGCCGAATTCTCAGACCTGGCCGATCACGATCGGGCCAGCGCCGAGAAt                                                                                  >  6:220149/1‑148 (MQ=255)
                                                                      gCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGa                                                                             <  7:97645/148‑1 (MQ=255)
                                                                      gCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGa                                                                             <  4:168922/148‑1 (MQ=255)
                                                                               gtcATGATGGAGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGGCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGAt                                                                   <  2:225457/149‑1 (MQ=255)
                                                                                tcATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATg                                                                  >  1:188084/1‑149 (MQ=255)
                                                                                tcATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATg                                                                  >  7:217906/1‑149 (MQ=255)
                                                                                    gatgGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGagt                                                                                                                                   <  8:136538/80‑1 (MQ=255)
                                                                                    gatgGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGagt                                                                                                                                   >  7:136538/1‑80 (MQ=255)
                                                                                           gCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGttt                                                       >  1:230367/1‑149 (MQ=255)
                                                                                               ggTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTgg                                                    >  4:108710/1‑148 (MQ=255)
                                                                                                   ggAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCAtg                                                <  3:108710/148‑1 (MQ=255)
                                                                                                     aGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCAtggtg                                             <  2:188084/149‑1 (MQ=255)
                                                                                                     aGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCAtggtg                                             >  3:81996/1‑149 (MQ=255)
                                                                                                       ttCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATggtggt                                           <  2:230367/149‑1 (MQ=255)
                                                                                                             aGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTc                                      >  3:132967/1‑148 (MQ=255)
                                                                                                               gATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCAAGGGCGCCGGAGTAGTCGAATTCGAAGGCCTGTCCGATCACGATCGGGCCAGCGCCGAGAATCAGGAGGCTTTAGAAGTCAGTACGTTTTGGCATGGTGGTCACTCaa                                    >  4:260685/1‑148 (MQ=255)
                                                                                                                    cgcgGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCgg                              <  4:132967/149‑1 (MQ=255)
                                                                                                                     gcgGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCggg                             >  6:18475/1‑149 (MQ=255)
                                                                                                                                     cgcAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCg             <  3:107362/149‑1 (MQ=255)
                                                                                                                                         gggCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCttt         <  4:81996/149‑1 (MQ=255)
                                                                                                                                          ggCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTg                                                     >  1:116677/1‑104 (MQ=255)
                                                                                                                                          ggCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTg                                                     <  2:116677/104‑1 (MQ=255)
                                                                                                                                            cTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGaa      <  4:226005/149‑1 (MQ=255)
                                                                                                                                            cTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGaa      <  5:18475/149‑1 (MQ=255)
                                                                                                                                                 cAGGCCTTGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAAt  <  3:154641/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGCTCTTTCTCGATGATCTTGGCCACCGATTGCCACTTGATCGGCTCGATGTAGGTGGCGTCAGCCATGGCCGGGTCGGTCATGATGGTGGCTGGGTTGGAGTTCACCAGGATGACGCGGAAACCTTCCTCGCGCAGGGCTTTACAGGCCTGGGCGCCGGAGTAGTCGAATTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAAT  >  NC_002947/5371682‑5371974

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: ‑