Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,972,909 (C)6→7 coding (954/1194 nt) lplT ← lysophospholipid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,972,9031.C95.8% 82.3 / ‑2.6 24coding (960/1194 nt)lplTlysophospholipid transporter
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (12/11);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GCGAGTAAATGCCCAGCATCAACAACATGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCAAAATCCCGGCTGGCATACAGCGTGACACGGTTTCCAGC  >  NC_000913/2972773‑2973029
                                                                                                                                   |                                                                                                                              
gCGAGTAAATGCCCAGCATCAACAACATGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGccccccca                                                                                                                        <  2:712957/139‑3 (MQ=255)
 cGAGTAAATGCCCAGCATCAACAACATGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGcccccccat                                                                                                                       >  2:566696/1‑136 (MQ=255)
        aTGCCCAGCATCAACAACATGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCa                                                                                                                 <  2:678348/139‑1 (MQ=255)
                cATCAACAACATGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCAt                                                                                                         <  2:116469/139‑1 (MQ=255)
                          aTGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCa                                                                                               >  2:52954/1‑139 (MQ=255)
                              cgcTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCaa                                                                                                                >  1:639154/1‑118 (MQ=255)
                              cgcTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCaa                                                                                                                <  2:639154/118‑1 (MQ=255)
                                             ggTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCgt                                                                            >  1:529611/1‑139 (MQ=255)
                                             ggTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCgt                                                                            >  1:384803/1‑139 (MQ=255)
                                                   gTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCa                                                                      <  1:52954/139‑1 (MQ=255)
                                                        gCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGTGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCACGTGGCGCAGGGaa                                                                 >  2:486886/1‑139 (MQ=255)
                                                         cAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGaaa                                                                >  2:616665/1‑139 (MQ=255)
                                                                      cgccccgACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTaccac                                                   >  1:69918/1‑139 (MQ=255)
                                                                          ccgACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACg                                               <  1:566696/139‑1 (MQ=255)
                                                                             aCGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCTACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCa                                            >  2:433654/1‑139 (MQ=255)
                                                                                          cccGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCAAAATCCCgg                               <  2:529611/139‑1 (MQ=255)
                                                                                                   gTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCAAAATCCCGGCTGGCATACa                     >  1:41752/1‑139 (MQ=255)
                                                                                                         aCGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAg                                                                                                            <  2:211441/47‑1 (MQ=255)
                                                                                                         aCGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAg                                                                                                            >  1:211441/1‑47 (MQ=255)
                                                                                                          cGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCa                                        >  1:663583/1‑114 (MQ=255)
                                                                                                          cGCATTGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCa                                        <  2:663583/114‑1 (MQ=255)
                                                                                                              ttGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTg                                                                        >  1:111863/1‑78 (MQ=255)
                                                                                                              ttGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTg                                                                        <  2:111863/78‑1 (MQ=255)
                                                                                                               tGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATcaaca                                                                                                    >  1:516250/1‑49 (MQ=255)
                                                                                                               tGAGCGGAACGACAAAAAAGCCCCCCCATCACGCCAATCAGCATcaaca                                                                                                    <  2:516250/49‑1 (MQ=38)
                                                                                                                       aCGACAAAAAAGCCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCAAAATCCCGGCTGGCATACAGCGTGACACGGTTTCCAGc  <  1:616665/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                              
GCGAGTAAATGCCCAGCATCAACAACATGGCGCTGTTTTCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCTGGCAGCAGCTCGTGTTGCAGGGAAAAAATCAGTACCACCACGCCAATCAAAATCCCGGCTGGCATACAGCGTGACACGGTTTCCAGC  >  NC_000913/2972773‑2973029

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