Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,454,612 2 bp→AC coding (33‑34/438 nt) unknown → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,454,6120CA87.0% 65.6 / ‑4.5 23I11I (ATC→ATA) unknownlipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (15/5);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 3.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*11,454,6130TC100.0% 72.5 / NA 23S12P (TCT→CCT) unknownlipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/8);  total (15/8)

TAAGCCAATACTGTGTAATAAAGCGGTTTACTATAACGACCCTTATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTCTGCTTGTTTTTTAAGCGCGT‑GCAACTCGTCACCGCCGGAACCGCCACAAATTAATGAAGACGCGCCGGCAATCTATTTAAATAATCAGGCGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCA  >  1/1454369‑1454845
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                         
caaaccgattatatgtAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAg                                                                                                                                                                                                                           <  1:41541/236‑1 (MQ=255)
            atgtAATAAGGCAGTGTATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGt                                                                                                                                                                                                            >  1:412919/2‑251 (MQ=255)
               taataaggcagtgtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCa                                                                                                                                                                                                          >  1:165664/20‑250 (MQ=255)
               taataaggcagtgtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCa                                                                                                                                                                                                          >  1:565257/20‑250 (MQ=255)
                aataaggcagtgtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCAcc                                                                                                                                                                                                        <  1:646724/233‑1 (MQ=255)
                           gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAg                                                                                                                                                                                                                           >  1:466134/8‑222 (MQ=255)
                           gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCAc                                                                                                                                                                                                         <  1:504059/232‑1 (MQ=255)
                           gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAAc                                                                                                                                                                                                >  1:643527/8‑248 (MQ=255)
                           gtattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCg                                                                                                                                                                                              >  1:659897/8‑250 (MQ=255)
                            tattacaaCGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGc                                                                                                                                                                                             <  1:817973/244‑1 (MQ=255)
                                       tccTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAAttgata                                                                                                                                                                                 >  1:274096/2‑247 (MQ=255)
                                            tttttataGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTGTTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGAc                                                                                                                                                                            <  1:497178/250‑1 (MQ=255)
                                              tttataGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAg                                                                                                                                                                                                                           <  1:233148/203‑1 (MQ=255)
                                                                                                                                                             ctgtcagtccaataaattaaatcaacaAAAACTGTGAAAAACCAGGAGAAAAATAAAAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTCCTTTTAAGGCG‑GTGGCAGCTCGACACCACCAAAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATc                                                                        <  1:408089/226‑1 (MQ=255)
                                                                                                                                                                      caatcaatttaatCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                  >  1:393289/12‑247 (MQ=255)
                                                                                                                                                                                           aCTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTAtgt                                         >  1:58895/1‑251 (MQ=255)
                                                                                                                                                                                            cTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTAtgtg                                        >  1:398930/1‑251 (MQ=255)
                                                                                                                                                                                                       cAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGa                             >  1:22401/1‑251 (MQ=255)
                                                                                                                                                                                                            ggaTAATATAAAAAAAATTTTACTTATAAATATGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCCCCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATg                                                                                          <  1:461020/185‑1 (MQ=255)
                                                                                                                                                                                                                               ttACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTg      >  1:94806/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCaa                                                                      >  1:224595/1‑183 (MQ=255)
                                                                                                                                                                                                                                  cTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGca  >  1:892947/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACtt        >  1:499998/1‑241 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                         
TAAGCCAATACTGTGTAATAAAGCGGTTTACTATAACGACCCTTATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTCTGCTTGTTTTTTAAGCGCGT‑GCAACTCGTCACCGCCGGAACCGCCACAAATTAATGAAGACGCGCCGGCAATCTATTTAAATAATCAGGCGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCA  >  1/1454369‑1454845

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