Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,167,184 A→G noncoding (71/87 nt) cyaR → small regulatory RNA CyaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,167,1840AG97.9% 160.7 / ‑4.4 48noncoding (71/87 nt)cyaRsmall regulatory RNA CyaR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (27/20);  total (28/20)
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.99e-01

ATGGAAAATTCTTAGAAACCGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGC  >  NC_000913/2167053‑2167320
                                                                                                                                   |                                                                                                                                        
aTGGAAAATTCTTAGAAACCGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCt                                                                                                                                  >  1:899731/1‑140 (MQ=255)
           ttAGAAACCGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTAcg                                                                                                                       >  2:528984/1‑140 (MQ=255)
             aGAAACCGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTAcgcg                                                                                                                     <  1:989853/140‑1 (MQ=255)
                aaCCGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTaa                                                                                                                  <  1:95248/140‑1 (MQ=255)
                   cGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTAcgcg                                                                                                                     >  2:196104/1‑134 (MQ=255)
                   cGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTAcgcg                                                                                                                     <  1:196104/134‑1 (MQ=255)
                                cTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACaa                                                                                                             <  2:483265/129‑1 (MQ=255)
                                cTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACaa                                                                                                             >  1:483265/1‑129 (MQ=255)
                                cTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAg                                                                                                  >  2:152930/1‑140 (MQ=255)
                                cTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAg                                                                                                  >  1:961277/1‑140 (MQ=255)
                                cTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAg                                                                                                  >  1:949101/1‑140 (MQ=255)
                                cTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAg                                                                                                  >  1:695731/1‑140 (MQ=255)
                                    aTTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGatt                                                                                              >  1:1035379/1‑140 (MQ=255)
                                    aTTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCA‑‑‑‑CTTATTTTTTACGCGTAATACAATGAAATAAAAGGatttatt                                                                                          >  2:653268/1‑140 (MQ=255)
                                        attaAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATaaaa                                                                                                   <  1:95956/131‑1 (MQ=255)
                                        attaAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATaaaa                                                                                                   >  2:95956/1‑131 (MQ=255)
                                         ttaAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGatttattt                                                                                         <  1:364995/140‑1 (MQ=255)
                                                tATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCa                                                                                  <  2:1131839/140‑1 (MQ=255)
                                                 atcatcCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCAc                                                                                 >  1:945822/1‑140 (MQ=255)
                                                 atcatcCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCAc                                                                                 >  2:198048/1‑140 (MQ=255)
                                                    atcCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTc                                                                              <  2:811154/140‑1 (MQ=255)
                                                        gTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCcaca                                                                          >  2:355948/1‑140 (MQ=255)
                                                              cTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGa                                                                    >  1:11108/1‑140 (MQ=255)
                                                                  aaaaCATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACcac                                                                <  2:1035379/140‑1 (MQ=255)
                                                                      cATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCg                                                            >  1:805060/1‑140 (MQ=255)
                                                                                    tGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCt                                              >  2:919772/1‑140 (MQ=255)
                                                                                     gCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTc                                             <  1:919772/140‑1 (MQ=255)
                                                                                     gCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTc                                             <  1:1061958/140‑1 (MQ=255)
                                                                                          cTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTg                                        >  1:900251/1‑140 (MQ=255)
                                                                                            gTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAg                                      <  2:562571/140‑1 (MQ=255)
                                                                                             tACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAgg                                     >  1:131308/1‑140 (MQ=255)
                                                                                              aCCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAggg                                    >  1:578765/1‑140 (MQ=255)
                                                                                               ccAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAgggg                                   >  2:101070/1‑140 (MQ=255)
                                                                                                cAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGc                                  <  2:900251/140‑1 (MQ=255)
                                                                                                   gAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAg                                      <  2:300502/133‑1 (MQ=255)
                                                                                                   gAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAg                                      >  1:300502/1‑133 (MQ=255)
                                                                                                   gAACCACCTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTc                               <  2:945822/140‑1 (MQ=255)
                                                                                                     accaccTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAGTAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCt                             >  2:254707/1‑140 (MQ=255)
                                                                                                     accaccTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACt                                         <  1:1042652/128‑1 (MQ=255)
                                                                                                     accaccTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACt                                         >  2:1042652/1‑128 (MQ=255)
                                                                                                     accaccTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCt                             <  1:1129643/140‑1 (MQ=255)
                                                                                                       caccTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGATAAAAAAGCCCCTCGACTGAGGGGCTTCCTGt                           >  2:940044/1‑140 (MQ=255)
                                                                                                       caccTCCTTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGt                           >  2:951533/1‑140 (MQ=255)
                                                                                                         cctcctTAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGttt                         <  1:528984/140‑1 (MQ=255)
                                                                                                              ttAGCCTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAAt                    >  2:1035822/1‑140 (MQ=255)
                                                                                                                  ccTGTGTAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACa                <  2:344158/140‑1 (MQ=255)
                                                                                                                       gtAATCTCCCTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCcac           >  1:522925/1‑140 (MQ=255)
                                                                                                                                cTTGCACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTtgc  <  2:961277/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                        
ATGGAAAATTCTTAGAAACCGATCACATACAGCTGCATTTATTAAGGTTATCATCCGTTTCGCTGAAAAACATAACCCATAAAATGCTAGCTGTACCAGGAACCACCTCCTTAGCCTGTGTAATCTCCCTTACACGGGCTTATTTTTTACGCGTAATACAATGAAATAAAAGGATTTATTTCTGGTCACGTCCACACATTGACCACATCGACAAAAAAGCCCCTCGACTGAGGGGCTTCCTGTTTCTAATTACATCCACATAATTTGC  >  NC_000913/2167053‑2167320

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