Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 584,121 G→T 77.2% intergenic (+205/‑67) rplA → / → rplJ ribosomal protein L1/ribosomal protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730584,1210GT77.2% 95.8 / 33.8 61intergenic (+205/‑67)rplA/rplJribosomal protein L1/ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base T (22/25);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

AAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTT  >  CP000730/583982‑584256
                                                                                                                                            |                                                                                                                                        
aaGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGt                                                                                                                                         <  2:525482/141‑1 (MQ=255)
aaGTAATTTCCAAGTTATATTACGTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGt                                                                                                                                         <  2:13027/141‑1 (MQ=255)
  gTAATTTCTAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTgg                                                                                                                                       <  1:182752/141‑1 (MQ=255)
       ttCCAAGTTTTATTTCTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGa                                                                                                                                  <  2:256261/141‑1 (MQ=255)
       ttCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGa                                                                                                                                  <  2:520260/141‑1 (MQ=255)
        tCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGaa                                                                                                                                 >  2:335975/1‑141 (MQ=255)
         ccAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGAAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGaaa                                                                                                                                >  2:155614/1‑141 (MQ=255)
          cAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAg                                                                                                                               <  2:184958/141‑1 (MQ=255)
             gTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGc                                                                                                                            >  1:328445/1‑141 (MQ=255)
             gTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGc                                                                                                                            >  1:99504/1‑141 (MQ=255)
               tataTTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGCTGCTGTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCtt                                                                                                                         >  1:568342/1‑141 (MQ=255)
                atatTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTCGGTAGAAAGTGCtct                                                                                                                         <  2:559250/141‑3 (MQ=255)
                          ttGTGATTATTTTACCTAAGCCAGTAGGAGTTATTTATAACTTAAAATTTATCCTTCCGAGGCTAAAATTGACTTGCACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAt                                                                                                              <  1:1671/141‑1 (MQ=255)
                          ttGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAt                                                                                                              <  1:503715/141‑1 (MQ=255)
                          ttGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAAt                                                                                                              <  1:1655/141‑1 (MQ=255)
                              aatcatctaaCCTAAAACAATAGGAGTTATATAAAACTAAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATAATCTTTCAAGCACTTTTTGCCTTAGGTAGAAAGTGCTTTTTTTATTAATTTTa                                                                                                          <  2:209965/132‑1 (MQ=255)
                                     ttACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATTATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAg                                                                                                   <  2:589690/141‑1 (MQ=255)
                                     ttACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAg                                                                                                   >  1:411319/1‑141 (MQ=255)
                                     ttACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAg                                                                                                   >  1:461521/1‑141 (MQ=255)
                                      tACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGc                                                                                                  >  1:96787/1‑141 (MQ=255)
                                      tACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTACTTTTTTTATTAATTTTAAAAAAAGc                                                                                                  >  1:350734/1‑141 (MQ=255)
                                                   aGGAGTTATTTATAACCTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTaa                                                                                     >  1:304174/1‑141 (MQ=255)
                                                    ggAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTaaa                                                                                    >  1:358903/1‑141 (MQ=255)
                                                    ggAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTaaa                                                                                    <  1:375747/141‑1 (MQ=255)
                                                      aGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTACTTTTAAAAAAGGCACCAAAAATTTAAATg                                                                                  >  1:44494/1‑141 (MQ=255)
                                                             tataACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGgtgt                                                                           <  1:392944/141‑1 (MQ=255)
                                                              ataACTTAAAATTTCTCCTGCCGAGGCTAAAATCGTCCTGAACGTGTTCAACTATGATCTTTCAAGGACAGGTATGCTGGTCTCGAGTTGTTCTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTc                                                                          <  1:143565/141‑1 (MQ=255)
                                                               taACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCt                                                                         <  1:520076/141‑1 (MQ=255)
                                                                       aaTTTATCCTGCCGAGGCTAAAATTGACCTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCt                                                                 <  1:424269/141‑1 (MQ=255)
                                                                        aTTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCAATGATCTTTCCAGCACTTTTTGCCGTGGGTAGAAAACGCTTTTTTTATTAATGTTAAAAAAATCACGAAAAATATAAATGGCGGGGTCTGAATGTCTg                                                                >  2:198917/1‑141 (MQ=255)
                                                                            aTCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTAt                                                            >  1:513557/1‑141 (MQ=255)
                                                                            aTCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTAt                                                            >  1:432382/1‑141 (MQ=255)
                                                                              ccTGCCGAGGCTAAAATTGACTTGAACGTGATGAACTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTAAAAAATAAAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCa                                                          <  1:89074/141‑1 (MQ=255)
                                                                                tGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCAtt                                                        >  2:26514/1‑141 (MQ=255)
                                                                                  ccGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAACGCAAAAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGa                                                      >  1:567989/1‑141 (MQ=255)
                                                                                   cggggCTTAAATTCTCGTGAACGCGTTGATCTATGGTCGTCCAAGGAGTTTGTGCCTGTTGTAGAAAGTGCTTTTTTTATTAATGTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGaa                                                     <  2:419974/138‑1 (MQ=255)
                                                                                   cGAGGCTAAAATTGACTTGAACGTGATGAGCTATGATCTTTCAATCACATTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTGATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGaa                                                     <  2:96371/141‑1 (MQ=255)
                                                                                     aGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGc                                                   >  1:431516/1‑141 (MQ=255)
                                                                                        cTAAAATTGACTTGAACGTGCTGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTaa                                                >  1:371270/1‑141 (MQ=255)
                                                                                        cTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCGTGATGTCTGCTATCGTTGAAGCTaa                                                >  1:77025/1‑141 (MQ=255)
                                                                                          aaaaTTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTaaaa                                              <  2:411319/141‑1 (MQ=255)
                                                                                                gACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTTTTTTTTTTATAAAGTGCTTTTTTTATTAATTTTAAATAAAGCACCAAAAATTTAAATGGCGGTGTCTGTATGTCTGATATCAATGAAGCTAAAAaacaac                                        >  2:472586/1‑141 (MQ=255)
                                                                                                    tGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGt                                    <  1:565681/141‑1 (MQ=255)
                                                                                                    tGAACGTGATGATCCATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGt                                    <  1:565682/141‑1 (MQ=255)
                                                                                                     gAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTACATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAGAAAACAACTAGtg                                   >  1:57009/1‑140 (MQ=255)
                                                                                                        cGTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGATAAAAAACAACTAGTtgat                                <  2:234076/141‑1 (MQ=255)
                                                                                                         gTGATGATCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTtgatg                               >  2:604236/1‑141 (MQ=255)
                                                                                                          tgatgaTCTATGATCTTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAACAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTtgatga                              <  1:491685/141‑1 (MQ=255)
                                                                                                           tatgatctatgatctTTTACGGGTGTAGTGTCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGGTGaa                             <  2:426804/140‑1 (MQ=255)
                                                                                                           gatgatCTATGATCTTTGAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAATTTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGaa                             <  1:235736/141‑1 (MQ=255)
                                                                                                            atgatctatgatctTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGaaa                            >  2:497581/1‑141 (MQ=255)
                                                                                                                  tatgatctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGAAACAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTg                      >  2:217955/1‑141 (MQ=255)
                                                                                                                      atctTTCAAGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGt                  <  1:335975/141‑1 (MQ=255)
                                                                                                                      atctTTCAAGCACTGTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAACTTTAAATGGAGGTGTCTGAATGTCTGCGATCATGGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGt                  <  1:360098/141‑1 (MQ=255)
                                                                                                                             aaGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCa           <  1:364752/141‑1 (MQ=255)
                                                                                                                             aaGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTTTTATTTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCACTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTCTCa           >  2:230652/1‑141 (MQ=255)
                                                                                                                              aGCACTTTTTGCCTTTGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          <  1:15666/141‑1 (MQ=255)
                                                                                                                              aGCACTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCaa          <  1:540996/141‑1 (MQ=255)
                                                                                                                                  cTTTTTGCCTTGGGTAGAAAGTGCTTTTTTTATTTTTTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTc      >  2:470571/1‑141 (MQ=255)
                                                                                                                                       ttGCCTTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGTTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCTGtt  >  1:321087/1‑141 (MQ=255)
                                                                                                                                       ttGCCGTGGGTAGAAAGTGCTTTTTTTTTTTTTTTTAAAAAAAGCACCAAAAATTTAAATTGAGATGGCTGAATGTCTGCTATCATTGAAGCTAAACAACACGTAGTTGTTTAAACTGCTGAGGTACGATCTAAtttattt  >  2:103668/1‑136 (MQ=255)
                                                                                                                                            |                                                                                                                                        
AAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACTTAAAATTTATCCTGCCGAGGCTAAAATTGACTTGAACGTGATGATCTATGATCTTTCAAGCACTTTTTGCCGTGGGTAGAAAGTGCTTTTTTTATTAATTTTAAAAAAAGCACCAAAAATTTAAATGGAGGTGTCTGAATGTCTGCTATCATTGAAGCTAAAAAACAACTAGTTGATGAAATTGCTGAGGTACTATCAAATTCAGTT  >  CP000730/583982‑584256

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