Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,742,575 G→A intergenic (+14/+26) cfa → / ← ribC cyclopropane fatty acyl phospholipid synthase, SAM‑dependent/riboflavin synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,742,5750GA97.6% 125.2 / ‑6.8 42intergenic (+14/+26)cfa/ribCcyclopropane fatty acyl phospholipid synthase, SAM‑dependent/riboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/20);  minor base T (1/0);  total (21/21)
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.95e-01

TGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGA  >  NC_000913/1742449‑1742679
                                                                                                                              |                                                                                                        
tGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCATGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgg                                                                                              >  2:257897/1‑139 (MQ=255)
    tACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCAGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGc                                                                                          <  1:621503/139‑1 (MQ=255)
         attatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgggg                                                                                            <  2:134895/132‑1 (MQ=255)
         attatCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCgggg                                                                                            >  1:134895/1‑132 (MQ=255)
            atCTGAATGCCTGTGCAGGTGCTTTCAGCGCCCGTGATATTCGGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAAAGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCTCACCCTCTCCGGGGGCGATTTCAg                                                                                  >  2:403283/1‑139 (MQ=255)
                 aaTGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAg                                                                             >  2:499074/1‑139 (MQ=255)
                    gCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCt                                                                          >  1:357632/1‑139 (MQ=255)
                     ccTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCATGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCtt                                                                         <  1:257897/139‑1 (MQ=255)
                       tgtgCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCt                                                                       >  2:610713/1‑139 (MQ=255)
                       tgtgCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCt                                                                       <  1:312390/139‑1 (MQ=255)
                           cAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGc                                                                   <  2:606645/139‑1 (MQ=255)
                              gTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTg                                                                >  1:288940/1‑139 (MQ=255)
                                gCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGt                                                              <  2:472778/139‑1 (MQ=255)
                                   ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGt                                                              <  1:872935/136‑1 (MQ=255)
                                   ttCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGt                                                              >  2:872935/1‑136 (MQ=255)
                                        cgcCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCAt                                                      >  2:593720/1‑139 (MQ=255)
                                        cgcCCGTGATATTCAGCTCTCGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCAt                                                      >  1:677769/1‑139 (MQ=255)
                                                     cAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGt                                         <  1:499074/139‑1 (MQ=255)
                                                       gCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTgc                                       <  2:357632/139‑1 (MQ=255)
                                                           tGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                      <  1:700117/72‑1 (MQ=255)
                                                           tGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCAcccc                                                                                                      >  2:700117/1‑72 (MQ=255)
                                                            ggCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCa                                  <  2:288940/139‑1 (MQ=255)
                                                              cAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGc                                <  1:286885/139‑1 (MQ=255)
                                                                ggTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCCCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGcac                              <  2:524168/139‑1 (MQ=255)
                                                                  tCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCggggg                                                                                           <  2:838226/76‑1 (MQ=255)
                                                                  tCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCggggg                                                                                           >  1:838226/1‑76 (MQ=255)
                                                                   cGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACg                           <  1:610713/139‑1 (MQ=255)
                                                                     tgtTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                         >  2:395207/1‑139 (MQ=255)
                                                                        tctcACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                      <  1:593720/139‑1 (MQ=255)
                                                                            aCGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGt                  >  1:205037/1‑139 (MQ=255)
                                                                               cgGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATc               >  1:867170/1‑139 (MQ=255)
                                                                                    ttGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCa          >  1:461688/1‑139 (MQ=255)
                                                                                        aaaCGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGc      <  1:395207/139‑1 (MQ=255)
                                                                                         aaCGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGcc     <  2:205037/139‑1 (MQ=255)
                                                                                            ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                         >  1:637025/1‑116 (MQ=255)
                                                                                            ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGtt                         <  2:637025/116‑1 (MQ=255)
                                                                                            ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                      <  2:373464/119‑1 (MQ=255)
                                                                                            ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTa                      >  1:373464/1‑119 (MQ=255)
                                                                                            ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGa  >  2:768603/1‑139 (MQ=255)
                                                                                            ggCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGa  >  1:791990/1‑139 (MQ=255)
                                                                                                                        tCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGAATCATGGCATTTTCTCGTGCCGCCAGc                                <  2:455869/81‑1 (MQ=255)
                                                                                                                        tCTATCACCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGAATCATGGCATTTTCTCGTGCCGCCAGc                                >  1:455869/1‑81 (MQ=255)
                                                                                                                              |                                                                                                        
TGTTTACCTATTATCTGAATGCCTGTGCAGGTGCTTTCCGCGCCCGTGATATTCAGCTCTGGCAGGTCGTGTTCTCACGCGGTGTTGAAAACGGCCTTCGAGTGGCTCGCTAAAGGCTATTCTATCGCCCCCTCTCCGGGGGCGATTTCAGATCAGGCTTCTGTGCCTGGTTGATTCATGGCATTTTCTCGTGCCGCCAGCACACGTTCTACCGTATCTACCACTGCCTGA  >  NC_000913/1742449‑1742679

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