Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,636,499	G→A	N286N (AAC→AAT)	der ←	50S ribosomal subunit stability factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,636,499	0	G	A	97.9%	160.9 / NA	48	N286N (AAC→AAT)	der	50S ribosomal subunit stability factor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/22);  total (26/22)
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

CCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACGTTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCATTGGG  >  NC_000913/2636372‑2636628                                                                                                                                |                                                                                                                                  ccTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGc                                                                                                                              >  2:425873/1‑133 (MQ=255)     aCTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATAGCTTCGCGCGCATAATTGACTAACATCACCATGTTGGCGTCt                                                                                                                          >  2:278113/1‑133 (MQ=255)           gCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATg                                                                                                                    >  2:31978/1‑133 (MQ=255)                cccACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTg                                                                                                               >  1:318749/1‑133 (MQ=255)                     ttATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACg                                                                                                          >  1:445001/1‑133 (MQ=255)                          caccacAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTg                                                                                                     >  1:635904/1‑133 (MQ=255)                          caccacAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTg                                                                                                     >  2:147060/1‑133 (MQ=255)                           accacAATAACAAGTCAGCGCCGACTATTGTGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGa                                                                                                    <  1:196748/133‑1 (MQ=255)                            ccacAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGAt                                                                                                   >  1:448174/1‑133 (MQ=255)                             cacaATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTAGCGTCTTGAATGGCATGCAACATTTTGAtt                                                                                                  >  1:185297/1‑133 (MQ=255)                              acaATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTa                                                                                                 >  2:334846/1‑133 (MQ=255)                                  tGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGga                                                                                             >  1:566642/1‑133 (MQ=255)                                  tGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGga                                                                                             >  1:64393/1‑133 (MQ=255)                                      aaGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAAt                                                                                         >  1:372557/1‑133 (MQ=255)                                        gTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAAttt                                                                                       <  2:30921/133‑1 (MQ=255)                                          gAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTtct                                                                                     >  1:7398/1‑133 (MQ=255)                                               cccACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACa                                                                                <  1:402437/133‑1 (MQ=255)                                                  aCTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCa                                                                             >  2:130406/1‑133 (MQ=255)                                                            aTAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGCGAttt                                                                   >  1:89471/1‑133 (MQ=255)                                                              aaaGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTg                                                                 >  1:89946/1‑133 (MQ=255)                                                              aaaGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTg                                                                 >  2:45493/1‑133 (MQ=255)                                                              aaaGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTg                                                                 >  2:589673/1‑133 (MQ=255)                                                                aGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGcc                                                               <  2:635904/133‑1 (MQ=255)                                                                aGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAACCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTCTTGATTACTGAGAATTTCTCTACAGAATCGGTGATTTTgac                                                               >  2:286785/1‑131 (MQ=255)                                                                               gaTCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGcc                                                <  2:89946/133‑1 (MQ=255)                                                                                 tCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAg                                              <  1:600947/133‑1 (MQ=255)                                                                                   cTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCCGTGATTTGGCCCCGTTTACGTACGCCAGCg                                            >  1:275971/1‑133 (MQ=255)                                                                                         ggAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCa                                      <  1:589673/133‑1 (MQ=255)                                                                                          gAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCaa                                     <  2:346047/133‑1 (MQ=255)                                                                                          gAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCaa                                     <  2:543972/133‑1 (MQ=255)                                                                                          gAAATACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCaa                                     <  2:448174/133‑1 (MQ=255)                                                                                           aaaTACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAAt                                    <  1:147060/133‑1 (MQ=255)                                                                                           aaaTACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAAt                                    <  1:45493/133‑1 (MQ=255)                                                                                           aaaTACCTTCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAAt                                    <  2:445001/133‑1 (MQ=255)                                                                                            aaTACCTTCGTGCGCTTCAATCACAAACGTCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATg                                   <  2:228998/133‑1 (MQ=255)                                                                                                 cTTCGCGCGCATCAATCACTAATATCACCATATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGcac                              <  1:286785/133‑1 (MQ=255)                                                                                                  ttCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGcaca                             <  1:334846/133‑1 (MQ=255)                                                                                                   tCGCGCGCATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACAt                            >  1:21766/1‑133 (MQ=255)                                                                                                     gcgcgcATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCTGTATCAATGAGCACATAc                          >  1:228998/1‑133 (MQ=255)                                                                                                     gcgcacATCAATCACTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTGTGCCGCGTCTACGTACACCAGCGGTGTCAAGGAGCACATAc                          >  2:168933/1‑133 (MQ=255)                                                                                                            tcaatcaCTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTc                   <  1:130406/133‑1 (MQ=255)                                                                                                              aatcaCTAACCTCACCACATTGGCTTCTCCAATGGCCTGCATCGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCa                 <  2:185297/133‑1 (MQ=255)                                                                                                              aatcaCTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCa                 <  2:566642/133‑1 (MQ=255)                                                                                                              aatcaCTAACATCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCa                 <  2:64393/133‑1 (MQ=255)                                                                                                                         tCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCAt      <  2:372557/133‑1 (MQ=255)                                                                                                                         tCACCACATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCAt      >  1:356098/1‑133 (MQ=255)                                                                                                                          caccacATTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCAtt     >  1:142138/1‑133 (MQ=255)                                                                                                                           accacATTGGCGTCTACAATGGCCGGCAACGTATTGATTACTGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCATTg    <  2:248162/133‑1 (MQ=255)                                                                                                                             cacaTTGGCGTCTTCAATCGCCTGCAACGTGTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCATTggg  <  1:278113/133‑1 (MQ=255)                                                                                                                                |                                                                                                                                  CCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAATAAAGCCCAGCAGCGAGAGATCCTGATCGGAAATACCTTCGCGCGCATCAATCACTAACATCACCACGTTGGCGTCTTCAATGGCCTGCAACGTTTTGATTACGGAGAATTTCTCTACAGCATCGGTGATTTTGCCGCGTTTACGTACGCCAGCGGTGTCAATGAGCACATACTCACGTCCATCGCGTTCCATTGGG  >  NC_000913/2636372‑2636628

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