Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,203 T→C A305A (GCA→GCG lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,2030TC93.5% 85.8 / ‑1.7 31A305A (GCA→GCGlplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (15/14);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTC  >  NC_000913/4622964‑4623449
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
cGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGTAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTAc                                                                                                                                                                                                                                               >  1:278554/1‑249 (MQ=255)
                                                                                                                            tttaatcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAgct                                                                                                                  >  1:48782/15‑248 (MQ=255)
                                                                                                                                  cagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGt                                                                                                           >  1:283879/9‑251 (MQ=255)
                                                                                                                                   agcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGtt                                                                                                          >  1:575533/8‑251 (MQ=255)
                                                                                                                                     cgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTCTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTc                                                                                                         >  1:565100/6‑250 (MQ=255)
                                                                                                                                       agttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATc                                                                                                      >  1:66009/4‑251 (MQ=255)
                                                                                                                                         ttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCCGCTAAAGCGTTCATccc                                                                                                    >  2:286819/2‑250 (MQ=255)
                                                                                                                                             cTTACAGCCCACGCTGTCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCCACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGAGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGGGGTATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCCCACGCGTTCATCcaaca                                                                                                 >  2:546301/1‑247 (MQ=255)
                                                                                                                                               taccgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCcagcag                                                                                                <  1:405603/245‑1 (MQ=255)
                                                                                                                                                accgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGAt                                                                                              <  1:199255/247‑1 (MQ=255)
                                                                                                                                                accgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGAGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATg                                                                                             >  1:656407/4‑251 (MQ=255)
                                                                                                                                                 tcgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGATCCGAGAAAGCAACCAACAGCGCTTCGCACTCTTGTTGCAGAATATCCGCGCGGTCCAGGCAGCCTTGCAGTCGCCCGGCGCGGTCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAACCACCTGGGCGCGGGCGATCTGGCCTTTTTCAACGTCGCACTACCGTTCCACGCCGCCCCAGCTAAAGCTATCATCcagccgttgg                                                                                            >  2:496789/3‑245 (MQ=255)
                                                                                                                                                                acgtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGcc                                                                             >  1:293981/9‑251 (MQ=255)
                                                                                                                                                                 cgtcgacaccttccGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCCGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCt                                                                            <  1:473921/239‑1 (MQ=255)
                                                                                                                                                                 cgtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACCCCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCt                                                                            <  1:206825/244‑1 (MQ=255)
                                                                                                                                                                            cccGTAGTTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCGTGCAGTCGCCCGGCGAGGGCTTCCAGCGGAGCGGGGTTGAGGCTGTCGGTAAACACCTGTGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTc                                                                 <  2:596114/251‑1 (MQ=255)
                                                                                                                                                                              cGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTcc                                                                >  1:545451/1‑250 (MQ=255)
                                                                                                                                                                                 aGCTCTTTTTCCTTGTCCGGGAAGTCATCCCACCGCGCCTCGCCCTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTc                                                            <  2:451381/251‑1 (MQ=255)
                                                                                                                                                                                     ctTTGTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTccc                                                          <  2:334180/249‑1 (MQ=255)
                                                                                                                                                                                        tttCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCt                                                      >  1:439276/1‑250 (MQ=255)
                                                                                                                                                                                                 ccGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACgg                                            <  1:345929/251‑1 (MQ=255)
                                                                                                                                                                                                         tCAACCAACAGCGCTTCGCACTCTGGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGTCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCAACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGCGCCTGACCGAATTTCGATTCCCAGCTACTCTGACGGGCAAAGGt                                    >  2:221577/1‑251 (MQ=255)
                                                                                                                                                                                                           aaccaacAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAgggtt                                  >  1:381111/1‑248 (MQ=255)
                                                                                                                                                                                                             ccaccAGCCCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGTGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCTGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCAGCCAGCAGATGCGAGAAGGCCGGAGCCTGACCGAAGTTCCATTCCCAGTTATTCTGACGGGCAAATGTTTCg                                <  2:269690/247‑1 (MQ=255)
                                                                                                                                                                                                                acAGCGCTGCGCACTCGGGTTGCAGCATACCCGCGCGGGACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCATCCAGCGGCGCGGGGTTGAGGCTGTCGGTAACCACCTGGGCGCTGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAACTAAAGCGTTCATCCAGCATATGCGAGAATGCCGTGGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCg                             <  2:52185/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgtCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGttt                       <  1:361422/249‑1 (MQ=255)
                                                                                                                                                                                                                        tCGCACTCGTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGttt                       <  2:283879/249‑1 (MQ=255)
                                                                                                                                                                                                                         cGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTgg                     <  1:267694/250‑1 (MQ=255)
                                                                                                                                                                                                                           taCTCCTGTTGCAGCATATCTGCGCGGTGCAAGCTGCCTTGCAGTCGTCCGCCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTGAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGCAAAGCAGTTCCACGTCACACCAGCTAAAGCGTTCATCCAGCAGATGAGAGACTGCCGGAGCCTGCCCTAATATCCATTCCCGGCTCCTCTGGCGGGCAAAGGTGTCGGCgtcgtttttgt                   >  2:21174/2‑240 (MQ=255)
                                                                                                                                                                                                                                 tGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTg             <  2:155507/250‑1 (MQ=255)
                                                                                                                                                                                                                                           tatCCGCGCGGTCCAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAAATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAATCTGGCGTTTTGtcc  >  2:376925/1‑249 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
CGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTC  >  NC_000913/4622964‑4623449

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