Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,747,020 G→A A259V (GCG→GTG)  ffh ← Signal Recognition Particle (SRP) component with 4.5S RNA (ffs)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,747,0200GA97.4% 118.6 / ‑4.2 39A259V (GCG→GTG) ffhSignal Recognition Particle (SRP) component with 4.5S RNA (ffs)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/21);  total (16/23)
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.57e-01

ATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCGCCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAATCGAC  >  NC_000913/2746878‑2747166
                                                                                                                                              |                                                                                                                                                  
atCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCTCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGcaccgcaccg                                                                                                                                            >  2:521541/1‑151 (MQ=255)
atCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCGCCGCAcc                                                                                                                                             <  2:452897/150‑1 (MQ=255)
     cGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGATCGGCTCCAGCGCTTCTGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTGGCCAGTGATGTGACGAATAGAGAGCACCGCAccgccg                                                                                                                                         <  1:117390/149‑1 (MQ=255)
     cGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACcgccgc                                                                                                                                        <  2:328663/150‑1 (MQ=255)
         cAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGc                                                                                                                                    >  2:564297/1‑150 (MQ=255)
                aGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGt                                                                                                                            >  2:203808/1‑151 (MQ=255)
                aGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGt                                                                                                                            >  2:319618/1‑151 (MQ=255)
                             tGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCa                                                                                                                >  1:320270/1‑150 (MQ=255)
                                       acgcgacgcgATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCg                                                                                                     <  2:100645/151‑1 (MQ=255)
                                        cgcgacgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCgg                                                                                                    <  2:507473/151‑1 (MQ=255)
                                          cgacgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGt                                                                                                   >  2:170509/1‑150 (MQ=255)
                                           gacgcgaTGCGGTCCGGATGGAACGGCTCCAGCTCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCTCTTTGGTCAACACTACGCCGGTa                                                                                                  <  1:678873/150‑1 (MQ=255)
                                             cgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAg                                                                                                <  1:338005/150‑1 (MQ=255)
                                                          ggATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGATGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCTCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCtt                                                                                     <  1:321930/148‑1 (MQ=255)
                                                              ggAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGaa                                                                              <  2:559276/151‑1 (MQ=255)
                                                               gAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAAt                                                                             <  1:235868/151‑1 (MQ=255)
                                                                aaCGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCCCCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATg                                                                            <  1:563096/151‑1 (MQ=255)
                                                                     cTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCtt                                                                         <  1:649797/149‑1 (MQ=255)
                                                                         aGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCg                                                                   <  2:487185/151‑1 (MQ=255)
                                                                           cgcCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTa                                                                 <  2:596892/151‑1 (MQ=255)
                                                                           cgcCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGAGCCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCAGt                                                                  <  1:436798/150‑1 (MQ=255)
                                                                                cAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTg                                                              <  2:165888/149‑1 (MQ=255)
                                                                                      tctcGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCAt                                                       <  2:338721/150‑1 (MQ=255)
                                                                                       ctcGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATcc                                                     >  1:153604/1‑151 (MQ=255)
                                                                                         cgcaaccacCGAGAAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGCGCCCCGCACCGCCGCGGGCATCTCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTg                                                   <  1:70797/145‑1 (MQ=255)
                                                                                         cGCCAACACCGAGGACCTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCt                                                    <  2:186390/150‑1 (MQ=255)
                                                                                         cGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCt                                                    <  1:640686/150‑1 (MQ=255)
                                                                                                 ccGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTc                                            >  2:462200/1‑150 (MQ=255)
                                                                                                      gAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCg                                      >  1:605697/1‑151 (MQ=255)
                                                                                                      gAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCg                                      >  2:48653/1‑151 (MQ=255)
                                                                                                        aCTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTc                                    >  1:122086/1‑151 (MQ=255)
                                                                                                            gATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGTATCCTGACCGGTCATGGCGTCa                                   <  1:122995/148‑1 (MQ=255)
                                                                                                                ggTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACaa                            <  1:672545/151‑1 (MQ=255)
                                                                                                                    tGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAAc                          >  2:131294/1‑149 (MQ=255)
                                                                                                                           gATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTc                  <  1:126663/150‑1 (MQ=255)
                                                                                                                            aTGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCa                 >  2:588608/1‑150 (MQ=255)
                                                                                                                             tgtgACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAAc               >  1:77308/1‑151 (MQ=255)
                                                                                                                                  cGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGt          >  1:491336/1‑151 (MQ=255)
                                                                                                                                       agagagCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAATc     >  2:651834/1‑151 (MQ=255)
                                                                                                                                          gagCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAATCGAc  >  2:352751/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCGCCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAATCGAC  >  NC_000913/2746878‑2747166

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: