Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,860 A→G 10.5% intergenic (‑102/+13) 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,8600AG10.5% 134.2 / 7.3 49intergenic (‑102/+13)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (2/3);  total (24/25)
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.23e-01

CATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTG  >  NC_000913/1987716‑1988002
                                                                                                                                                 |                                                                                                                                              
caTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAg                                                                                                                                            <  4:214452/149‑1 (MQ=255)
   cAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGtt                                                                                                                                         <  2:337660/149‑1 (MQ=255)
      aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGGtttt                                                                                                                                       >  4:237749/1‑148 (MQ=255)
      aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGGtttt                                                                                                                                       <  3:237749/148‑1 (MQ=255)
       aaCCAATTCCTGATCGGCATTCAGCGGCTAGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAATAGAACCTGC‑CGGAGCAGGtttttt                                                                                                                                     <  2:24654/149‑1 (MQ=255)
          cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCATCGTCATCAGGACACAAAAAAACGTGCTCCG‑GCAGGttttttt                                                                                                                                    <  3:364963/147‑2 (MQ=255)
          cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGGttttttt                                                                                                                                    >  4:364963/1‑146 (MQ=255)
                cTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGa                                                                                                                            >  1:521033/1‑149 (MQ=255)
                         aTTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGc                                                                                                                   >  1:51074/1‑149 (MQ=255)
                             aGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTgg                                                                                                               >  2:27719/1‑149 (MQ=255)
                                 gCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGt                                                                                                           >  4:456523/1‑149 (MQ=255)
                                 gCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGt                                                                                                           <  1:214338/149‑1 (MQ=255)
                                      gACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTc                                                                                                      >  2:3391/1‑149 (MQ=255)
                                             tAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                               >  2:280930/1‑149 (MQ=255)
                                                 cTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATAtt                                                                                                                     >  1:95998/1‑123 (MQ=255)
                                                 cTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATAtt                                                                                                                     <  2:95998/123‑1 (MQ=255)
                                                  tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCG‑GCAGGttttttt                                                                                                                                    >  4:214413/1‑105 (MQ=38)
                                                  tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGC‑TCCGGCAGGttttttt                                                                                                                                    <  3:214413/106‑2 (MQ=38)
                                                    cTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                        <  1:27719/149‑1 (MQ=255)
                                                        tAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                    <  1:280930/149‑1 (MQ=255)
                                                            aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                               <  2:268200/134‑1 (MQ=255)
                                                            aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                               >  1:268200/1‑134 (MQ=255)
                                                            aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTa                                                                                >  4:23781/1‑149 (MQ=255)
                                                                          tataGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACtttt                                                                  <  3:390693/149‑1 (MQ=255)
                                                                            tagtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                <  1:3391/149‑1 (MQ=255)
                                                                               tagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                             >  4:403039/1‑149 (MQ=255)
                                                                                         acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                               <  2:489428/105‑1 (MQ=255)
                                                                                         acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                               >  1:489428/1‑105 (MQ=255)
                                                                                           aTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                 >  2:162057/1‑149 (MQ=255)
                                                                                            tCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGttt                                                >  3:568061/1‑149 (MQ=255)
                                                                                            tCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGttt                                                >  4:508053/1‑149 (MQ=255)
                                                                                               tCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCAt                                             <  2:51074/149‑1 (MQ=255)
                                                                                                     aTAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGt                                                  <  2:203097/138‑1 (MQ=255)
                                                                                                     aTAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGAGCAGt                                                  >  1:203097/1‑138 (MQ=255)
                                                                                                             aTCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                               >  3:476542/1‑149 (MQ=255)
                                                                                                                  ccGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAttt                          <  1:162057/149‑1 (MQ=255)
                                                                                                                   cGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTg                         <  3:403039/149‑1 (MQ=255)
                                                                                                                        tcaGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTtata                    <  3:456523/149‑1 (MQ=255)
                                                                                                                        tcaGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTtata                    <  4:476542/149‑1 (MQ=255)
                                                                                                                               acaAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTc                                                                                            <  4:396088/70‑1 (MQ=255)
                                                                                                                               acaAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTc                                                                                            >  3:396088/1‑70 (MQ=255)
                                                                                                                                 aaaaaaaCCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt               <  2:520956/145‑1 (MQ=255)
                                                                                                                                 aaaaaaaCCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt               >  1:520956/1‑145 (MQ=255)
                                                                                                                                  aaaaaaCCTGC‑TCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                             <  2:31974/98‑1 (MQ=39)
                                                                                                                                  aaaaaaCCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt          <  2:521033/149‑1 (MQ=255)
                                                                                                                                    aaaaCCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg        >  1:449161/1‑149 (MQ=255)
                                                                                                                                        ccTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCtt    >  1:589517/1‑149 (MQ=255)
                                                                                                                                         cTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCttt   <  3:508053/149‑1 (MQ=255)
                                                                                                                                          tGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTg  >  4:230065/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTG  >  NC_000913/1987716‑1988002

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: ‑