Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,075,234 C→T V1217I (GTA→ATA)  putA ← fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,075,2340CT94.3% 88.5 / ‑5.7 35V1217I (GTA→ATA) putAfused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/15);  total (19/16)
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.69e-01

ATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAA  >  NC_000913/1074987‑1075480
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
atatACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCTAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGCTGCAGCGCGTCATCCGGCCACAGTATCTg                                                                                                                                                                                                                                                     >  1:269428/1‑251 (MQ=255)
                aaaGGATATTGCTTTCGCCCCGGGAAAATCCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCACCTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGc                                                                                                                                                                                                                                     <  1:36294/250‑1 (MQ=255)
                  cggATATTGATTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGGGCCTTCCCGCGCGGCAACTGCTTCACCCAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCAc                                                                                                                                                                                                                                   >  2:40528/2‑251 (MQ=255)
                    gATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGGTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCGGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAcgg                                                                                                                                                                                                                                 >  2:56813/1‑251 (MQ=255)
                                         aaaaCCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGTAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCAt                                                                                                                                                                                                            >  1:151131/1‑251 (MQ=255)
                                                                    aTCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACacaaca                                                                                                                                                                                 <  1:119992/251‑1 (MQ=255)
                                                                           gcgGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTc                                                                                                                                                                          >  1:102244/1‑251 (MQ=255)
                                                                             ggCAACTGCTTCACACAAGGGGCGAGGCGGATCCGAATCACCGTGGAAGATCCCCGCACAAAACGGTTGGGCGGGAATATTTTCCGCTTTCGCCATTTGAATACTGTCGTTGACTGCCGATGGCAATGCCTTCACAAACTGACGATGCACCGCGCCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGTCGGCGAGCTGATTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCTCTCAc                                                                                                                                                                        <  2:261268/251‑1 (MQ=255)
                                                                              gcaccTGCTTCACACAATGCGCGAAGCGGATCCGAATCACCGGGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACGAACGGACGATGCAGCGCGTCACCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCAcg                                                                                                                                                                       <  2:102244/247‑1 (MQ=255)
                                                                                    gCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGgc                                                                                                                                                                   >  1:77691/1‑249 (MQ=255)
                                                                                           acaAGGCGGGAAGCTGAGCCGAATCCCCGTGGAAGATCCCCGCATCAACCGTTGGAGCGGTTATATTTTCCGCTTTCGCAAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTc                                                                                                                                                          <  2:197910/251‑1 (MQ=255)
                                                                                             aaTGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTcc                                                                                                                                                         >  1:182122/1‑250 (MQ=255)
                                                                                                          gATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTc                                                                                                                                            >  1:47325/1‑250 (MQ=255)
                                                                                                                     ccGTGGAAGACCACCCCAACAAACGGTTGAGGGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACGGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGccc                                                                                                                                <  2:1686/251‑1 (MQ=255)
                                                                                                                       gTGGAAGTTCCCCGCATCAACCGGTTGGGCGTTTATATTTGCCGCTTTCGCCGGTTGAATACGTTCGTTGACTGCCGATGGCAAGGCCTTCACTACCTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCg                                                                                                                              <  2:6551/251‑1 (MQ=255)
                                                                                                                       gTGGAAGATCACCGCATCACTCGGTTGAGCGGTTATATTTTCCGATTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATATGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGAGCTGCGTTAACACGTCGGCCCCg                                                                                                                              >  2:265781/1‑251 (MQ=255)
                                                                                                                       gTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCg                                                                                                                              >  1:172714/1‑251 (MQ=255)
                                                                                                                                   cGCATCAAACTGTGGAGCGGTGATATTTTCCCCTTTCGCCAGTTGAACACGTTCGCTGACTGCCGTGGTCAAGGCCTGCACCAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCGCGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAAGACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGGGTTGCGTTCACCCGTCGGCCCCGGCAGCATTCGtt                                                                                                                  <  2:116149/251‑1 (MQ=255)
                                                                                                                                             ggTTGAGCGGTTATAGTTTCCGCTGTCGCTGGTTAAATCCGTTCGCTGACTGCTGATGGCAATTCCTTCACTAACGGCCGATGCAGCGTGTCTTCCGGCCACATTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCTGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGGGTTcctg                                                                                                         <  2:145510/250‑1 (MQ=255)
                                                                                                                                                   gCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGATGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTgcgc                                                                                                  >  1:252877/1‑251 (MQ=255)
                                                                                                                                                              ttCCGCTTTCGCCAGTTGAATCCTTGCTCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCCCAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCa                                                                                        <  1:56813/250‑1 (MQ=255)
                                                                                                                                                                     gccgctatcggaatccgctccgtGACTGCCTATGGCAATGCCTTCCTTAACCGACGAGGCTGCGCGCCTCCCGGCCACAGTTTCTGGCTCCCCACCGCCCGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCAGCGTCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTGGCGTGCACCCGTCGGCCCCGTCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGc                                                                                  <  2:269428/227‑1 (MQ=255)
                                                                                                                                                                           aGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCAc                                                                                                                                                                        >  2:90300/1‑157 (MQ=255)
                                                                                                                                                                           aGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCAc                                                                                                                                                                        <  1:90300/157‑1 (MQ=255)
                                                                                                                                                                                 aTACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTc                                                                                                                                     >  1:172298/1‑186 (MQ=255)
                                                                                                                                                                                            gACTGCCGATGGTAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCGGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGAc                                                         <  2:62661/251‑1 (MQ=255)
                                                                                                                                                                                              cTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACg                                                        >  1:58499/1‑250 (MQ=255)
                                                                                                                                                                                                             gCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTc                                         >  1:116317/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaCTGACGATGCAGCGCGCCATCCGGCCACCGTATCTGGCTGCCCACCGCCAGAACGGCGGCGAGCTGAGTCAGCGAATCCTGCTCATACTCGGCAATACACAACACGCGCTAACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCAGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCCTATTGCGTACCTAACGCCTGCAATTATGGACGATTTACTGCCCATTCCCACAGTGc                                >  2:192392/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGc                                >  1:189600/1‑250 (MQ=255)
                                                                                                                                                                                                                          taCGATGCAGCGCTTCACCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTTCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCAtt                             <  1:226255/249‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCGCTTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGGCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCTTCCAGGTGTTTCGTTCACCCGTCGGCCCCGGCATCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCt                     <  1:268449/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cATCCAGCCCCAGTATCTGGCTGCCAACCGCCAGCACGCCGGCGAGCTGAATCAGCGCATCCTGCTCATCATCGGCAATACACCACACGCGCTCACGCGGCCGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCCTTTAGCGGCTGAGTAAAt             >  2:239195/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    aGTATCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCa   >  2:183797/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    aGTACCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATTCACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCaa  <  1:183946/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
ATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCACGGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAA  >  NC_000913/1074987‑1075480

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