Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,352,815 G→A 100% A92V (GCA→GTA)  rpsI ← ribosomal protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,352,8150GA92.6% 75.9 / ‑0.4 27A92V (GCA→GTA) rpsIribosomal protein S9
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (10/15);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTAATGATT  >  CP000730/2352670‑2352959
                                                                                                                                                 |                                                                                                                                                
tttGTATATATCGTCCGGCAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCtt                                                                                                                                               >  1:266122/1‑149 (MQ=255)
tttGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCtt                                                                                                                                               >  1:176150/1‑149 (MQ=255)
tttGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACtt                                                                                                                                               <  2:297277/149‑1 (MQ=255)
  tGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCt                                                                                                                                             <  2:122358/149‑1 (MQ=255)
        tatCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTaataat                                                                                                                                       >  2:181259/1‑149 (MQ=255)
                        aCGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATCATGCACGAGCGATTCCGt                                                                                                                       >  2:2789/1‑149 (MQ=255)
                        aCGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGt                                                                                                                       >  1:81420/1‑149 (MQ=255)
                                 gaATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATag                                                                                                              <  1:282912/149‑1 (MQ=255)
                                 gaATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATag                                                                                                              <  2:158474/149‑1 (MQ=255)
                                                   aGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGa                                                                                            <  1:99310/149‑1 (MQ=255)
                                                             aGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAaccaccacc                                                                                  >  2:282267/1‑149 (MQ=255)
                                                                     ttttttACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGt                                                                          <  1:111279/149‑1 (MQ=255)
                                                                                     aTACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATa                                                          >  1:90100/1‑149 (MQ=255)
                                                                                       aCGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGt                                                        >  2:186546/1‑149 (MQ=255)
                                                                                         gTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTa                                                      <  2:14633/149‑1 (MQ=255)
                                                                                         gTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTa                                                      >  1:226060/1‑149 (MQ=255)
                                                                                         gTGGGTCACGAGTAAGTAATCCAGCGCGTTATAAAGAACCTCTGGATGCAGGATCTACTGCTAATAATGGACGAGCGAGTCCGTGACGGATAGCTTAAGCTAGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTa                                                      >  1:44423/1‑149 (MQ=255)
                                                                                            ggTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCt                                                   <  2:143605/149‑1 (MQ=255)
                                                                                            ggTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCt                                                   <  1:268871/149‑1 (MQ=255)
                                                                                                        gtaaTCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACAACCATGAACGTTAACTAAAACATCATAGTTACCTCTAGTTTCAGtt                                       >  2:103082/1‑149 (MQ=255)
                                                                                                                  gcgTTTTAAAGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTACAACATCATAGTTACCTTTAGTTTCAGTTACATCAAAtg                             >  2:94072/1‑149 (MQ=255)
                                                                                                                         aaaGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATggttggtt                      <  1:282267/149‑1 (MQ=255)
                                                                                                                         aaaGAACCTCTGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATggttggtt                      <  2:282023/149‑1 (MQ=255)
                                                                                                                                ctctGTATTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCt               <  2:29062/149‑1 (MQ=255)
                                                                                                                                      aTTCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAAtt         <  2:269044/149‑1 (MQ=255)
                                                                                                                                        tCAGGATCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTaa       <  2:185515/149‑1 (MQ=255)
                                                                                                                                             aTCTACTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTAATGAtt  <  1:143840/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGTTCTTTCATACGTGGGTCACGAGTAAGTAATCCAGCGCGTTTTAAAGAACCTCTGTATTCAGGATCTGCTTCTAATAATGCACGAGCGATTCCGTGACGGATAGCTTGAGCTTGTCCAGTGAAACCACCACCATGAACGTTAACTAAAACATCATAGTTACCTTTAGTTTCAGTTACATCAAATGGTTGGTTTAAGTCTAAAATTAATGATT  >  CP000730/2352670‑2352959

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