Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 507,271 A→C 6.7% noncoding (68/84 nt) chiX → sRNA antisense regulator of OM chitoporin ChiP(YbfM), Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913507,2710AC 6.7% 199.8 / 5.5 60noncoding (68/84 nt)chiXsRNA antisense regulator of OM chitoporin ChiP(YbfM), Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base A (34/22);  new base C (2/2);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CTGATGGCGGTCCAAAAAAAAGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATCAGC  >  NC_000913/507127‑507412
                                                                                                                                                |                                                                                                                                             
cTGATGGCGGTCCAAAAAAAAGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATAtt                                                                                                                                           >  3:83579/1‑149 (MQ=255)
 tGATGGCGGTCCAAAAAAAAGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATAtt                                                                                                                                           <  1:433126/148‑1 (MQ=255)
     ggCGGTCCAAAAAAAAGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATATCGCTATTGGcc                                                                                                                                       <  2:841168/148‑1 (MQ=255)
           ccAAAAAAAAGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCAttttt                                                                                                                                 >  2:598471/1‑148 (MQ=255)
              aaaaaaaGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAg                                                                                                                             >  3:455826/1‑149 (MQ=255)
              aaaaaaaGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAg                                                                                                                             >  3:537235/1‑149 (MQ=255)
                         cATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACAttt                                                                                                                   >  2:459851/1‑148 (MQ=255)
                         cATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACAttt                                                                                                                   <  2:233946/148‑1 (MQ=255)
                          aTCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACACTCGAGTTACACCGTCGCTTAAAGTGACGGCATAGTAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTTGCCATTTTTTTAGCGCAACATTAGc                                                                                                                 <  4:561679/149‑1 (MQ=255)
                          aTCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTg                                                                                                                  <  2:539241/148‑1 (MQ=255)
                             ttGCCTAAGAGTATTGGAAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGCCGGCATAATAATAAGAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGc                                                                                                              >  3:544868/1‑149 (MQ=255)
                                ccTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCaaa                                                                                                           <  3:536996/149‑1 (MQ=255)
                                   aaGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAAtt                                                                                                         <  1:308204/148‑1 (MQ=255)
                                    agagTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTcc                                                                                                       <  3:521233/149‑1 (MQ=255)
                                     gagTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTccc                                                                                                      >  3:882259/1‑149 (MQ=255)
                                     gagTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTccc                                                                                                      <  4:482912/149‑1 (MQ=255)
                                     gagTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTcc                                                                                                       <  1:884696/148‑1 (MQ=255)
                                     gagTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTcc                                                                                                       <  2:820307/148‑1 (MQ=255)
                                         aTTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTtc                                                                                                   >  2:287069/1‑148 (MQ=255)
                                            ggCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCa                                                                                               >  4:318401/1‑149 (MQ=255)
                                            ggCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCa                                                                                               >  4:383818/1‑149 (MQ=255)
                                                 gATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAg                                                                                          >  4:520214/1‑149 (MQ=255)
                                                           ttGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGc                                                                                 >  1:372575/1‑148 (MQ=255)
                                                             gAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCacc                                                                              <  3:520214/149‑1 (MQ=255)
                                                               gCGACAATCTAGTTACACCGTCGCTTAATGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCaccga                                                                            <  3:764310/149‑1 (MQ=255)
                                                                cGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGACCAATAGCGATAATGGCCATTTTTTTAGCGCAACATTTGCGTCAAATTCCCTTCTCCATACAGGTGTAGTGCaccga                                                                            >  2:98930/1‑148 (MQ=255)
                                                                    aaTCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACcgcg                                                                       >  4:76783/1‑149 (MQ=255)
                                                                          gTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACcac                                                                  >  2:81489/1‑148 (MQ=255)
                                                                              cacCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACatat                                                              >  2:500835/1‑148 (MQ=255)
                                                                                   tCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGccc                                                         >  2:276477/1‑148 (MQ=255)
                                                                                       ttAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCcggcg                                                    >  3:892836/1‑149 (MQ=255)
                                                                                           aGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCa                                                 <  2:372575/148‑1 (MQ=255)
                                                                                                     taataataaAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTaa                                       <  1:598471/148‑1 (MQ=255)
                                                                                                        taataaAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTa                                   >  3:91814/1‑149 (MQ=255)
                                                                                                         aataaAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTa                                   >  1:161130/1‑148 (MQ=255)
                                                                                                          ataaAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAg                                                                                          <  4:201720/92‑1 (MQ=255)
                                                                                                          ataaAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAg                                                                                          >  3:201720/1‑92 (MQ=255)
                                                                                                            aaaaaaaTGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATcc                                >  1:538702/1‑148 (MQ=255)
                                                                                                            aaaaaaaTGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATcc                                >  1:856545/1‑148 (MQ=255)
                                                                                                              aaaaaTGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACa                                                                                           <  2:392902/87‑1 (MQ=255)
                                                                                                              aaaaaTGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACa                                                                                           >  1:392902/1‑87 (MQ=255)
                                                                                                                aaaTGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGcc                            >  1:509348/1‑148 (MQ=255)
                                                                                                                     aaaTTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCga                      <  4:455826/149‑1 (MQ=255)
                                                                                                                       aTTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCgag                     >  2:538097/1‑148 (MQ=255)
                                                                                                                         tCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGAt                   >  2:394728/1‑148 (MQ=255)
                                                                                                                            tctTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCg               >  4:413414/1‑149 (MQ=255)
                                                                                                                            tctTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCg               >  4:172603/1‑149 (MQ=255)
                                                                                                                             ctTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCg               >  1:338839/1‑148 (MQ=255)
                                                                                                                             ctTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCg               <  1:276477/148‑1 (MQ=255)
                                                                                                                             ctTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCg               >  2:492092/1‑148 (MQ=255)
                                                                                                                                tGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGa           <  4:83579/149‑1 (MQ=255)
                                                                                                                                  aCGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTATTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATg         >  4:552073/1‑149 (MQ=255)
                                                                                                                                     ggCCAATATCGCTATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCAt      >  3:57526/1‑149 (MQ=255)
                                                                                                                                       ccAATATCGCTATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTg            >  2:548454/1‑141 (MQ=255)
                                                                                                                                       ccAATATCGCTATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTg            <  1:548454/141‑1 (MQ=255)
                                                                                                                                       ccAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATc     >  2:455048/1‑148 (MQ=255)
                                                                                                                                       ccAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATCa    <  3:76783/149‑1 (MQ=255)
                                                                                                                                        cAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATCa    <  1:459851/148‑1 (MQ=255)
                                                                                                                                          aTAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATCAGc  <  1:287069/148‑1 (MQ=255)
                                                                                                                                          aTAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATCAGc  >  1:174856/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CTGATGGCGGTCCAAAAAAAAGAGTCATCTTGCCTAAGAGTATTGGCAGGATGGTGAGATTGAGCGACAATCGAGTTACACCGTCGCTTAAAGTGACGGCATAATAATAAAAAAATGAAATTCCTCTTTGACGGGCCAATAGCGATATTGGCCATTTTTTTAGCGCAACATTTGCGGCAAATTCCCTTCTCCATACAGGTGTAGTGCACCGACCGCGACCACATATCGCCCCGGCGGCATGGCGCGTAATTTATCCCGCCAGGCGAGATTTCGCTGATGCATCAGC  >  NC_000913/507127‑507412

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