Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,339,749 A→G V106A (GTC→GCC)  mlaF ← intermembrane phospholipid transport system, ATP binding subunit MlaF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,339,7490AG97.7% 135.1 / ‑4.2 44V106A (GTC→GCC) mlaFintermembrane phospholipid transport system, ATP binding subunit MlaF
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/25);  total (18/26)
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.53e-01

CCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCT  >  NC_000913/3339626‑3339875
                                                                                                                           |                                                                                                                              
ccGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCaa                                                                                                                       >  2:338317/1‑133 (MQ=255)
  ggAAAGTTCTGAAGGTATTAGTTTCGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAAt                                                                                                                     <  2:272481/133‑1 (MQ=255)
  ggAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAAt                                                                                                                     <  2:375706/133‑1 (MQ=255)
    aaaGTTCAGAAGGCATTAGTTTAGGCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATAc                                                                                                                   <  2:279827/133‑1 (MQ=255)
     aaGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACg                                                                                                                  >  2:637886/1‑133 (MQ=255)
     aaGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTGGGGTATGTTCGCCCCGTTGATAGGCTGCGTGTTCAAATagg                                                                                                                  >  2:311233/1‑131 (MQ=255)
      aGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGt                                                                                                                 >  2:643211/1‑133 (MQ=255)
        ttCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTc                                                                                                               >  1:622344/1‑133 (MQ=255)
         tCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCa                                                                                                              <  2:407898/133‑1 (MQ=255)
                       tttAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGaacaac                                                                                                >  1:608454/1‑133 (MQ=255)
                               cTCCACGCCGCCCCAGGGACTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGa                                                                                        <  1:308807/133‑1 (MQ=255)
                               cTCCACGCAGCCCCACGGCCTCCACCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGa                                                                                        <  2:265860/133‑1 (MQ=255)
                                tCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGTCGTTGTCAGACAAGTTGATATCAGTGAACACCGCCCCGGAc                                                                                       >  2:169109/1‑133 (MQ=255)
                                    cGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATAACAGTGAACAACGCTGCGGACTGaa                                                                                   >  2:308807/1‑133 (MQ=255)
                                      cAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAAt                                                                                 >  2:48502/1‑133 (MQ=255)
                                             aCGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGc                                                                          <  1:48502/133‑1 (MQ=255)
                                                    ccAGCTTCATCATCAGCGTACTATACTACAATCGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCg                                                                   <  2:312824/133‑1 (MQ=255)
                                                         ttCATCATCACCGTAGTAGGCAACAATGGCTCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTgc                                                              <  1:222977/133‑1 (MQ=255)
                                                         ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTgc                                                              <  1:455897/133‑1 (MQ=255)
                                                          tcatcatcaTCGTACTATGCAACAATGGCGCGGGAAGTTGGTTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTgcg                                                             <  2:304131/133‑1 (MQ=255)
                                                          tcatcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTgcg                                                             <  1:548138/133‑1 (MQ=255)
                                                            atcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCa                                                           <  2:622344/133‑1 (MQ=255)
                                                            atcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCa                                                           <  1:338317/133‑1 (MQ=255)
                                                            atcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATATTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCa                                                           <  1:643211/133‑1 (MQ=255)
                                                            atcatcaCCGGACTATGCAACAATGGCGTGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCa                                                           <  2:274649/133‑1 (MQ=255)
                                                                           tGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGcgcg                                            >  1:378162/1‑133 (MQ=255)
                                                                                  aTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACgaga                                     >  2:132611/1‑133 (MQ=255)
                                                                                     gcgcGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACAt                                  >  2:289033/1‑133 (MQ=255)
                                                                                     gcgcGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACAt                                  >  1:439014/1‑133 (MQ=255)
                                                                                        cgGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGc                               >  1:448546/1‑133 (MQ=255)
                                                                                        cgGGAAGTTGGGTATGTTAGCGGAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGc                               <  1:311233/133‑1 (MQ=255)
                                                                                         gggAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGcc                              <  1:73525/133‑1 (MQ=255)
                                                                                          ggAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCg                                <  1:426507/130‑1 (MQ=255)
                                                                                          ggAAGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCg                                >  2:426507/1‑130 (MQ=255)
                                                                                            aaGTTGGGTATGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGa                           <  2:448546/133‑1 (MQ=255)
                                                                                                     aTGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGa                           >  2:77893/1‑124 (MQ=255)
                                                                                                     aTGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGa                           <  1:77893/124‑1 (MQ=255)
                                                                                                     aTGTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCa                  <  1:289033/133‑1 (MQ=255)
                                                                                                      ttttCGCGCTGTGGACAGGCGGCGTTGTAAAATACGTTCATGTCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCAc                 <  2:303605/131‑1 (MQ=255)
                                                                                                       gTTCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGTCGGAATGTTCTCAcc                >  1:303605/1‑133 (MQ=255)
                                                                                                         tCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCTGAATATTCTCACCAt              >  2:19344/1‑133 (MQ=255)
                                                                                                         tCGCGCAGTGGATAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCACCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCAt              <  1:526444/133‑1 (MQ=255)
                                                                                                                    aTAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAtc   <  2:439014/133‑1 (MQ=255)
                                                                                                                     tAGGCGGCGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAtct  <  2:608454/133‑1 (MQ=255)
                                                                                                                           |                                                                                                                              
CCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCT  >  NC_000913/3339626‑3339875

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