Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,435,075 T→A 15.8% G144G (GGT→GGA dapD → 2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,435,0750TA15.8% 131.5 / 2.8 61G144G (GGT→GGAdapD2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/27);  new base A (9/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAAAGGT  >  CP000730/1434936‑1435213
                                                                                                                                           |                                                                                                                                          
gTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGATTTCGAATGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                                                                           >  1:851246/1‑141 (MQ=255)
   cGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATAGGCGCAACAATTAATATTGGCGCAGCCGTTGGCGAAGGTACCATGATTGATATGAATGCGACTCTCGGTGGTCGGGCTACAACTGGAaaaa                                                                                                                                        >  1:200950/1‑141 (MQ=255)
             agaACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTATAATTGATAAAAATGTAGATGt                                                                                                                              >  2:246171/1‑141 (MQ=255)
             agaACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACCACTGGAAAAAAAGTACATga                                                                                                                              >  2:670430/1‑140 (MQ=255)
                  aacacatatttGAATATGGTCCTGTCGTTATTAATCGGGAAACAATTATTTTTGGCGCAGTATTTGGCGAAGGTATTATGATTGATAGGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGGACATGTAtggg                                                                                                                         <  2:101788/132‑1 (MQ=255)
                  aaGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgggg                                                                                                                         <  1:181776/141‑1 (MQ=255)
                  aaGCCATAATTGAAGATGGTGTTGTCGTTATGAGGGGCGCAACAGTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgggg                                                                                                                         <  2:757670/141‑1 (MQ=255)
                    accATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCt                                                                                                                       <  1:815367/140‑1 (MQ=255)
                              taGATGGTGTGGTCGGTATGATGTGGGCAACACTTAATAGTGGAGCAGTCGTTCGCGAAGGTACAATGATTGATATGATTGCTACTCTCGGTGGTCGAGAAAAAAAAAAAAAAAATGTACATGTAGGGGCTGGCGCAGTAt                                                                                                             <  2:65885/140‑1 (MQ=255)
                               aGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                                                                            >  1:333602/1‑141 (MQ=255)
                                  tGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                                                                                         <  1:715179/141‑1 (MQ=255)
                                    gcccTGTCGCTAAGATGGTGGTAACAAATAGTATTGGCGCCGTCGTTGGCGAAGGGACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTAAAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAg                                                                                                       <  2:462030/138‑1 (MQ=255)
                                           tggtgtgaggGGCGCAACGATTTGTATTTGGGCAGTTGTTGGGGAAGGTACAATGATTGATATGAATGTTACTCTCGGTGGTCGTGCTACTACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                                                                <  2:426776/136‑1 (MQ=255)
                                           cgttgggatgGCCGCAACAATTAAAATTGGCGCAGTTGGTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                                                                <  2:498433/135‑1 (MQ=255)
                                           cGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                                                                                >  1:694196/1‑141 (MQ=255)
                                            gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTCAAAATGTACATGTAGGGGCTGGCGCAGTCTTAGCAGGTGTGAtt                                                                                               >  1:58551/1‑141 (MQ=255)
                                              tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTAGAAAAAAAGAACAAGTAGGGGCGGGCGACGGATTACCAGGAGTGATGGa                                                                                             >  2:575101/1‑141 (MQ=255)
                                              tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTATCAGGAGCGATGGa                                                                                             >  2:211395/1‑141 (MQ=255)
                                                 gatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGCATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGAGCAGTATTAGCAGGTGTGATTGAAcc                                                                                          >  1:530087/1‑141 (MQ=255)
                                                      gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTa                                                                                     >  1:621834/1‑141 (MQ=255)
                                                       cacaacaaTTAATGTTCGCGCAGCCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                                                                    <  2:622065/139‑1 (MQ=255)
                                                        gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTATCATTAAAAAAAAAGTACATGTAGGGGCTGGAGCAGTATTAGAAGATGAGATTGAAGCCCATAGt                                                                                   >  2:345375/1‑141 (MQ=255)
                                                        ccaacaaTTAGTATTGGCGCAGTCGTTCGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                                                                                   <  2:158608/140‑1 (MQ=255)
                                                         caacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAAGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                                                                  >  1:34192/1‑141 (MQ=255)
                                                           acaaTTCATACTGGCGCAGTCGTTGGCGAGGGTACAATGATTGATATGTATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTTGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCt                                                                                <  1:821434/141‑1 (MQ=255)
                                                            aaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCtt                                                                               <  1:554711/140‑1 (MQ=255)
                                                               ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAATAAAAAAGCGCATGAGGGGACTGCCGAAGCATCAGCAGGTGTAGTTGAAccacatcatcgctcac                                                                            >  2:93165/1‑127 (MQ=255)
                                                               ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAGTACAAGTAGGGGCTGGCGCAGTATTAGCAGGTGGGATTGAACCCGCTAGTTCTTCAc                                                                            >  2:455297/1‑141 (MQ=255)
                                                                   tatTGGCGCAGTCGTGGGCGAGGGTACAATGATGGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                                                                        <  2:23685/141‑1 (MQ=255)
                                                                    atTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTTGTGCTTCACCGGtt                                                                       >  1:90616/1‑141 (MQ=255)
                                                                       ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCTAGTGCTTCACCGGTTATaa                                                                   >  1:802903/1‑141 (MQ=255)
                                                                         cgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAt                                                                  >  1:419786/1‑141 (MQ=255)
                                                                          ggagtagtagtccaaGGTAGAATGATTTGTATGAATGCTACGGGCGTTGGTCGTGCTACAAGTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTGGCAGGTGTGATTGAGCCCCCTAGGGCTTCACCGGTTATAATc                                                                 <  2:29419/128‑1 (MQ=255)
                                                                             gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCCCGGGGTATAGACGAg                                                              >  1:323148/1‑141 (MQ=255)
                                                                                  tGGCGAAGGGACAATGATTGATATGAATGTTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatga                                                         <  2:40343/141‑1 (MQ=255)
                                                                                   ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAGAAAAAAAAAGCATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCGACTAGTGCTTCACCGGTTATAATCGAGgatgag                                                        >  2:657250/1‑140 (MQ=255)
                                                                                    accaaGGGACAAAGAATGATATAAAAGCTACTCTCGGTCGTCGTGCTACAACTGGAAAAAATGTACATGTAGGGGCCGGCGCAGTATAAGCAGGTGTGATTGAACCCCCTAGTGCGTCACCGGTTATAATCCAGgaaaatg                                                       <  2:148443/138‑6 (MQ=255)
                                                                                             cAATGATTGATATGAATGCTACTCTCGGTGGTAGTGTCCCAACAGGTACACATTTACGTGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTGTAATCGAGGATGATGTATTAATCg                                              >  1:856280/1‑141 (MQ=255)
                                                                                               aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGt                                            <  1:449844/141‑1 (MQ=255)
                                                                                               aTGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAGAAAAAATATACATGGAGGGGATGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTCGTGCTTCACCGGTTATAATCGTGGACGATGTATTacgcggc                                            >  2:773004/1‑135 (MQ=255)
                                                                                                      atatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATg                                     <  2:621834/141‑1 (MQ=255)
                                                                                                      atatGAATGCGACTCTCGGTGGTCGGGCTACAACTGGTAAAAATGTCCATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATg                                     <  2:318943/141‑1 (MQ=255)
                                                                                                       tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGCGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCaaaaga                                    >  1:280860/1‑138 (MQ=255)
                                                                                                       aatGAATGCTACTCGCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGATTCACCGGTCATAATCGAGGATGATGTATTAATCGGTGCAAATGc                                    <  1:250900/140‑1 (MQ=255)
                                                                                                          gAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGt                                 <  1:670430/141‑1 (MQ=255)
                                                                                                                tACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAACAGGAGTGATTGAACACCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTCAtttt                           >  2:210907/1‑141 (MQ=255)
                                                                                                                   tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGa                        >  1:488616/1‑141 (MQ=255)
                                                                                                                       ggtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTCTTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATCATTGAGGAGGTTTTATTAATCGGTGCAATTGCAGTTTTTTTAGAAggg                    >  1:586788/1‑140 (MQ=255)
                                                                                                                        gtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGgtg                   <  2:851246/141‑1 (MQ=255)
                                                                                                                        gcggtggtGCTACAACTGGTAAAAATGTACATGTGGGGGTGCGGTCACTGTTAGGGTGGGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGgtg                   <  2:303404/139‑1 (MQ=255)
                                                                                                                             cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACgt              >  1:190546/1‑141 (MQ=255)
                                                                                                                              gTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACgtg             >  1:97328/1‑141 (MQ=255)
                                                                                                                               tGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACgtgt            >  1:101366/1‑141 (MQ=255)
                                                                                                                                 ttACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTg          <  1:861580/140‑1 (MQ=255)
                                                                                                                                 cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTg          <  1:218582/141‑1 (MQ=255)
                                                                                                                                 cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGAGCAAATGCAGTTATTTTAGAAGGTGTAGGTGTTg          >  2:481265/1‑141 (MQ=255)
                                                                                                                                  caCAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTgg         <  1:72040/140‑1 (MQ=255)
                                                                                                                                   aCAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGt        <  1:674163/141‑1 (MQ=255)
                                                                                                                                    cAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTa       <  1:254213/141‑1 (MQ=255)
                                                                                                                                     aaCTGTGAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAACGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTaa      <  2:630272/141‑1 (MQ=255)
                                                                                                                                         ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGGGCTTAACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAAAGGt  >  2:257179/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
GTGCGTTTATTAGAGAACAAGCCATTATTGAAGATGGTGCTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTATTAATCGGTGCAAATGCAGTTATTTTAGAAGGTGTACGTGTTGGTAAAGGT  >  CP000730/1434936‑1435213

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