Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,656 T→A 11.8% intergenic (‑33/+151) 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,6560TA11.8% 86.2 / 5.9 34intergenic (‑33/+151)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 = 7.88e-02

CAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGT  >  NC_000913/3101524‑3101798
                                                                                                                                    |                                                                                                                                              
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)
                                                                                                                                    ttGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt  >  1:127226/1‑143 (MQ=255)
                                                                                                                                    |                                                                                                                                              
CAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGT  >  NC_000913/3101524‑3101798

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: