Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,134,891	A→T	8.2%	I494K (ATA→AAA)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,134,891	0	A	T	8.2%	146.8 / 4.7	49	I494K (ATA→AAA)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base T (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTG  >  NC_000913/3134747‑3135035                                                                                                                                                 |                                                                                                                                                 aTCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATcc                                                                                                                                               >  5:241568/1‑148 (MQ=255)  tCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATc                                                                                                                                                >  1:364520/1‑146 (MQ=255)   cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAg                                                                                                                                             >  6:205194/1‑148 (MQ=255)    aGCCTCTGGAACAGGTGTTGGTTGATGCAGAATAAGCTATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCGCGATCGCCAGGACCGGTCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGt                                                                                                                                            <  3:328214/148‑1 (MQ=255)    aGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGt                                                                                                                                            <  7:355718/148‑1 (MQ=255)       cACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTa                                                                                                                                           >  4:280914/1‑146 (MQ=255)                aGGTGTTGGTAGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGTAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCccac                                                                                                                                   >  3:128143/1‑145 (MQ=255)                   tgtTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCaga                                                                                                                               >  8:115957/1‑146 (MQ=255)                   tgtTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCaga                                                                                                                               >  2:8082/1‑146 (MQ=255)                   tgtTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCACGATCGCCAGTACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCaga                                                                                                                               >  2:80968/1‑146 (MQ=255)                   tgtTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCACGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCACCTGCAATGCTATCCAGTACAGCCCACCaga                                                                                                                               >  6:145171/1‑146 (MQ=255)                        gttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAgaaagaaa                                                                                                                         <  2:80030/147‑1 (MQ=255)                        gttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGaaaa                                                                                                                        <  2:106929/148‑1 (MQ=255)                         ttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGaaaaa                                                                                                                       <  7:115957/148‑1 (MQ=255)                                aGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgccgc                                                                                                                <  2:364520/148‑1 (MQ=255)                                          aTTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagta                                                                                                      <  6:106887/148‑1 (MQ=255)                                          aTTACGGATTTTATGACCGCGCGGGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagta                                                                                                      <  2:287636/148‑1 (MQ=255)                                           ttACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagtaa                                                                                                     <  6:241568/148‑1 (MQ=255)                                                 aTTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATAccc                                                                                               <  4:128143/148‑1 (MQ=255)                                                            acgcgTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAgg                                                                                    <  5:145171/147‑1 (MQ=255)                                                                  ggTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagta                                                                                                      <  3:132398/124‑1 (MQ=255)                                                                  ggTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagta                                                                                                      >  4:132398/1‑124 (MQ=255)                                                                   gTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTg                                                                              <  3:280914/147‑1 (MQ=255)                                                                    ttATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGt                                                                            >  6:209323/1‑148 (MQ=255)                                                                     tATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTTTCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTg                                                                              <  1:201544/145‑1 (MQ=255)                                                                     tATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTTTCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTg                                                                              >  2:207757/1‑145 (MQ=255)                                                                     tATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTTTCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTg                                                                              >  2:201544/1‑145 (MQ=255)                                                                     tATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTTTCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTg                                                                              <  1:207757/145‑1 (MQ=255)                                                                     tATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATAccc                                                                                               <  2:353880/128‑1 (MQ=255)                                                                     tATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATAccc                                                                                               >  1:353880/1‑128 (MQ=255)                                                                            aGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgccgcc                                                                                                               <  7:288743/105‑1 (MQ=255)                                                                            aGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgccgcc                                                                                                               >  8:288743/1‑105 (MQ=255)                                                                            aGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagtaa                                                                                                     >  5:10494/1‑115 (MQ=255)                                                                            aGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagtaa                                                                                                     <  6:10494/115‑1 (MQ=255)                                                                                    aTGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGt                                                            <  5:205194/148‑1 (MQ=255)                                                                                            gCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCAcc                                                    >  7:133178/1‑148 (MQ=255)                                                                                            gCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCAcc                                                    >  7:66003/1‑148 (MQ=255)                                                                                            gCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCAcc                                                    >  4:182232/1‑148 (MQ=255)                                                                                                    gCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTcc                                             <  8:66003/147‑1 (MQ=255)                                                                                                      cAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTCAGCTGTAACCCACCGCCGTccac                                           >  5:47890/1‑147 (MQ=255)                                                                                                         gACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccaccacc                                       <  1:80968/148‑1 (MQ=255)                                                                                                         gACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTccaccacc                                       <  1:8082/148‑1 (MQ=255)                                                                                                          aCCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTAc                                                                 >  2:39564/1‑121 (MQ=255)                                                                                                          aCCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTAc                                                                 <  1:39564/121‑1 (MQ=255)                                                                                                               gCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTc                                 <  5:209323/148‑1 (MQ=255)                                                                                                                        ggAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCt                                                                               <  5:96441/93‑1 (MQ=255)                                                                                                                        ggAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCt                                                                               >  6:96441/1‑93 (MQ=255)                                                                                                                                    cAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGagag             >  7:21009/1‑147 (MQ=255)                                                                                                                                              gCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGagttg  >  3:320651/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                 ATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTG  >  NC_000913/3134747‑3135035

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