Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,547,395 A→G intergenic (‑226/‑6) yddG ← / → fdnG aromatic amino acid exporter/formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,547,3950AG92.5% 161.7 / 1.0 54intergenic (‑226/‑6)yddG/fdnGaromatic amino acid exporter/formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (23/26);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGAAGCAA  >  NC_000913/1547303‑1547490
                                                                                            |                                                                                               
gtgGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAAt                                                                                          >  2:517166/1‑100 (MQ=255)
 tgGCAAACAGATGTTCTTTTTTATTTCGTGCAAAAAGAATCAGAATTTTACTGTTAGTTTCCTCGCCCAGTTATACCCCTGAGTAATGTCGAAGGCAATg                                                                                         >  2:124345/1‑100 (MQ=255)
 tgGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATg                                                                                         >  2:330519/1‑100 (MQ=255)
  ggCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATgg                                                                                        >  2:1228895/1‑100 (MQ=255)
    cAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGAc                                                                                      <  2:898948/100‑1 (MQ=255)
     taaCAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACg                                                                                     <  1:25157/99‑1 (MQ=255)
      aaCAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACga                                                                                    >  2:25157/1‑99 (MQ=255)
            tGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCg                                                                              >  2:582729/1‑100 (MQ=255)
              ttCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCa                                                                            <  1:826237/100‑1 (MQ=255)
                cTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGa                                                                          >  1:647467/1‑100 (MQ=255)
                 tttttGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGAc                                                                         <  2:19841/100‑1 (MQ=255)
                       aTTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAtttt                                                                   >  1:710696/1‑100 (MQ=255)
                        tttCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGTGAAGCAATGGACGTCAGTCGCAGACAAttttt                                                                  <  1:1498817/100‑1 (MQ=255)
                         ttCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAtttttt                                                                 <  2:746929/100‑1 (MQ=255)
                         ttCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGGAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAtttttt                                                                 <  2:90847/100‑1 (MQ=255)
                           cgcgCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               <  1:529560/100‑1 (MQ=255)
                           cgcgCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               <  1:206206/100‑1 (MQ=255)
                           cgcgCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               >  2:365574/1‑100 (MQ=255)
                            gcggAAAAAGATTCAGAATGTTACTGTTAGTGTCCTCGCGCAGTAATACCCCTCAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAtttttt                                                                 <  2:1276158/97‑1 (MQ=255)
                            gcgcAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaaa                                                              <  1:1083507/100‑1 (MQ=255)
                             cgcAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaaaa                                                             >  1:720449/1‑100 (MQ=255)
                               cAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATc                                                           >  2:835156/1‑100 (MQ=255)
                               cAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATc                                                           >  2:1571304/1‑100 (MQ=255)
                                   aaGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGTAAGCAATGGATGTCAGTCGCAGACAATTTTTTAAAATCTgtg                                                       <  1:1066517/100‑3 (MQ=255)
                                   aaGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTgcg                                                       >  1:582490/1‑100 (MQ=255)
                                       ttCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcggg                                                   >  1:122421/1‑100 (MQ=255)
                                         cAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcgggcg                                                 <  1:1473420/100‑1 (MQ=255)
                                             aTTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTAt                                             <  2:763953/100‑1 (MQ=255)
                                             aTTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTAt                                             <  2:622532/100‑1 (MQ=255)
                                                taaCTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGGAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTAGGGc                                          >  1:1463918/3‑100 (MQ=255)
                                                   cTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTgg                                       <  2:51949/100‑1 (MQ=255)
                                                   cTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTgg                                       <  2:1463464/100‑1 (MQ=255)
                                                    tGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGa                                      <  1:822611/100‑1 (MQ=255)
                                                     gTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaa                                     >  1:246435/1‑100 (MQ=255)
                                                      ttAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaac                                    <  1:242582/100‑1 (MQ=255)
                                                         gTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaacaac                                 <  1:365574/100‑1 (MQ=255)
                                                          tttCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAacaaca                                <  1:321656/100‑1 (MQ=255)
                                                          tttCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAacaaca                                <  1:138700/100‑1 (MQ=255)
                                                                 gcgcAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCgg                         <  1:253226/100‑1 (MQ=255)
                                                                 gcgcAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCgg                         <  2:190644/100‑1 (MQ=255)
                                                                   gcAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCa                       <  1:491446/100‑1 (MQ=255)
                                                                    cAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCAt                      >  1:1297830/1‑100 (MQ=255)
                                                                      gTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTg                    >  1:986998/1‑100 (MQ=255)
                                                                       tAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTgg                   >  1:269137/1‑100 (MQ=255)
                                                                          tACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCt                <  2:828831/100‑1 (MQ=255)
                                                                          tACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCt                <  1:835161/100‑1 (MQ=255)
                                                                            ccccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCttt              >  2:1195602/1‑100 (MQ=255)
                                                                              ccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGc            >  2:190977/1‑100 (MQ=255)
                                                                               cTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcc           >  1:84278/1‑100 (MQ=255)
                                                                                tGAAAAAAGAGGGAAGCAACGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGccc          <  1:709655/100‑1 (MQ=255)
                                                                                 gAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcccc         <  1:1312508/100‑1 (MQ=255)
                                                                                  aaaaaaGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCg        >  1:1131888/1‑100 (MQ=255)
                                                                                       agagGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACGACAGTAGCGGCATTGGGCTTTGCCCCGaagca   <  1:124345/100‑1 (MQ=255)
                                                                                        gagGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGaagcaa  >  1:1453239/1‑100 (MQ=255)
                                                                                            |                                                                                               
GTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGAAGCAA  >  NC_000913/1547303‑1547490

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