Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,327,216	A→G	18.6%	D122G (GAT→GGT)	rluB →	23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,327,216	0	A	G	18.6%	150.4 / 31.4	78	D122G (GAT→GGT)	rluB	23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base G (4/10);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

GTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTG  >  NC_000913/1327122‑1327314                                                                                               |                                                                                                   gtgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                                >  2:2386201/1‑99 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                         <  2:895515/100‑4 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                         <  1:80764/100‑4 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                         <  1:78980/100‑4 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                         <  2:208886/100‑4 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                         <  1:1601321/100‑4 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                         >  2:1080232/1‑100 (MQ=255)       gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                         >  1:1079967/1‑100 (MQ=255)        aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagca                                                                                        <  1:845602/100‑5 (MQ=255)        aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagca                                                                                        <  2:1794337/100‑5 (MQ=255)        aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                        <  1:1239907/100‑1 (MQ=255)        aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTAGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                        <  1:1666696/100‑1 (MQ=255)          cGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaa                                                                                      >  2:363460/1‑100 (MQ=255)          cGTCGGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGGCGCTTGGACGTTAATAGCTGTGGTCTGCTGCTGCTCACCATCGGTGGTGAAcagcagc                                                                                      <  2:102756/100‑7 (MQ=255)           gTCTGCCAAAACTGCGTGGCGCACGCTGGATTTCCGTGTGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCGCCATCGGTGGTGAAcagcag                                                                                       <  2:1536682/98‑6 (MQ=255)           gTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTTGGTCGGTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaaa                                                                                     <  1:720780/100‑1 (MQ=255)            tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagcag                                                                                    <  2:316865/100‑9 (MQ=38)            tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAc                                                                                    <  1:2334034/100‑1 (MQ=255)               gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                                 >  1:110073/1‑100 (MQ=255)               gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                                 >  2:259052/1‑100 (MQ=255)               gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGCGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                                 >  2:741096/1‑100 (MQ=255)                ccAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGtt                                                                                >  1:1733389/1‑100 (MQ=255)                 cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                               <  1:619241/100‑1 (MQ=255)                 cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                               <  2:1591597/100‑1 (MQ=255)                  aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTa                                                                              <  2:1126622/100‑1 (MQ=255)                  aaaaCTGCGTGGCGCACGCTGGATCGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGTGGGTGAACTGGCAAACCGtt                                                                                <  2:1453132/98‑1 (MQ=255)                     aCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                           <  2:1648569/100‑1 (MQ=255)                         cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                       >  1:322748/1‑100 (MQ=255)                         cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                       <  2:1709708/100‑1 (MQ=255)                            ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCaa                                                                    >  1:1187969/1‑100 (MQ=255)                              cgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGc                                                                  >  1:23775/1‑100 (MQ=255)                                  cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                              <  1:1065171/100‑1 (MQ=255)                                  cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                              >  1:2035678/1‑100 (MQ=255)                                   gCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                             >  1:1776696/1‑100 (MQ=255)                                    cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGaa                                                            >  2:2242541/1‑100 (MQ=255)                                    cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGaa                                                            >  2:2332938/1‑100 (MQ=255)                                       gATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGtt                                                         <  1:819782/100‑1 (MQ=255)                                       gATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGtt                                                         <  2:2288676/100‑1 (MQ=255)                                       gATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGtt                                                         >  1:1359902/1‑100 (MQ=255)                                           gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaan                                                     <  1:393249/100‑2 (MQ=255)                                           gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                     <  1:1926435/100‑1 (MQ=255)                                            ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGCACTGGCAATACATTTATTTCACCCAAACAGTGAAGTTGAACg                                                    >  2:1673961/1‑100 (MQ=255)                                            ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                    <  2:2007207/100‑1 (MQ=255)                                            ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                    <  1:1869381/100‑1 (MQ=255)                                             cGTGTGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTGTAATGCACCCAAGCCGTGAAGTTGaa                                                      >  2:2268618/1‑97 (MQ=255)                                                 ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGaa                                                <  2:2141763/99‑1 (MQ=255)                                                  gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCCAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                              >  2:1223679/1‑100 (MQ=255)                                                     gCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAACGCACCCAAGCCGTGAAGTTGAACGTGAATATGc                                           <  2:717004/100‑1 (MQ=255)                                                      cttcgaggatAATACCTCTGGTCTCCTGCTGTTCGCCACCGATGGTGAAGTGGCAATCCGTTTAACGCTCCCAAGCCGGGAAGTTGATCGTGAATATGcc                                          <  2:1745559/91‑1 (MQ=255)                                                      cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGcc                                          >  1:990065/1‑100 (MQ=255)                                                        tGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgt                                        <  2:455419/100‑1 (MQ=255)                                                            cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                     >  1:1660015/1‑99 (MQ=255)                                                             gTTAATAGCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcgtg                                   <  2:1550627/100‑1 (MQ=255)                                                             gTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcgtg                                   <  2:520788/100‑1 (MQ=255)                                                              ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCgtgt                                  >  1:1316699/1‑100 (MQ=255)                                                                 aTACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                               >  1:1006419/1‑100 (MQ=255)                                                                   aCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                             >  2:1823223/1‑100 (MQ=255)                                                                    ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                            <  2:749631/100‑1 (MQ=255)                                                                        tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAATCCGTTTAATGCACCCAAGTCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTc                          <  2:391275/98‑1 (MQ=255)                                                                               ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg                 >  2:2291238/1‑100 (MQ=255)                                                                               ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg                 >  1:2286020/1‑100 (MQ=255)                                                                                 gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac               >  1:1756012/1‑100 (MQ=255)                                                                                 gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac               >  2:1019745/1‑100 (MQ=255)                                                                                 ccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac               >  1:774235/4‑100 (MQ=255)                                                                                  ctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg              >  1:2013507/1‑100 (MQ=255)                                                                                  cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg              >  1:551622/3‑100 (MQ=255)                                                                                   tgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCATGTATTTGGTCAGGTTGACGAcgc             <  2:322746/100‑1 (MQ=255)                                                                                   agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  1:1785604/2‑100 (MQ=255)                                                                                   agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  1:1308258/2‑100 (MQ=255)                                                                                    gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg            <  1:2155643/100‑1 (MQ=255)                                                                                    gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg            >  1:1077073/1‑100 (MQ=255)                                                                                      tCACCACCGATGGTGAACCGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaa          >  2:537256/1‑100 (MQ=255)                                                                                         ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACt       >  2:1617896/1‑100 (MQ=255)                                                                                          caccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg      >  2:1151924/1‑100 (MQ=255)                                                                                          caccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTg      <  1:2332959/100‑1 (MQ=255)                                                                                           accGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc     >  1:386653/1‑100 (MQ=255)                                                                                             cGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt   <  2:1079962/100‑1 (MQ=255)                                                                                              gATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTg  <  1:1556012/100‑1 (MQ=255)                                                                                               |                                                                                                   GTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTG  >  NC_000913/1327122‑1327314

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