Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,289,658 T→C S138S (AGT→AGC)  rssA → putative patatin‑like phospholipase RssA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,289,6580TC96.9% 106.8 / NA 32S138S (AGT→AGC) rssAputative patatin‑like phospholipase RssA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/16);  minor base G (0/1);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGTTGACCTGCAGCACGA  >  NC_000913/1289525‑1289792
                                                                                                                                     |                                                                                                                                      
gcgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCa                                                                                                                                  >  1:104015/1‑140 (MQ=255)
 cgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAg                                                                                                                                 >  2:402681/1‑140 (MQ=255)
 cgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAg                                                                                                                                 >  1:595737/1‑140 (MQ=255)
   aaaTAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGa                                                                                                                               <  1:402681/140‑1 (MQ=255)
           ccGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGc                                                                                                                       <  2:143255/140‑1 (MQ=255)
           ccGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGc                                                                                                                       <  2:104015/140‑1 (MQ=255)
                              aTTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGc                                                                                                    <  1:455/140‑1 (MQ=255)
                                     ccGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActggc                                                                                             <  2:451739/140‑1 (MQ=255)
                                                     gTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActgg                                                                                              >  1:77005/1‑123 (MQ=255)
                                                     gTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActgg                                                                                              <  2:77005/123‑1 (MQ=255)
                                                                  tAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAAtt                                                                <  1:400190/140‑1 (MQ=255)
                                                                        cGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTAtt                                                          >  2:149320/1‑140 (MQ=255)
                                                                          ggACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAGCGGCTACTCGCTTGATGGTGGAGCATTCGTTAACCCAATTCCTATTTc                                                        >  2:339337/1‑140 (MQ=255)
                                                                           gACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTcc                                                       >  2:509849/1‑140 (MQ=255)
                                                                             cGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCt                                                     >  2:667037/1‑140 (MQ=255)
                                                                              gTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTc                                                    >  2:486894/1‑140 (MQ=255)
                                                                               tGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCa                                                   <  1:617064/140‑1 (MQ=255)
                                                                                        gTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCAt                                          <  1:693508/140‑1 (MQ=255)
                                                                                           tACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAATCCAATTCCTATTTCCCTCACGCGTGCATTgg                                       >  2:507338/1‑140 (MQ=255)
                                                                                            aCTGAAGGCGATCTCCATCTTGCTATTCGTGCATCATCCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTggg                                      <  1:288603/140‑1 (MQ=255)
                                                                                            aCTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTggg                                      <  1:486894/140‑1 (MQ=255)
                                                                                              tGAAGGCGATCTCCATCTTGCTCTTCGCGCAACATGCAGGATTCCAGGTCTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTggggg                                    <  2:270911/140‑1 (MQ=255)
                                                                                              tGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTggggg                                    <  1:667037/140‑1 (MQ=255)
                                                                                                aaGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTggg                                      >  1:155269/1‑136 (MQ=255)
                                                                                                aaGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTggg                                      <  2:155269/136‑1 (MQ=255)
                                                                                                             tCTTGCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAg                     >  1:618014/1‑140 (MQ=255)
                                                                                                                 gCTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGt                 >  2:381071/1‑140 (MQ=255)
                                                                                                                  cTATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTCCATCGGCGGCTGAGATTGTGATAGCGGtt                >  1:304959/1‑140 (MQ=255)
                                                                                                                   tATTCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGCGGCTGATATTGTGATAGCGGTTg               <  2:175342/140‑1 (MQ=255)
                                                                                                                     ttCGCGCATCATGCAGCATTCCAGGACTCATGGCACCTGTTGTACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGTTGAc             >  2:10760/1‑140 (MQ=255)
                                                                                                                             tcatGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGTTGACCTgcagca     <  1:31132/140‑1 (MQ=255)
                                                                                                                              catGCAGCATTCCAGTACTCATCGCTCCTTTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGTTGACCTGCAGCAc    <  1:339337/140‑1 (MQ=255)
                                                                                                                                tGCAGCATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGTTGACCTGCAGCACGa  <  1:509849/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
GCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTCGTTAACCCAATTCCTATTTCCCTCACGCGTGCATTGGGGGCTGATATTGTGATAGCGGTTGACCTGCAGCACGA  >  NC_000913/1289525‑1289792

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