Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,618,274 C→T G212D (GGC→GAC)  marC ← UPF0056 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,618,2740CT97.9% 174.5 / ‑5.8 47G212D (GGC→GAC) marCUPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/18);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCAC  >  NC_000913/1618129‑1618410
                                                                                                                                                 |                                                                                                                                        
gtgGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGtt                                                                                                                                      <  2:229704/150‑1 (MQ=255)
  gggCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTtaat                                                                                                                                   >  2:771926/1‑151 (MQ=255)
      gcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTtaataata                                                                                                                               <  2:48858/151‑1 (MQ=255)
        ggTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTtaataataa                                                                                                                              >  2:170934/1‑150 (MQ=255)
         gTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATaaa                                                                                                                             >  2:563092/1‑150 (MQ=255)
              tGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGcac                                                                                                                       >  1:134906/1‑151 (MQ=255)
               gCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTTGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCC‑TTTTTAATGGTACGTTTTAATGATGTCCAGGATGTCGTTAATAATAAACTGcac                                                                                                                       <  1:436686/149‑1 (MQ=255)
                  tttGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACAccc                                                                                                                   >  1:467855/1‑151 (MQ=255)
                    tGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCcat                                                                                                                 >  2:208526/1‑151 (MQ=255)
                     gCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCcat                                                                                                                 >  1:638285/1‑150 (MQ=255)
                        gcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCcatacat                                                                                                             <  1:380251/151‑1 (MQ=255)
                         cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACTCCcatacata                                                                                                            <  1:386739/151‑1 (MQ=255)
                             aaTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACca                                                                                                         <  2:481115/150‑1 (MQ=255)
                                ttGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACcagc                                                                                                       >  1:318311/1‑149 (MQ=255)
                                      aaTCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATcc                                                                                               <  2:617707/151‑1 (MQ=255)
                                          aTTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                           <  1:466132/151‑1 (MQ=255)
                                           ttATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                           >  2:41938/1‑150 (MQ=255)
                                           ttATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                           >  2:153203/1‑150 (MQ=255)
                                            tataTTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATc                                                                                           <  2:716078/149‑1 (MQ=255)
                                                tttCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACgg                                                                                     >  2:646689/1‑151 (MQ=255)
                                                   cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGgag                                                                                  <  2:73568/151‑1 (MQ=255)
                                                    gccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACggg                                                                                    >  1:617384/1‑148 (MQ=255)
                                                         tGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGc                                                                             >  1:138636/1‑150 (MQ=255)
                                                          ggCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCtt                                                                           >  1:413492/1‑151 (MQ=255)
                                                              aGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAAt                                                                       >  2:624828/1‑151 (MQ=255)
                                                                  acacTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGc                                                                     >  2:560099/1‑149 (MQ=255)
                                                                  acacTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCa                                                                   <  2:388164/151‑1 (MQ=255)
                                                                         aagaaCAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGc                                                            <  1:508698/151‑1 (MQ=255)
                                                                                gACGCAATAGTTGAAAGGACCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTGCCAGGATGTCGTTAGTAATAAACTGCACGCCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTGAATGCCACCCTTGCcctcca                                                      <  1:673725/150‑1 (MQ=255)
                                                                                     aaTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGc                                                 >  2:420200/1‑150 (MQ=255)
                                                                                      aTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                               <  1:132050/151‑1 (MQ=255)
                                                                                         gTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAAt                                            >  2:504903/1‑151 (MQ=255)
                                                                                                         gggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGt                            >  1:714098/1‑151 (MQ=255)
                                                                                                           gCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAg                          <  1:321903/151‑1 (MQ=255)
                                                                                                             cTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGc                        >  1:229650/1‑151 (MQ=255)
                                                                                                              ttttttAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCt                       >  2:469271/1‑151 (MQ=255)
                                                                                                               tttttAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCtt                      >  1:645197/1‑151 (MQ=255)
                                                                                                               tttttAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCtt                      >  2:268258/1‑151 (MQ=255)
                                                                                                                ttttAATGGTACGTTGTAATGATTTCCAGGATGCCATTAATAATAAACTGCTCCCCCATACAAACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCTCCCTTGCCCCCCACCCGCATAATTGCGCCGGAGCTGCGTAGGCt                       <  1:152551/149‑1 (MQ=255)
                                                                                                                 tttAATGGTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTcc                    >  1:348169/1‑151 (MQ=255)
                                                                                                                       ggTACGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaa              <  2:374847/151‑1 (MQ=255)
                                                                                                                          aCGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaaaa            <  2:318932/150‑1 (MQ=255)
                                                                                                                           cGTTTTAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAAt           >  1:338090/1‑150 (MQ=255)
                                                                                                                               ttAATGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCg      >  1:9638/1‑151 (MQ=255)
                                                                                                                                 aaTGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcc    <  2:265704/151‑1 (MQ=255)
                                                                                                                                   tGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGccac  >  2:696784/1‑151 (MQ=255)
                                                                                                                                   tGATTTCCAGGATGTCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGccac  >  1:121814/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
GTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCAC  >  NC_000913/1618129‑1618410

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