Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,305 C→T R395C (CGT→TGT)  icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,3050CT98.3% 161.1 / NA 60R395C (CGT→TGT) icdisocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/25);  minor base A (1/0);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGT  >  NC_000913/1196167‑1196443
                                                                                                                                          |                                                                                                                                          
gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGACGCGGCTGACTTAATTGTTCAAGGTATGCAAGGCGCACTCCACGCGAAAACCGTCACCTATGACTTCGCGTGTCt                                                                                                                                        >  1:66991/1‑143 (MQ=255)
   aaaGTAAGTCCTGGCTCTGTTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgat                                                                                                                                     <  1:253533/143‑1 (MQ=255)
      gTAAATCCTGGCTCTATTATTCTCTCGGCTGAGATGATGCTGAGCCACATGGGCTGGACCGAAGCGGCTGACTTAATTGTTCAAGGCATGGAAGGCCCAATCACCGCCAAAACAGTAACATATGACTGCGAGAGGCTgatggg                                                                                                                                  >  1:483534/1‑142 (MQ=255)
      gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAGCCGTAACCTATGACTTCGGGTGTCTgatgga                                                                                                                                  >  2:660000/1‑143 (MQ=255)
         aaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGatagatgg                                                                                                                               >  1:218038/1‑143 (MQ=255)
           tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                                                                                                             >  2:871252/1‑143 (MQ=255)
                 ctctATTCTTCTGTCCGCTGATATGATGCTGCGCCACAGGGGTTGGACCGAAGCGGCTGACGTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAAc                                                                                                                       <  1:186389/143‑1 (MQ=21)
                  tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGTACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAGGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAAct                                                                                                                      >  1:441873/1‑143 (MQ=255)
                      ttgtgcgcgccGCTGAGATGATGCGGCGCCGCAGGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAActgctg                                                                                                                  <  2:432777/134‑1 (MQ=21)
                       tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGACGCGGCTGCCTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGATAAACTGCTGa                                                                                                                 >  1:425737/1‑143 (MQ=255)
                       tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCCGa                                                                                                                 >  1:303355/1‑143 (MQ=255)
                             ccccgctcgggtgacgCTGCGCCACATGGGGCGGACCGAACCGGCTGACTTAATTCTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                                                                                                           <  2:651838/132‑1 (MQ=18)
                              tcCGCTGCGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTc                                                                                                          <  2:212073/143‑1 (MQ=255)
                              tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTc                                                                                                          <  2:718730/143‑1 (MQ=255)
                               ccGCTGAGATGAGGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                                                                                                         <  2:880907/143‑1 (MQ=255)
                                cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCag                                                                                                        >  2:128257/1‑143 (MQ=255)
                                   tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGCGTGTCTGATGGATGGCGCTAAACTGCTCACATGTTCAGAGt                                                                                                     >  1:774352/1‑143 (MQ=255)
                                     agaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt                                                                                                   <  2:650317/143‑1 (MQ=255)
                                     agaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt                                                                                                   <  2:149623/143‑1 (MQ=255)
                                      gatgggttttgcccACGTGCGGTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGTGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                                                                                  <  1:116863/131‑1 (MQ=17)
                                        ggctgctgcGCCACAGTGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTTGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGt                                                                                                <  1:973769/140‑1 (MQ=18)
                                         gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTg                                                                                               >  1:994712/1‑143 (MQ=255)
                                         gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTg                                                                                               >  2:1054207/1‑143 (MQ=255)
                                                    acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                                                                    <  2:1012585/143‑1 (MQ=255)
                                                        gggTTGGGCCGTAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                <  1:507010/143‑1 (MQ=255)
                                                        gggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGGTTGGTGACGCGATCATCGaaa                                                                                >  2:459166/1‑143 (MQ=255)
                                                           ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCggacctc                                                                             >  2:319051/1‑137 (MQ=255)
                                                            tGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACAt                                                                            <  2:278812/143‑1 (MQ=255)
                                                                ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                                                                        >  2:727216/1‑143 (MQ=255)
                                                                 cGAAGCGGCTGGCTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAt                                                                       <  1:517987/143‑1 (MQ=255)
                                                                  gAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGACAACCTGTAATg                                                                      >  2:225363/1‑143 (MQ=255)
                                                                  gAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGGTTGGTGGCGCGATCATCGAAAACATGTAATg                                                                      >  2:135366/1‑143 (MQ=255)
                                                                   aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGAAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGc                                                                     >  2:722705/1‑143 (MQ=255)
                                                                        gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAg                                                                >  2:521342/1‑143 (MQ=255)
                                                                            aCTTAATTCCTAAAAGTATGGGAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGAGGGCGCTAAACGGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTg                                                            <  2:692376/143‑1 (MQ=18)
                                                                              ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTCACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAATTTGGTGACGCGATCATCGAAACAAGGTTATGCCGTGGTCTGtt                                                          >  1:646893/1‑143 (MQ=18)
                                                                              ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAACACATGTAATGCCGGAGTTTGtt                                                          >  1:576269/1‑143 (MQ=255)
                                                                              ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                                          >  2:1043770/1‑143 (MQ=255)
                                                                              ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGAATTCGAGTGTCTGATGGATGGCGCTAAACTGCTGACATGTTCAGAGATTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTgac                                                          >  1:166272/1‑141 (MQ=255)
                                                                                 aTTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa                                                       <  2:250441/143‑1 (MQ=255)
                                                                                 aTTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa                                                       <  2:852998/143‑1 (MQ=255)
                                                                                    gttatttgTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTGGTTAAAttg                                                    >  1:1045845/8‑142 (MQ=255)
                                                                                    gTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAttt                                                    >  2:542113/1‑143 (MQ=255)
                                                                                     ttAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATtta                                                   >  1:285858/1‑143 (MQ=255)
                                                                                      tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACACGTCATGTCGTAGTTTCTTAAAtgcag                                                  >  2:277851/1‑139 (MQ=18)
                                                                                      tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGAATTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTAAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATttat                                                  >  1:611304/1‑143 (MQ=255)
                                                                                       aaaGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTCATGCCGTAGTGTGTTAAATttatt                                                 >  1:721988/1‑143 (MQ=255)
                                                                                         aGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAACGCCCGAGTTAATAAATTtaataaa                                              >  2:752052/1‑142 (MQ=17)
                                                                                         aGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAAATATGATAACGGTTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTaa                                               >  1:238936/1‑143 (MQ=21)
                                                                                            tATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                            <  1:221114/143‑1 (MQ=255)
                                                                                                      gcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                  <  2:994712/143‑1 (MQ=255)
                                                                                                           caccgcgAAACCCGTAACCTATGCCTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTccc                             <  2:348443/140‑1 (MQ=255)
                                                                                                             acgcgcagtccGTATCCGAGGACTTCGAGTGTCTGATGGATGGCGCTTAGCTGCTGGAATGTTCAGCGTGTGGTGACGCGATCATCGAAAACATGTAGTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                           <  1:20805/134‑1 (MQ=255)
                                                                                                                     cccGTACCCTATGACTTCGAGTGTCTGACGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaataaat                   <  2:127322/142‑8 (MQ=25)
                                                                                                                     aCCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGtttttt                   <  1:727216/143‑1 (MQ=255)
                                                                                                                       cGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGt                 >  2:34834/1‑143 (MQ=255)
                                                                                                                           aCCTATGACTTCGAGTGTCTGATGGCTGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg             <  2:630428/143‑1 (MQ=255)
                                                                                                                              tATGACTTCGAGTGTCTGATGGATGGCGATGAACTGCTGAAATGTTCAGTGTTTGGAGCCGCCACCAAAGACCAGTGATAATGCCGTAGTTTGTTTAATTTATTAACGGGAGCGTAACCCTCCCGTTGTTGTTTGTTAGgccg          >  1:1090840/1‑141 (MQ=255)
                                                                                                                                  aCTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTCTTTTGTTAGGCTGCTaa      >  1:1109428/1‑143 (MQ=255)
                                                                                                                                      cGAGTGTCTGATGGATGGCGCTAAACTGCTGCAATGTTAAGAGTTTGGTGAAGCGATCATCGAAAACCTGTAATGCCGTAGTTTGTTAAATTTATTGACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGGTGCTAACggc  >  2:754735/1‑142 (MQ=37)
                                                                                                                                          |                                                                                                                                          
GACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGT  >  NC_000913/1196167‑1196443

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