Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,178,296	A→T	intergenic (+8/‑151)	nusG → / → rplK	transcription termination factor/50S ribosomal subunit protein L11

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,178,296	0	A	T	92.1%	243.4 / 10.5	76	intergenic (+8/‑151)	nusG/rplK	transcription termination factor/50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (43/27);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTA  >  NC_000913/4178200‑4178395                                                                                                 |                                                                                                    aCGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCCTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCa                                                                                                  >  2:926210/1‑100 (MQ=255)   gagaAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCaaa                                                                                                <  1:1210264/100‑1 (MQ=255)   gagaAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCaaa                                                                                                >  2:750638/1‑100 (MQ=255)         tctcGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAgcg                                                                                          <  2:513411/100‑1 (MQ=255)          ctcGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAgcgg                                                                                         >  1:675068/1‑100 (MQ=255)            cGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAgcggcg                                                                                       <  2:1171636/100‑1 (MQ=255)             ntCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGa                                                                                      >  1:1639045/2‑100 (MQ=255)               cTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGAtt                                                                                    >  1:572088/1‑100 (MQ=255)                tGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGAttt                                                                                   >  1:808725/1‑100 (MQ=255)                     gtgtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCg                                                                              <  2:1041326/100‑1 (MQ=255)                       gtCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAAACCAGCGTTCAAAAAAGCGGCGATTTAATCGtt                                                                            >  1:596539/1‑100 (MQ=255)                        tCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAAAAAAAAATAAAAAAAAGCGGCGATTTAATCGTTg                                                                           <  1:547250/100‑1 (MQ=255)                         cTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGc                                                                          >  2:955403/1‑100 (MQ=255)                         cTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGc                                                                          >  1:1739377/1‑100 (MQ=255)                         cTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGc                                                                          <  2:1124618/100‑1 (MQ=255)                          tGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGca                                                                         <  2:1530649/100‑1 (MQ=255)                          tGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGca                                                                         >  1:68882/1‑100 (MQ=255)                           gTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGcac                                                                        >  2:567421/1‑100 (MQ=255)                            tttCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGcaca                                                                       >  2:373140/1‑100 (MQ=255)                                 aTCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCg                                                                  >  2:352538/1‑100 (MQ=255)                                 aTCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCg                                                                  >  2:111420/1‑100 (MQ=255)                                   cTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTg                                                                >  1:146463/1‑100 (MQ=255)                                   cTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTg                                                                <  2:295786/100‑1 (MQ=255)                                     tCGGTCGTGCGACCGCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTga                                                               <  2:1912126/99‑1 (MQ=255)                                     tCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTgag                                                              <  1:567419/100‑1 (MQ=255)                                      cGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTgaga                                                             <  1:1384088/100‑1 (MQ=255)                                         tCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTg                                                          <  2:401610/100‑1 (MQ=255)                                           gTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTggn                                                        <  2:419907/100‑2 (MQ=255)                                           gTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGa                                                        >  1:1030291/1‑100 (MQ=255)                                           gTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGa                                                        <  1:953961/100‑1 (MQ=255)                                            tGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGaa                                                       <  1:538010/100‑1 (MQ=255)                                             gCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAAt                                                      >  1:1061507/1‑100 (MQ=255)                                               gcccccGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACCAGGCGTGAGATTGGAAtcc                                                    >  1:825972/3‑98 (MQ=255)                                                 ccccGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACaa                                                  >  2:1320177/1‑100 (MQ=255)                                                  cccGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAAt                                                 <  2:515701/100‑1 (MQ=255)                                                   ccGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAAtt                                                >  2:185757/1‑100 (MQ=255)                                                      gTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTcg                                             <  1:373138/100‑1 (MQ=255)                                                      gTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTcg                                             <  1:1539966/100‑1 (MQ=255)                                                       tAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTcgc                                            <  2:275015/100‑1 (MQ=255)                                                         gagCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCgcgc                                          >  2:278793/1‑100 (MQ=255)                                                          agCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGcc                                         >  2:1352528/1‑100 (MQ=255)                                                             tGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCttt                                      >  2:778482/1‑100 (MQ=255)                                                             tGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCttt                                      >  2:1254072/1‑100 (MQ=255)                                                              ggACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCtttt                                     >  1:1633233/1‑100 (MQ=255)                                                                 cTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAGTTTCGCGCCTTTTGtt                                  >  1:1137228/1‑100 (MQ=255)                                                                 cTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGtt                                  >  2:738793/1‑100 (MQ=255)                                                                  ttCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGttt                                 <  2:517870/100‑1 (MQ=255)                                                                   tCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGtttt                                >  1:1867910/1‑100 (MQ=255)                                                                   tCAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGtttt                                >  1:36643/1‑100 (MQ=255)                                                                    cAGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGttttt                               >  1:259964/1‑100 (MQ=255)                                                                     aGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTa                              <  1:1453155/100‑1 (MQ=255)                                                                     aGCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTa                              <  1:860178/100‑1 (MQ=255)                                                                      gCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                             >  2:1216267/1‑100 (MQ=255)                                                                      gCCAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                             <  2:1591065/100‑1 (MQ=255)                                                                       ccAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATg                            >  1:1075450/1‑100 (MQ=255)                                                                       ccAGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATg                            <  2:23512/100‑1 (MQ=255)                                                                         aGGTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATggg                          <  1:1446879/100‑1 (MQ=255)                                                                          ggTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGc                         >  2:1119921/1‑100 (MQ=255)                                                                          ggTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGc                         >  1:835774/1‑100 (MQ=255)                                                                          ggTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGc                         >  1:1874545/1‑100 (MQ=255)                                                                           gTTGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGcc                        >  2:1491886/1‑100 (MQ=255)                                                                             tGAAAAAGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCtt                      >  1:960008/1‑100 (MQ=255)                                                                               aaaaaGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGc                    >  1:1141774/1‑100 (MQ=255)                                                                                aaaaGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                   >  2:1739883/1‑100 (MQ=255)                                                                                aaaaGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                   <  1:1631162/100‑1 (MQ=255)                                                                                   aGCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGt                >  2:583805/1‑100 (MQ=255)                                                                                    gCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTa               >  1:1543314/1‑100 (MQ=255)                                                                                    gCCTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTa               >  2:1166615/1‑100 (MQ=255)                                                                                     ccTAACCCAGCTTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaa              <  2:377879/100‑1 (MQ=255)                                                                                      cTAACCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaa             <  2:1639055/100‑1 (MQ=255)                                                                                        aaCCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAAc           >  2:631153/1‑100 (MQ=255)                                                                                        aaCCCAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAAc           <  2:986061/100‑1 (MQ=255)                                                                                          cccAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGa         >  2:1189559/1‑100 (MQ=255)                                                                                          cccAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGa         >  2:1191472/1‑100 (MQ=255)                                                                                           ccAGCGTTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAt        >  1:1684428/1‑100 (MQ=255)                                                                                                gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt   >  2:4328/1‑100 (MQ=255)                                                                                                 aTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTta  <  1:1637695/100‑1 (MQ=255)                                                                                                 |                                                                                                    ACGAGAAATCTCGTCTGAAAGTGTCTGTTTCTATCTTCGGTCGTGCGACCCCGGTAGAGCTGGACTTCAGCCAGGTTGAAAAAGCCTAACCCAGCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTA  >  NC_000913/4178200‑4178395

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