Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,052,261 G→A 9.4% intergenic (+21/+28) gnsA → / ← yccM putative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,052,2610GA 9.4% 162.3 / 7.5 53intergenic (+21/+28)gnsA/yccMputative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (2/3);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAG  >  NC_000913/1052121‑1052405
                                                                                                                                            |                                                                                                                                                
gcgcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTg                                                                                                                                           >  1:761605/1‑148 (MQ=255)
  gcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGAt                                                                                                                                         <  2:761605/148‑1 (MQ=255)
                           aaTACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGt                                                                                                                                       >  3:473498/1‑125 (MQ=255)
                           aaTACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGt                                                                                                                                       <  4:473498/125‑1 (MQ=255)
                                 gggAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGt                                                                                                                                            <  3:289325/114‑1 (MQ=255)
                                 gggAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGt                                                                                                                                            >  4:289325/1‑114 (MQ=255)
                                   gAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgttg                                                                                                                            >  2:1091665/1‑128 (MQ=255)
                                   gAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgttg                                                                                                                            <  1:1091665/128‑1 (MQ=255)
                                   gAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAt                                                                                                        >  4:900915/1‑148 (MQ=255)
                                           tctcTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                                                                        <  4:927997/108‑1 (MQ=255)
                                           tctcTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                                                                        >  3:927997/1‑108 (MQ=255)
                                              ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAg                                                                                                                                              >  1:665420/1‑99 (MQ=255)
                                              ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAg                                                                                                                                              <  2:665420/99‑1 (MQ=255)
                                              ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgt                                                                                                                              >  4:306408/1‑115 (MQ=255)
                                              ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgt                                                                                                                              <  3:306408/115‑1 (MQ=255)
                                                        cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCg                                                                                   >  2:684214/1‑148 (MQ=255)
                                                         ttCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                                                                 <  3:604708/117‑1 (MQ=255)
                                                         ttCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                                                                 >  4:604708/1‑117 (MQ=255)
                                                             ccGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCg                                                                                   <  4:447191/143‑1 (MQ=255)
                                                             ccGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCg                                                                                   >  3:447191/1‑143 (MQ=255)
                                                                cACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACa                                                                           <  1:684214/148‑1 (MQ=255)
                                                                 aCGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACaa                                                                          <  3:900915/148‑1 (MQ=255)
                                                                    agaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAATATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCa                                                                                                             <  2:633423/110‑1 (MQ=255)
                                                                    agaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAATATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCa                                                                                                             >  1:633423/1‑110 (MQ=255)
                                                                                  ggCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTa                                                                                       <  3:93506/118‑1 (MQ=255)
                                                                                  ggCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTa                                                                                       >  4:93506/1‑118 (MQ=255)
                                                                                          aGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCAt                                                 >  3:1749/1‑148 (MQ=255)
                                                                                             tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCt                                                                                                                                 <  3:409217/65‑1 (MQ=255)
                                                                                             tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCt                                                                                                                                 >  4:409217/1‑65 (MQ=255)
                                                                                               tgatgTCAAAATCAAAATAATGTGATTTTGTGCACATCACCCCGTGCGAGGTGATGTTCCGCttg                                                                                                                               <  1:320169/65‑1 (MQ=255)
                                                                                               tgatgTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttg                                                                                                                               >  2:320169/1‑65 (MQ=255)
                                                                                                gttcacaaaatcacattATTTTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGt                                                                              <  4:38615/108‑1 (MQ=38)
                                                                                                     caaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgt                                                                                                                              >  2:810891/1‑60 (MQ=255)
                                                                                                     caaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgt                                                                                                                              <  1:810891/60‑1 (MQ=255)
                                                                                                          tcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcc                                                                      <  3:303421/111‑1 (MQ=255)
                                                                                                          tcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcc                                                                      >  4:303421/1‑111 (MQ=255)
                                                                                                          tcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                 <  4:1749/148‑1 (MQ=255)
                                                                                                              aataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACaa                                                                          <  1:242130/103‑1 (MQ=255)
                                                                                                              aataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACaa                                                                          >  2:242130/1‑103 (MQ=255)
                                                                                                              aataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCa                                                               <  2:956754/114‑1 (MQ=255)
                                                                                                              aataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCa                                                               >  1:956754/1‑114 (MQ=255)
                                                                                                              aataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCt                                                        >  1:680179/1‑121 (MQ=255)
                                                                                                              aataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCt                                                        <  2:680179/121‑1 (MQ=255)
                                                                                                                             tgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACaa                                                                          <  1:239639/88‑1 (MQ=255)
                                                                                                                             tgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACaa                                                                          >  2:239639/1‑88 (MQ=255)
                                                                                                                               aaCATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                 >  1:698109/1‑127 (MQ=255)
                                                                                                                               aaCATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                 <  2:698109/127‑1 (MQ=255)
                                                                                                                                 cATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGc                                                                                    <  2:489552/74‑1 (MQ=255)
                                                                                                                                 cATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGc                                                                                    >  1:489552/1‑74 (MQ=255)
                                                                                                                                   tCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACAtt                                  <  3:881951/122‑1 (MQ=255)
                                                                                                                                   tCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAAt        >  3:633839/1‑148 (MQ=255)
                                                                                                                                       cccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATcccc    <  4:633839/148‑1 (MQ=255)
                                                                                                                                         cGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAg  >  2:662001/1‑148 (MQ=255)
                                                                                                                                            gCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAt                                                                                                        >  4:86143/1‑43 (MQ=255)
                                                                                                                                            gCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAt                                                                                                        <  3:86143/43‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAG  >  NC_000913/1052121‑1052405

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑