Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,308,318 G→C R110G (CGT→GGT)  clsA ← cardiolipin synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,308,3180GC97.0% 112.2 / ‑6.8 33R110G (CGT→GGT) clsAcardiolipin synthase 1
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (25/7);  total (26/7)
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

TCCGCCATGCCGCCGGGCTGCCAGATATAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGCGCGGCGTTTGCCTAAATGG  >  NC_000913/1308170‑1308463
                                                                                                                                                    |                                                                                                                                                 
tCCGCCATGCCGCCGGGCTGCCAGATATAAAACCCCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTgaccac                                                                                                                                                 <  2:221513/151‑1 (MQ=255)
       tGCCGCCGGGCTGCCAGATATAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCa                                                                                                                                          >  2:47152/1‑151 (MQ=255)
              gggCTGCCAGATATAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTg                                                                                                                                   >  1:307488/1‑151 (MQ=255)
                        atatAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGGTCCCCTGACCACGCTCGCAAAGCTTGAATAATGgcg                                                                                                                         >  1:366030/1‑151 (MQ=255)
                        atatAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATgg                                                                                                                           <  2:131704/149‑1 (MQ=255)
                        atatAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGgcg                                                                                                                         >  1:345517/1‑151 (MQ=255)
                        atatAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGgcg                                                                                                                         >  1:149261/1‑151 (MQ=255)
                        aaataAAACACCATCTCAATATTTGGGCGCGCGATTTGGGTTTTGCGGATTAACACCTGCACCACACCACCTGACTCGGGTATCAGTTGTAGCTGATTTCCTTGGACTCCAGCGATTCCCTGACGCCGCTCCCAACGCTTaaaacatgcg                                                                                                                          >  1:156135/3‑143 (MQ=255)
                         tataAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGgcg                                                                                                                         >  2:125053/1‑150 (MQ=255)
                                 acCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGacac                                                                                                                <  2:172711/151‑1 (MQ=255)
                                     tctcAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTg                                                                                                              >  1:305738/1‑149 (MQ=255)
                                     tctcAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTg                                                                                                              >  2:63795/1‑149 (MQ=255)
                                      ctcAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTa                                                                                                           >  2:93840/1‑151 (MQ=255)
                                          atatTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTAtttt                                                                                                       >  2:113403/1‑151 (MQ=255)
                                          atatTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTAtttt                                                                                                       >  2:381712/1‑151 (MQ=255)
                                           tattatGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttc                                                                                                      <  1:356024/151‑1 (MQ=255)
                                           tattatGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttc                                                                                                      <  2:342315/151‑1 (MQ=255)
                                            attatGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTAtttt                                                                                                       >  1:374262/1‑149 (MQ=255)
                                                                  tCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAgg                                                                               >  2:354432/1‑151 (MQ=255)
                                                                       gATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTa                                                                          >  2:80207/1‑151 (MQ=255)
                                                                       gATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTa                                                                          >  2:351276/1‑151 (MQ=255)
                                                                        aTTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTAAAGGCTTTaa                                                                         >  2:291084/1‑151 (MQ=255)
                                                                         ttAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTCCAGGCTTTAAg                                                                        >  2:190036/1‑151 (MQ=255)
                                                                          tAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTGAAgg                                                                       <  1:6608/151‑1 (MQ=255)
                                                                                  gCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAg                                                               >  2:132105/1‑151 (MQ=255)
                                                                                  gCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAg                                                               >  2:343587/1‑151 (MQ=255)
                                                                                         atcatcTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCAtttt                                                        <  1:78977/151‑1 (MQ=255)
                                                                                           catcTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTg                                                       >  1:327757/1‑150 (MQ=255)
                                                                                           catcTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTg                                                       >  1:359252/1‑150 (MQ=255)
                                                                                                             gTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCg                                     >  1:267759/1‑150 (MQ=255)
                                                                                                                                      aGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGCGCGGCGTTTg           >  2:63605/1‑151 (MQ=255)
                                                                                                                                       gCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGAAACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGCGCGGCGTTTg           >  2:167561/1‑150 (MQ=255)
                                                                                                                                             cccTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGCGCGGCGTTTGCCTAAAt    >  1:76657/1‑151 (MQ=255)
                                                                                                                                               cTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGCGCGGCGTTTGCCTAAATgg  <  1:221940/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCCGCCATGCCGCCGGGCTGCCAGATATAAAACACCATCTCAATATTATGGCGCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCGGTGGAAGGCCACATCGCTCTGGCGCGCTCAGCGCGGCGTTTGCCTAAATGG  >  NC_000913/1308170‑1308463

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