Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,632,303 (T)6→5 intergenic (‑18/‑587) ydfJ ← / → ynfT putative transporter YdfJ/Qin prophage; protein YnfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,632,2980T.95.7% 171.2 / 2.7 46intergenic (‑13/‑592)ydfJ/ynfTputative transporter YdfJ/Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (22/22);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TACGCGCCACATAACCTGCACCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACT  >  NC_000913/1632159‑1632433
                                                                                                                                           |                                                                                                                                       
tACGTGCCACATAACCTGCACCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGgtgt                                                                                                                                         <  2:726643/140‑1 (MQ=255)
     gcCACATAGCCTGCACTTTAGGCTCCCATTTCCTGTATTAACGCCATTTCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGtttttc                                                                                                                                    <  2:401086/140‑2 (MQ=255)
     gcCACATAACCTGCACCGTAGGTTCCCATTGTCAGCATTGACGCCATTTCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGtttttc                                                                                                                                    <  2:315469/140‑2 (MQ=255)
       cacaTAACCTGCACCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTccc                                                                                                                                  >  2:826934/1‑140 (MQ=255)
       cacaTAACCTGCACCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTccc                                                                                                                                  >  2:650755/1‑140 (MQ=255)
                 gCACCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGt                                                                                                                       >  2:515322/1‑140 (MQ=255)
                   aCCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGttt                                                                                                                     >  2:735046/1‑140 (MQ=255)
                   aCCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGttt                                                                                                                     >  1:12504/1‑140 (MQ=255)
                       tAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCtt                                                                                                                 >  1:673271/1‑140 (MQ=255)
                       tAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCtt                                                                                                                 >  1:739767/1‑140 (MQ=255)
                       tAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCtt                                                                                                                 >  1:432271/1‑140 (MQ=255)
                        aGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTg                                                                                                                >  2:539165/1‑140 (MQ=255)
                         ggTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTgg                                                                                                               <  1:735046/140‑1 (MQ=255)
                               cATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAAACCATAATTAACAAGTGTTTTTCCCAACCTGAGGTTTCCTTGGCGtttt                                                                                                         >  2:379461/1‑140 (MQ=255)
                                     cAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGt                                                                                                   >  1:448020/1‑140 (MQ=255)
                                        aaTTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGAGTTTCTCCCATCATGTAGTCTCCTTGGCGTTTTCTAGGtttt                                                                                                >  1:393623/1‑140 (MQ=255)
                                             aCGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGATTTTTTCAg                                                                                           >  2:502225/1‑140 (MQ=255)
                                                  aTTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGt                                                                                      >  2:795897/1‑140 (MQ=255)
                                                    tGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTg                                                                                    >  2:165433/1‑140 (MQ=255)
                                                        gTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCAtt                                                                                >  2:460879/1‑140 (MQ=255)
                                                               atTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaa                                                                         <  2:448020/140‑1 (MQ=255)
                                                                ttCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaa                                                                        <  1:826934/140‑1 (MQ=255)
                                                                ttCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaa                                                                        >  1:377049/1‑140 (MQ=255)
                                                                 tCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaaa                                                                       <  1:795897/140‑1 (MQ=255)
                                                                  cAGGAAAAAATATTTCTTGAAAGACTAACGCCGCGTCGAGCGAATATATCTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTAGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaaaa                                                                      <  2:414545/140‑1 (MQ=255)
                                                                  cAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaaaa                                                                      >  2:236171/1‑140 (MQ=255)
                                                                   aGGAAAAAATATTTCATGAAACTCTGACGGTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAg                                                                     <  1:379461/140‑1 (MQ=255)
                                                                   aGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAg                                                                     <  1:650755/140‑1 (MQ=255)
                                                                   aGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAg                                                                     <  1:142787/140‑1 (MQ=255)
                                                                   aGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAg                                                                     <  2:12504/140‑1 (MQ=255)
                                                                        aaaaTATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATcc                                                                <  2:673271/140‑1 (MQ=255)
                                                                          aaTATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTa                                                              <  1:165433/140‑1 (MQ=255)
                                                                            tatTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAg                                                            <  1:502225/140‑1 (MQ=255)
                                                                                 cATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaaaa                                                                      <  1:49601/125‑1 (MQ=255)
                                                                                 cATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTaaaaa                                                                      >  2:49601/1‑125 (MQ=255)
                                                                                  aTGAAACACTAACGCTGTGCCGCGCGAATATAACTTGAAGTCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTAGAt                                                      <  2:322511/140‑1 (MQ=255)
                                                                                    gAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGAtct                                                    <  2:504127/140‑1 (MQ=255)
                                                                                    gAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGAtct                                                    <  1:167173/140‑1 (MQ=255)
                                                                                         acTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAgtg                                               >  2:458274/1‑140 (MQ=255)
                                                                                           tAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTc                                             <  1:539165/140‑1 (MQ=255)
                                                                                                   cgcCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCa                                                  <  1:475064/127‑1 (MQ=255)
                                                                                                   cgcCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCa                                                  >  2:475064/1‑127 (MQ=255)
                                                                                                    gcCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTgggg                                    <  1:460344/140‑1 (MQ=255)
                                                                                                     ccGAGCGAATATAACTGAAAATCCATAATTCCCAGTTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGc                                   <  2:377049/140‑1 (MQ=255)
                                                                                                      cGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGcc                                  <  1:730514/140‑1 (MQ=255)
                                                                                                            aaTATAACTGGAAATCCATAATTCACAGGTGTTTTTCTCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTc                            <  1:236171/140‑1 (MQ=255)
                                                                                                                      gAAATCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCatat                  >  2:535062/1‑140 (MQ=255)
                                                                                                                       aaaTCCATAATTCACAGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCatata                 >  1:541143/1‑140 (MQ=255)
                                                                                                                                      aGGTGTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACt  >  2:441700/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TACGCGCCACATAACCTGCACCGTAGGTTCCCATTGCCAGAATTAACGCCATTGCCGTTGATGATTCAGGAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACT  >  NC_000913/1632159‑1632433

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