Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	34,120	(T)9→8	intergenic (+82/‑180)	carB → / → caiF	carbamoyl phosphate synthetase subunit beta/DNA‑binding transcriptional activator CaiF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	34,112	0	T	.	97.7%	170.7 / ‑2.0	44	intergenic (+74/‑188)	carB/caiF	carbamoyl phosphate synthetase subunit beta/DNA‑binding transcriptional activator CaiF
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (17/25);  total (17/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 = 7.53e-01

CCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCC  >  NC_000913/33983‑34238                                                                                                                                  |                                                                                                                               ccGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttt                                                                                                                             >  2:123891/1‑133 (MQ=255)  cGATGCGACTGAAAACGTGATTGCGGTGCACGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttt                                                                                                                            <  1:201196/133‑1 (MQ=255)    aTGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttt                                                                                                                          <  1:403302/133‑1 (MQ=255)    aTGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttt                                                                                                                          <  2:463306/133‑1 (MQ=255)     tGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtttttttt                                                                                                                         <  2:84891/133‑1 (MQ=255)          cTGAAAAAGTTATTTAGGTGCAGGAAATTCACGCACAGATCGAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGA‑CGGTTTTTTTTGCCCa                                                                                                                   <  2:201065/133‑1 (MQ=255)          cTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCa                                                                                                                   <  1:487491/133‑1 (MQ=255)          cTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCa                                                                                                                   <  2:298126/133‑1 (MQ=255)          cTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCa                                                                                                                   <  1:4748/133‑1 (MQ=255)                          ggTGCAGGAAATGCACGCACAGAGCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTTTGATGAGGTCCG‑TTTTTTTTGCCCAATAAGTAAATCTTTTg                                                                                                   >  2:220980/1‑133 (MQ=255)                           gTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTgg                                                                                                  >  1:288136/1‑133 (MQ=255)                            tGCAGGAATTGCACGTACGGATAAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTggg                                                                                                 <  1:235050/133‑1 (MQ=255)                                 gAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                                                            <  1:284327/133‑1 (MQ=255)                                               gATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatga                                                                              >  1:449627/1‑133 (MQ=255)                                                    aataatAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATaa                                                                         >  1:12090/1‑133 (MQ=255)                                                      taataGCGAGTCATGACACATACTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGc                                                                       <  1:144465/133‑1 (MQ=255)                                                          aGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTGAATCTTTTGGGGAATCGATATTTTTAATGAGATAAGCAGGa                                                                   >  2:144465/1‑133 (MQ=255)                                                          aGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTCAATCTTTTCGGGAATCGATATTTTTAATGACATAAGCAGGa                                                                   >  2:205752/1‑133 (MQ=255)                                                          aGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGTATCGATATTTTTGATGACATAAGCAGGa                                                                   >  1:422599/1‑133 (MQ=255)                                                          aGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGa                                                                   >  1:273002/1‑133 (MQ=255)                                                           gCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAt                                                                  <  2:112684/133‑1 (MQ=255)                                                            cGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAtt                                                                 >  1:390096/1‑133 (MQ=255)                                                             gtgtCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAttt                                                                >  2:23159/1‑133 (MQ=255)                                                              tgtCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTa                                                               >  1:298488/1‑133 (MQ=255)                                                              tgtCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTa                                                               >  2:406023/1‑133 (MQ=255)                                                                tCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGc                                                             >  1:86677/1‑133 (MQ=255)                                                                 cATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCt                                                            <  2:288136/133‑1 (MQ=255)                                                                         atatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTa                                                    >  1:127435/1‑133 (MQ=255)                                                                             ttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGa                                                <  1:23159/133‑1 (MQ=255)                                                                                  aTCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTg                                           <  2:298488/133‑1 (MQ=255)                                                                                  aTCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTg                                           <  1:393088/133‑1 (MQ=255)                                                                                          aTTTGATCTAATAACTAACACGCTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  1:220980/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATGAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  2:406899/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  1:353651/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  2:453529/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  2:57749/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  2:19442/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  1:406023/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  2:127435/133‑1 (MQ=255)                                                                                          aTTTGATCGAATAACCAATACGGTGCTCTGATAAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  2:186815/133‑1 (MQ=255)                                                                                           tttGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCa                                   <  1:107323/133‑1 (MQ=255)                                                                                               aTCGAATATCTAATATGGTTCTTTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTa                              <  1:205752/133‑1 (MQ=255)                                                                                                    aTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGat                         <  1:111869/133‑1 (MQ=255)                                                                                                    aTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTGGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGat                         <  1:330095/133‑1 (MQ=255)                                                                                                           aTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCca                  <  2:449627/133‑1 (MQ=255)                                                                                                           aTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCca                  >  2:478751/1‑133 (MQ=255)                                                                                                                 ttCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAAttt            >  2:480986/1‑133 (MQ=255)                                                                                                                 ttCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGAGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAATTTGGATACTATCAATATTTCCACAAttt            >  2:200637/1‑133 (MQ=255)                                                                                                                      tgatgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAata       >  1:426335/1‑133 (MQ=255)                                                                                                                           aGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGcc  <  2:398179/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                               CCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCC  >  NC_000913/33983‑34238

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