Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,655 T→A 11.8% intergenic (‑32/+152) yggN ← / ← yggL putative EcfF/protein YggL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,6550TA11.8% 82.8 / 5.9 34intergenic (‑32/+152)yggN/yggLputative EcfF/protein YggL
Reads supporting (aligned to +/- strand):  ref base T (17/13);  new base A (2/2);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

CAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCG  >  NC_000913/3101524‑3101793
                                                                                                                                   |                                                                                                                                          
cAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                                                                                                                                 <  1:468778/143‑3 (MQ=255)
   tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAgc                                                                                                                              >  1:230918/1‑142 (MQ=255)
   tCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACGGAGAGTGGCGGCGCCATCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                                                                                              >  1:575775/1‑143 (MQ=255)
   tCACACCGTCGCGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgatt                                                                                                                              <  1:558776/143‑6 (MQ=255)
       gtCGTCACGCGGCGTGGCGCTGCACTGGGAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCAACCAACAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                                                                                                          <  2:48844/142‑1 (MQ=255)
        tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                                                                                                         <  1:121899/143‑1 (MQ=255)
                              aCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATg                                                                                                   >  2:546127/1‑143 (MQ=255)
                                 ggTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCg                                                                                                >  1:305183/1‑143 (MQ=255)
                                 ggTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCg                                                                                                >  1:309007/1‑143 (MQ=255)
                                     gTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGAGCAGCCGGCAGCAGGATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                                                                                            >  1:397113/1‑143 (MQ=255)
                                      tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGCCCAGCATTTGCTGGCCCCAGAATGATACCAAAGTACGCGTTTTCCATGCCGGGTgc                                                                                           >  1:479004/1‑143 (MQ=255)
                                        ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcgg                                                                                         <  1:341690/143‑1 (MQ=255)
                                        ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcgg                                                                                         >  1:555638/1‑143 (MQ=255)
                                         gCGTGAGCGGTCGTTGCCGTCACTGAAAGTGCTGCCGCTAGCAGCATTTTGCGCATCATAGTCTTCCCTCATGAAAAAATCAGGCCAGCCTTCGCGGGCCCCAGATTGATAACAAAGTACGCGTTTCCCATGCCGGATgcggc                                                                                        <  2:565959/143‑1 (MQ=255)
                                          cGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTCGCTGGCCCCAGATTGATATCAAAGTACGCGGTTTCCATGCCGGATgcggcg                                                                                       <  2:302412/143‑1 (MQ=255)
                                                 ggTCATGGCCGTGACGGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGGTTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac                                                                                <  2:393231/143‑1 (MQ=255)
                                                       tGCTGTCACTGAACGTGCTGGCGCTAGCAGCATCTTGCGCTTCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACACAGTACGCGTCATCCGGGGCTGATTCGGCGTAAACACCTTATc                                                                          <  2:352007/143‑1 (MQ=255)
                                                         ccGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATccc                                                                        >  1:548606/1‑143 (MQ=255)
                                                            tCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAACTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACACAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGcc                                                                     >  2:96054/1‑143 (MQ=255)
                                                            tCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTTCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGGTGCGGCGTAAACAACTTACCCCGcc                                                                     >  2:145301/1‑143 (MQ=255)
                                                                 gAAAGTGCTGCGGTCAGCAGCCTTTGGCGCATCGTAGTCTAAAAAAAAAAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTc                                                                <  2:397577/143‑1 (MQ=37)
                                                                    atggctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCATATTGATAACCAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGAt                                                             <  2:567933/140‑1 (MQ=255)
                                                                    aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGAt                                                             <  1:546127/143‑1 (MQ=255)
                                                                       gctgcGACCAGCAGCATTTTGCGTATCATATTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGt                                                          <  2:179200/143‑1 (MQ=255)
                                                                               cagcagCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCttt                                                  >  1:326943/1‑143 (MQ=255)
                                                                                    gcaTTTTGCGTATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAgg                                             <  2:242988/143‑1 (MQ=255)
                                                                                            tgcATCTTAGTCTACCCTTACGCAAAAAGCGGGTCGGCATTTGCTGGCCCCAGATTGATTACAAAGTCCGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATaa                                     <  2:284131/142‑1 (MQ=37)
                                                                                              catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGc                                   >  1:548397/1‑143 (MQ=255)
                                                                                              catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGAGAAGc                                   >  1:61246/1‑143 (MQ=255)
                                                                                                tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGt                                 >  1:527865/1‑143 (MQ=255)
                                                                                                             tCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTACACACCTTATCCCGGCTACGCGATCGTGCATCCTGCCAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  2:278106/1‑143 (MQ=255)
                                                                                                                     atctggggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGAGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAAGCATGATAATCGTAGCGCATCAGAAAGTTTTGCa            >  2:295351/7‑143 (MQ=255)
                                                                                                                            gCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCACTTGTCa     >  1:75016/1‑143 (MQ=255)
                                                                                                                               aGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACCCCTTATCCCGCCTACTCGATCGTGCCTCTTTCCAGGCCTGCTAAGCGTAGCGCCGCAGGCAGTTTTGCATTTGTCATCg  >  1:90710/1‑143 (MQ=39)
                                                                                                                                   |                                                                                                                                          
CAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCG  >  NC_000913/3101524‑3101793

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: