Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,080,811 G→A 6.7% A151V (GCG→GTG)  sbmC ← DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,080,8110GA 6.7% 191.7 / 4.2 60A151V (GCG→GTG) sbmCDNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base A (2/2);  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 = 1.55e-01

TCCGACGGACAAACCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATC  >  NC_000913/2080663‑2080953
                                                                                                                                                    |                                                                                                                                              
tCCGACGGACAAACCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCg                                                                                                                                                <  2:337089/149‑1 (MQ=255)
 ccGACGGACAAACCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGc                                                                                                                                               <  2:102204/149‑1 (MQ=255)
         cAAACCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAacataca                                                                                                                                       <  2:367230/149‑1 (MQ=255)
         cAAACCGAGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGGATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAacataca                                                                                                                                       <  3:302135/149‑1 (MQ=255)
          aaaCCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAacatacat                                                                                                                                      >  3:299951/1‑149 (MQ=255)
           aaCCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACAtt                                                                                                                                     >  4:567410/1‑149 (MQ=255)
             ccGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTc                                                                                                                                   <  4:185501/149‑1 (MQ=255)
               gggAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGAGATGTTTTGGCTGCACCGCAACATACACTTCGa                                                                                                                                 <  3:437843/149‑1 (MQ=255)
                      gtgtAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATcc                                                                                                                            <  2:487171/147‑1 (MQ=255)
                      gtgtAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATcc                                                                                                                            >  1:487171/1‑147 (MQ=255)
                       tgtAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAg                                                                                                                         >  3:630157/1‑149 (MQ=255)
                                   gCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTcc                                                                                                             <  3:212044/149‑1 (MQ=255)
                                     aGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCcgc                                                                                                           <  2:420537/149‑1 (MQ=255)
                                       gAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCgcgc                                                                                                         >  4:142090/1‑149 (MQ=255)
                                               aaaTTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGTACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTc                                                                                                 <  4:557310/149‑1 (MQ=255)
                                               aaaTTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTc                                                                                                 <  3:312727/149‑1 (MQ=255)
                                                       aGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAAc                                                                                         <  1:432350/149‑1 (MQ=255)
                                                         gcgcATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAAc                                                                                         <  1:133030/147‑1 (MQ=255)
                                                         gcgcATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAAc                                                                                         >  2:133030/1‑147 (MQ=255)
                                                                ttACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCa                                                                                >  2:599863/1‑149 (MQ=255)
                                                                  aCTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATGCATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAgg                                                                              <  3:318818/149‑1 (MQ=255)
                                                                    tctTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCaaa                                                                                              >  3:425801/1‑131 (MQ=255)
                                                                    tctTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCaaa                                                                                              <  4:425801/131‑1 (MQ=255)
                                                                      ttttCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCtt                                                                          >  2:455814/1‑149 (MQ=255)
                                                                        ttCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTg                                                                        <  3:142090/149‑1 (MQ=255)
                                                                         tCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGa                                                                                                                                 <  4:657941/91‑1 (MQ=255)
                                                                         tCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGa                                                                                                                                 >  3:657941/1‑91 (MQ=255)
                                                                           gCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCa                                                                                <  3:39254/138‑1 (MQ=255)
                                                                           gCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCa                                                                                >  4:39254/1‑138 (MQ=255)
                                                                             tttCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAAc                                                                   >  4:326807/1‑149 (MQ=255)
                                                                             tttCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAAc                                                                   >  4:13306/1‑149 (MQ=255)
                                                                             tttCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAAc                                                                   >  4:186182/1‑149 (MQ=255)
                                                                              ttCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACa                                                                  <  4:299951/149‑1 (MQ=255)
                                                                                cAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCACAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACAtt                                                                >  2:387190/1‑149 (MQ=255)
                                                                                  aaaaaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTc                                                              >  1:499777/1‑149 (MQ=255)
                                                                                  aaaaaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTc                                                              <  2:520376/149‑1 (MQ=255)
                                                                                   aaaaaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCa                                                             >  3:97102/1‑149 (MQ=255)
                                                                                    aaaaaGAGTCGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTAACCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTATCCTTTctt                                                            >  2:180547/1‑147 (MQ=255)
                                                                                      aaaGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATaa                                                          <  4:577358/149‑1 (MQ=255)
                                                                                         gagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCa                                                                                                                          >  3:57779/1‑82 (MQ=255)
                                                                                         gagTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCa                                                                                                                          <  4:57779/82‑1 (MQ=255)
                                                                                          agTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCAc                                                      <  4:97102/149‑1 (MQ=255)
                                                                                           gTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACt                                                     <  3:567410/149‑1 (MQ=255)
                                                                                                 aTCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCaaa                                                                                              >  1:181922/1‑102 (MQ=255)
                                                                                                 aTCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCaaa                                                                                              <  2:181922/102‑1 (MQ=255)
                                                                                                      gTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAaga                                          >  4:159738/1‑149 (MQ=255)
                                                                                                      gTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAaga                                          <  1:447727/149‑1 (MQ=255)
                                                                                                      gTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACCGTCCTGCAaga                                          <  2:499777/149‑1 (MQ=255)
                                                                                                                    cTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAAc                            >  3:136884/1‑149 (MQ=255)
                                                                                                                       agaTGAATTAGTGATGTTTTGGCTGCACCACAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGa                         <  1:387190/149‑1 (MQ=255)
                                                                                                                          tGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATAc                      >  3:629833/1‑149 (MQ=255)
                                                                                                                              ttAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAgg                  <  3:326807/149‑1 (MQ=255)
                                                                                                                                 gTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGtt               <  1:599863/149‑1 (MQ=255)
                                                                                                                                  tGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGttt              >  4:180062/1‑149 (MQ=255)
                                                                                                                                   gATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGTTTa             <  4:630157/149‑1 (MQ=255)
                                                                                                                                       ttttGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaaa                                >  2:244903/1‑126 (MQ=255)
                                                                                                                                       ttttGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTaaa                                <  1:244903/126‑1 (MQ=255)
                                                                                                                                       ttttGGCTGCACCACAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAg                   >  2:187114/1‑139 (MQ=255)
                                                                                                                                       ttttGGCTGCACCACAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAg                   <  1:187114/139‑1 (MQ=255)
                                                                                                                                            gCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAatca    >  3:130731/1‑149 (MQ=255)
                                                                                                                                              tGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAatcatc  <  3:186182/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TCCGACGGACAAACCGGGAACAGTGTAATTTCCCGGCAGGAGCGAAAAAATTGAAAGGCGCATCTTACTCTTTTCGCTTTCAAAAAAAAGAGTGGTCATCGCGTTAACACACCGCCCTGAGATGAATTAGTGATGTTTTGGCTGCACCGCAACATACATTTCGATATCCCAGTACCCATCTTCCGCGCCATTGTTCAAATAAACCTCGAAGCAGGGCTTTGGTAACATTTCATAAGCACTGTCCTGCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATC  >  NC_000913/2080663‑2080953

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