Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,988,051 A→T 11.7% intergenic (‑92/‑171) azuC ← / → yecR acid‑inducible small membrane‑associated protein/lipoprotein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,988,0510AT11.7% 140.2 / 6.6 58intergenic (‑92/‑171)azuC/yecRacid‑inducible small membrane‑associated protein/lipoprotein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base T (2/5);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTT  >  NC_000913/1987903‑1988198
                                                                                                                                                    |                                                                                                                                                   
tttGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTACGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGGGATACaa                                                                                                                                                     <  1:767006/149‑1 (MQ=255)
 tagAACGTCTTGCAATAGTTATTGAACATACTATTCAGGATTTAGCGCAGTTTCATAGCGTCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATAAATCTTTATGTGATACaaa                                                                                                                                                    <  2:546864/147‑1 (MQ=255)
  tGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAATATCAATCTTTATGTGATACAAAt                                                                                                                                                   >  2:248565/1‑149 (MQ=255)
  tGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAAt                                                                                                                                                   >  1:522650/1‑149 (MQ=255)
        tCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATa                                                                                                                                             >  1:467082/1‑149 (MQ=255)
         cTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATaa                                                                                                                                            >  1:1110347/1‑149 (MQ=255)
         cTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATaa                                                                                                                                            >  1:483836/1‑149 (MQ=255)
          ttGGAATAGTTATTGAACATACTATTCAGGATTTTGCGCAGTTTCATCGCGGCACTGCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAACTCAATCTTTATGTGCTTTGTATCACATaaa                                                                                                                                           <  1:750455/149‑1 (MQ=255)
          ttGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTTATCGCGGCACTCCGACGATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATTATATGACAACCATCACAAAAATCAATCTTTATGTGTTACAAATCACATaat                                                                                                                                           <  2:277852/149‑2 (MQ=255)
          ttGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTACGACCATTTGTTATACAGGTTTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaa                                                                                                                                           <  2:197948/149‑1 (MQ=255)
             cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaagat                                                                                                                                        <  1:337927/149‑4 (MQ=255)
             cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaagat                                                                                                                                        <  1:652340/149‑4 (MQ=255)
             cAATACTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTAAAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaagat                                                                                                                                        <  1:303742/149‑4 (MQ=255)
              aaTAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATAcc                                                                                                                                       >  2:888514/1‑149 (MQ=255)
                       tGAACATACTTTTCCGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAATTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATg                                                                                                                              <  1:972480/149‑1 (MQ=255)
                         aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                            <  2:41429/149‑1 (MQ=255)
                         aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                            <  2:522650/149‑1 (MQ=255)
                         aaCATACTTTTCAGGATTTTGAGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGtt                                                                                                                            <  2:9519/149‑1 (MQ=255)
                           cATACTTGTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTtat                                                                                                                          >  1:81197/1‑149 (MQ=255)
                                  ttCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACATATCACATAAATACCCCTTTAATGTTATATAAAtt                                                                                                                   >  1:316896/1‑148 (MQ=255)
                                    cAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACAGAAATACCCCTTTAATGTTATAAAAATGAt                                                                                                                 >  2:696568/1‑149 (MQ=255)
                                     aGGATTTTGCGCAGTTTCATCTCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATa                                                                                                                <  1:1183265/149‑1 (MQ=255)
                                              cgcAGTTTCATCGCGGCGCTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATAACATAAATACCCCATTAAGGTTATAAAAATGATAATCaaaaaa                                                                                                       >  1:491575/1‑149 (MQ=255)
                                                   tttCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGcc                                                                                                  <  1:222276/149‑1 (MQ=255)
                                                      cATCGCGTCACTCCGAGCAGTTGTTATACAGGTGTTATTGTCTTTGCGCTTAAAATATGACAACCATCACAAAACTCAATCATTATGTGATACAAATCAAATAAATACCCGTTTAATGTTATAAAAATGATAATCAAAAAACAGacccc                                                                                               >  2:1048166/1‑149 (MQ=255)
                                                        tCGCGGCAATCCGAACATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGAGAACCATCAAAAAAATCAATCTTTATGTGTTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCt                                                                                             <  2:562285/149‑1 (MQ=255)
                                                            ggCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCAAttt                                                                                         <  2:390280/149‑1 (MQ=255)
                                                                cTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                     <  1:508238/149‑1 (MQ=255)
                                                                 tCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGac                                                                                    <  2:172667/149‑1 (MQ=255)
                                                                      ccATTTGTTATACAGGTGTTATTGTCTTTGCTCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTa                                                                               <  2:620512/149‑1 (MQ=255)
                                                                         tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAg                                                                            >  1:580608/1‑149 (MQ=255)
                                                                         tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAg                                                                            <  1:888514/149‑1 (MQ=255)
                                                                         tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAg                                                                            >  2:937152/1‑149 (MQ=255)
                                                                            gTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATg                                                                         >  2:151552/1‑149 (MQ=255)
                                                                             ttATACAGGTGATATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTAAGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATgg                                                                        <  2:489888/149‑1 (MQ=255)
                                                                                   aGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTTATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaa                                                                  <  1:922037/149‑1 (MQ=255)
                                                                                   aGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaa                                                                  >  1:943221/1‑149 (MQ=255)
                                                                                     gtgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAAt                                                                >  2:340440/1‑149 (MQ=255)
                                                                                     gtgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAAt                                                                >  2:36298/1‑149 (MQ=255)
                                                                                      tttttatGTCTTTGAGATTCTAATATCACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATa                                                               <  1:1092601/143‑1 (MQ=255)
                                                                                       tttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAg                                                              <  1:34737/148‑1 (MQ=255)
                                                                                        ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATTATAATCAAAAAACAGCCACCCTATTTCTGACACCTACAGATGAAAAGAAATAgc                                                             >  1:819541/1‑149 (MQ=255)
                                                                                        ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCACCAATTTCAGACACCTACAGATGGCAAGAAATTgc                                                             >  1:894924/1‑149 (MQ=255)
                                                                                         tATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAgcg                                                            >  2:50382/1‑149 (MQ=255)
                                                                                             gTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCATAAAACAGCCCCCCTATTACTGACACCTACAGAGGGCAAGAAATAGCGCCTg                                                        >  2:444000/1‑149 (MQ=255)
                                                                                                      cTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTc                                               <  2:580608/149‑1 (MQ=255)
                                                                                                         ataataTGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATTTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGAGCCTGCCAGGAGTCttt                                            >  1:909786/1‑149 (MQ=255)
                                                                                                           aataTGACAACCACCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTc                                          <  1:296746/149‑1 (MQ=255)
                                                                                                            atatGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAAGGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTCACACCTACAGATGGCAATAAATAGCGCCTGCCAGGAGACTTTTcc                                         >  2:1074995/1‑149 (MQ=255)
                                                                                                               tGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCTCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGc                                      <  2:483836/149‑1 (MQ=255)
                                                                                                                     ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCAAGGCGTCTTTTCCGGCCATTGt                                >  2:1152896/1‑149 (MQ=255)
                                                                                                                        tCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAATACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgc                             <  2:81197/149‑1 (MQ=255)
                                                                                                                          acaAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGATATAAAAATGAAAATCAAAAAACAGCCCCCCTATTACAGACCCCTACAGATGGAAAGAAATAGCGCCTGCCAGGCGTCTTTACCGGCCATTGTCgcag                           >  1:909748/1‑149 (MQ=255)
                                                                                                                                  caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTaa                   >  2:830902/1‑149 (MQ=255)
                                                                                                                                   ggtatttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAc                  >  2:248315/5‑149 (MQ=255)
                                                                                                                                       tttATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGt              <  2:1102700/149‑1 (MQ=255)
                                                                                                                                            gtgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAAt         >  1:537156/1‑149 (MQ=255)
                                                                                                                                             tgATACAAATCACATAAATACCCCTTTACTGTTATAAAAATGATAATCAAAAAACAGCCACACTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTGATCCGGCCATTGTCGCAGCACGGTAACGCGTAAAATa        >  2:922431/1‑149 (MQ=255)
                                                                                                                                             tgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATa        >  2:408761/1‑149 (MQ=255)
                                                                                                                                               aTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGt      <  2:467082/149‑1 (MQ=255)
                                                                                                                                                   aaaTCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCtt  >  1:633348/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGCTT  >  NC_000913/1987903‑1988198

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