Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,858 G→C 7.8% intergenic (‑100/+15) yecJ ← / ← azuC DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,8580GC 7.8% 176.9 / 5.3 51intergenic (‑100/+15)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base C (2/2);  total (25/26)
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.87e-01

CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATT  >  NC_000913/1987713‑1987995
                                                                                                                                                 |                                                                                                                                          
ccACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAg                                                                                                                                         >  4:372635/1‑149 (MQ=255)
      cAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtt                                                                                                                                  <  3:372635/149‑1 (MQ=255)
                ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATc                                                                                                                        <  4:843362/149‑1 (MQ=255)
                       ttggCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAACCTGCTCCGGCAGGttt                                                                                                                                 <  1:1072056/130‑1 (MQ=255)
                       ttggCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAACCTGCTCCGGCAGGttt                                                                                                                                 >  2:1072056/3‑132 (MQ=255)
                               cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGc                                                                                                                                        >  3:1072130/1‑119 (MQ=255)
                               cAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGc                                                                                                                                        <  4:1072130/119‑1 (MQ=255)
                                aGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTgg                                                                                                        >  3:1082144/1‑149 (MQ=255)
                                 gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTg                                                                                                             >  3:1453885/1‑143 (MQ=255)
                                 gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTg                                                                                                             <  4:1453885/143‑1 (MQ=255)
                                           cATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                             <  4:1082144/149‑1 (MQ=255)
                                                          aTAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                        >  4:924672/1‑139 (MQ=255)
                                                          aTAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTAGCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                        <  3:924672/139‑1 (MQ=255)
                                                                ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtttttt                                                                                                                              <  1:81011/95‑1 (MQ=255)
                                                                ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGtttttt                                                                                                                              >  2:81011/1‑95 (MQ=255)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTg                                                                                                                             >  3:153750/1‑85 (MQ=255)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTg                                                                                                                             <  4:153750/85‑1 (MQ=255)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGt                                                                                                    >  2:549763/1‑110 (MQ=255)
                                                                           gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGt                                                                                                    <  1:549763/110‑1 (MQ=255)
                                                                            ttATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCg                                                                                                                       <  2:427060/90‑1 (MQ=255)
                                                                            ttATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCg                                                                                                                       >  1:427060/1‑90 (MQ=255)
                                                                                            acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCaa                                                                               <  4:762737/114‑1 (MQ=255)
                                                                                            acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCaa                                                                               >  3:762737/1‑114 (MQ=255)
                                                                                                     gCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGcc                                                                                                           <  3:1331476/77‑1 (MQ=255)
                                                                                                     gCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGcc                                                                                                           >  4:1331476/1‑77 (MQ=255)
                                                                                                      cGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAg                                                                            <  4:1142504/107‑1 (MQ=255)
                                                                                                      cGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAg                                                                            >  3:1142504/1‑107 (MQ=255)
                                                                                                            tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATc                                                                                                                        <  2:566922/57‑1 (MQ=255)
                                                                                                            tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATc                                                                                                                        >  1:566922/1‑57 (MQ=255)
                                                                                                            tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                  >  4:1281544/1‑79 (MQ=255)
                                                                                                            tCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                  <  3:1281544/79‑1 (MQ=255)
                                                                                                              acaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                             >  2:935423/1‑98 (MQ=255)
                                                                                                              acaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                             <  1:935423/98‑1 (MQ=255)
                                                                                                                          atcaGTACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                            <  1:1384186/71‑1 (MQ=255)
                                                                                                                          atcaGTACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                            >  2:1384186/1‑71 (MQ=255)
                                                                                                                          atcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGt                                                                                                >  1:170363/1‑67 (MQ=255)
                                                                                                                          atcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGt                                                                                                <  2:170363/67‑1 (MQ=255)
                                                                                                                           tcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtttt                                                  <  3:142179/112‑1 (MQ=255)
                                                                                                                           tcaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtttt                                                  >  4:142179/1‑112 (MQ=255)
                                                                                                                             aGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTg                                                                                  >  3:1053696/1‑78 (MQ=255)
                                                                                                                             aGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTg                                                                                  <  4:1053696/78‑1 (MQ=255)
                                                                                                                                  acaAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                             >  4:93580/1‑78 (MQ=255)
                                                                                                                                  acaAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                             <  3:93580/78‑1 (MQ=255)
                                                                                                                                    aaaaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa    >  1:256935/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCTGCTCCGGCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                         <  3:323950/64‑1 (MQ=18)
                                                                                                                                     aaaaaaCCTGCTCCGGCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTGTGaa                                                                                         >  4:323950/1‑64 (MQ=17)
                                                                                                                                      aaaaaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt  <  2:256935/149‑1 (MQ=255)
                                                                                                                                         aaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTc                                                          <  2:1334553/90‑1 (MQ=255)
                                                                                                                                         aaCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTc                                                          >  1:1334553/1‑90 (MQ=255)
                                                                                                                                             tGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtt                                                    >  3:471477/1‑92 (MQ=255)
                                                                                                                                             tGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtt                                                    <  4:471477/92‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATT  >  NC_000913/1987713‑1987995

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