Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,966,830	G→A	P138P (CCG→CCA)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,830	0	G	A	98.0%	186.8 / NA	51	P138P (CCG→CCA)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/26);  minor base T (0/1);  total (24/27)
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

AACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAC  >  NC_000913/3966684‑3966970                                                                                                                                                   |                                                                                                                                             aaCCCCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACt                                                                                                                                            <  1:888847/149‑1 (MQ=255)   ccTCCCCACTGGTAATCCCTTCTCGGGTAAGTTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCAc                                                                                                                                        <  1:844292/151‑1 (MQ=255)         cgcTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCaaa                                                                                                                                   <  2:863898/148‑1 (MQ=255)             ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAActc                                                                                                                                <  2:135987/149‑1 (MQ=255)                gATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGt                                                                                                                            <  1:148318/150‑1 (MQ=255)                 aTACCATCTCTGGAAATATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGt                                                                                                                            <  2:465497/149‑1 (MQ=255)                    ccATCTCTGGAAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCg                                                                                                                       <  2:273809/151‑1 (MQ=255)                         tctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATgg                                                                                                                  >  1:76596/1‑151 (MQ=255)                              tAAAATCCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATgg                                                                                                                  <  1:96864/146‑1 (MQ=255)                                   tcggcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTa                                                                                                         <  2:523532/147‑1 (MQ=255)                                      gcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACg                                                                                                      >  2:97267/1‑150 (MQ=255)                                        ccgccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                                                                                                   <  2:39740/151‑1 (MQ=255)                                            cgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAc                                                                                               >  1:158832/1‑151 (MQ=255)                                              aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTg                                                                                             >  2:786250/1‑151 (MQ=255)                                                  aagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                                                                                            >  1:577139/1‑148 (MQ=255)                                                   agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                                                                                            >  1:701586/1‑147 (MQ=255)                                                   agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                                                                                            >  2:794390/1‑147 (MQ=255)                                                   agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt                                                                                        <  1:349278/151‑1 (MQ=255)                                                       gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTaa                                                                                    <  1:196308/151‑1 (MQ=255)                                                       gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTaa                                                                                    <  2:611141/151‑1 (MQ=255)                                                        aaCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                   >  1:72173/1‑151 (MQ=255)                                                             tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCt                                                                               >  2:489249/1‑150 (MQ=255)                                                             tATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTc                                                                              >  1:257584/1‑151 (MQ=255)                                                              aTTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTcg                                                                             <  1:728929/151‑1 (MQ=255)                                                                     gctgctGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACt                                                                       >  2:203245/1‑150 (MQ=255)                                                                      ctgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTgg                                                                     >  2:175122/1‑151 (MQ=255)                                                                      ctgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAAAAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTg                                                                      >  2:871263/1‑150 (MQ=255)                                                                      ctgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCAGCAACAAAATCGTCTTTGAGACCTTAGCCCGACTGCACTCAAACTCTCGTCTGCGTATGCAACGTGTTAGCGGTTCTAATGAAGATTTAACTGCTCGCGTACTg                                                                      <  2:139221/150‑1 (MQ=255)                                                                               gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCAt                                                            <  1:791452/151‑1 (MQ=255)                                                                                      aaGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTAt                                                     >  2:789570/1‑151 (MQ=255)                                                                                        gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCg                                                   >  1:236890/1‑151 (MQ=255)                                                                                        gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCg                                                   <  1:221075/151‑1 (MQ=255)                                                                                        gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCg                                                   <  1:389285/151‑1 (MQ=255)                                                                                          tAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                    >  1:517715/1‑148 (MQ=255)                                                                                               tCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTg                                             <  1:723892/150‑1 (MQ=255)                                                                                                  aCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTc                                          >  2:640817/1‑150 (MQ=255)                                                                                                  aCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCa                                         <  2:694030/151‑1 (MQ=255)                                                                                                     aaCCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCg                                      >  2:203030/1‑151 (MQ=255)                                                                                                       ccTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTg                                    >  2:707296/1‑151 (MQ=255)                                                                                                         tGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGt                                  >  1:696988/1‑151 (MQ=255)                                                                                                          gAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTgg                                   <  1:714364/149‑1 (MQ=255)                                                                                                           aaaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTc                                 >  2:420494/1‑150 (MQ=255)                                                                                                            aaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTg                               <  2:701537/151‑1 (MQ=255)                                                                                                                gCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtg                           >  2:237655/1‑151 (MQ=255)                                                                                                                 cccGCAACAAAATCCTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                          >  1:276132/1‑151 (MQ=255)                                                                                                                              ccTCTTTGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaaa             <  1:897876/151‑1 (MQ=255)                                                                                                                                  tttGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg         >  2:398599/1‑151 (MQ=255)                                                                                                                                    tGAGAACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg        <  1:93877/150‑1 (MQ=255)                                                                                                                                      agaTCTTAACCCCACTGCACGCAACCTCTCGTTGGCGTATGAAACGTGGTAACGGTGTTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCTCCTAGCGGTCGTGGTCAGCGTGGTCTAATTGTGGCACCGCCGAAAGCCGGTaa     <  1:275284/151‑1 (MQ=255)                                                                                                                                      agaACTTAACCCCACTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaa     <  1:616410/151‑1 (MQ=255)                                                                                                                                          cTTACCCCCTCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAc  <  2:108881/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAC  >  NC_000913/3966684‑3966970

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