Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,543,754	T→C	84.6%	T103A (ACG→GCG)	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,543,754	0	T	C	84.6%	196.0 / 20.9	83	T103A (ACG→GCG)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base C (30/39);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTT  >  NC_000913/1543658‑1543849                                                                                                 |                                                                                                cgataacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCa                                                                                              >  2:855970/1‑100 (MQ=255)   ataacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATc                                                                                            <  2:1320137/100‑1 (MQ=255)     aacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCgg                                                                                          <  2:542986/100‑1 (MQ=255)       cgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCg                                                                                        >  1:104160/1‑100 (MQ=255)        gataaATATCCCAAAAGGAGTATTTGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCg                                                                                        <  1:1475728/99‑1 (MQ=255)            aaTATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCAc                                                                                   >  1:1433102/1‑100 (MQ=255)            aaTATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCAc                                                                                   <  2:1695723/100‑1 (MQ=255)            aaTATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCAc                                                                                   >  2:1914568/1‑100 (MQ=255)             atatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCAcc                                                                                  <  1:2478840/100‑1 (MQ=255)             atatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCAcc                                                                                  <  2:1836011/100‑1 (MQ=255)             atatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCAcc                                                                                  <  1:1995201/100‑1 (MQ=255)                tCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGaa                                                                               >  2:2353129/1‑100 (MQ=255)                tCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGaa                                                                               >  2:336552/1‑100 (MQ=255)                 cccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAat                                                                              >  1:1095276/1‑100 (MQ=255)                 cccAAAAGAAGTATTCGGATTTTGCACGGCAATTACGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAat                                                                              <  1:1647173/100‑1 (MQ=255)                  ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATTAGGATTGCAGTACTAAAAACCGCCCATCAGCGTCCACCGAata                                                                             >  1:2201273/1‑100 (MQ=255)                  ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAata                                                                             >  2:1796979/1‑100 (MQ=255)                  ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAata                                                                             >  2:55115/1‑100 (MQ=255)                  ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAata                                                                             <  2:435361/100‑1 (MQ=255)                   cAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAatat                                                                            >  1:434102/1‑100 (MQ=255)                    aaaaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAAtata                                                                           <  2:1868265/100‑1 (MQ=255)                     aaaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAg                                                                          <  1:1917030/100‑1 (MQ=255)                     aaaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAg                                                                          <  1:2501284/100‑1 (MQ=255)                     aaaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAg                                                                          <  2:45565/100‑1 (MQ=255)                       aGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGc                                                                        <  1:416340/100‑1 (MQ=255)                        ggAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAgg                                                                         >  1:1092412/1‑98 (MQ=255)                        ggAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCa                                                                       <  1:776315/100‑1 (MQ=255)                          aGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCAcc                                                                     <  1:423804/100‑1 (MQ=255)                          aGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCAcc                                                                     >  2:569266/1‑100 (MQ=255)                          aGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCAcc                                                                     <  2:1443109/100‑1 (MQ=255)                           gTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCa                                                                    <  2:706787/100‑1 (MQ=255)                            tATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCAt                                                                   >  2:236403/1‑100 (MQ=255)                             aTTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATa                                                                  >  1:2400851/1‑100 (MQ=255)                              ttCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATaa                                                                 >  1:1612292/1‑100 (MQ=255)                                cGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAg                                                               >  2:388217/1‑100 (MQ=255)                                cGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAg                                                               >  1:977489/1‑100 (MQ=255)                                cGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAg                                                               <  2:970533/100‑1 (MQ=255)                                   aTTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAg                                                            <  2:785101/100‑1 (MQ=255)                                      ttGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAg                                                         <  2:498209/100‑1 (MQ=255)                                       tGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAgg                                                        <  1:1245961/100‑1 (MQ=255)                                        gCACGGCAATTCCGAGCCAGACGCAAGGAATAATCATGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAggg                                                       <  1:1426616/100‑1 (MQ=255)                                        gCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATAGGCGTCCACCGAATATAGGCACCATAAAGGAGTAggg                                                       <  1:137015/100‑1 (MQ=255)                                         cACGGCAATCCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGa                                                      >  2:926028/1‑100 (MQ=255)                                           cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAac                                                    <  2:1612301/100‑1 (MQ=255)                                             gCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATCTAGGCACCATAAAGGAGTAGGGAacac                                                  >  2:994644/1‑100 (MQ=255)                                             gCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAacac                                                  >  2:2281679/1‑100 (MQ=255)                                             gCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCAGCGTCCACCGAATATAAGCACAATAAAGGAGTACGGAac                                                    >  1:1881446/1‑98 (MQ=38)                                             gCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAacac                                                  >  2:2417804/1‑100 (MQ=255)                                                  tCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaat                                             <  2:886731/100‑1 (MQ=255)                                                   ccGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGACCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaata                                            >  1:40431/1‑100 (MQ=255)                                                   ccGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACAAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaata                                            >  1:929141/1‑100 (MQ=255)                                                        ccAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATgcg                                       <  1:682711/100‑1 (MQ=255)                                                        ccAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATgcg                                       <  1:797707/100‑1 (MQ=255)                                                              gCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGaa                                 <  1:536421/100‑1 (MQ=255)                                                               cAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGaaa                                >  2:2239385/1‑100 (MQ=255)                                                                   gAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGa                            <  2:903449/100‑1 (MQ=255)                                                                        atCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTa                       >  2:672930/1‑100 (MQ=255)                                                                          cAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAAt                     >  1:1238260/1‑100 (MQ=255)                                                                           aGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATg                    >  2:691314/1‑100 (MQ=255)                                                                           aGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATg                    <  2:929144/100‑1 (MQ=255)                                                                            ggATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGc                   <  1:253303/100‑1 (MQ=255)                                                                            ggATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGc                   >  2:1239122/1‑100 (MQ=255)                                                                                tGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGtt               <  2:940818/100‑1 (MQ=255)                                                                                tGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCAGAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGtt               <  1:926025/100‑1 (MQ=255)                                                                                 gCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTta              >  1:2433142/1‑100 (MQ=255)                                                                                    gTACTAAAAACCGCCCATCGGCGTCCACCGAATAGAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaat            <  1:2436733/99‑1 (MQ=255)                                                                                    gTAATAAAAACCGCCCATCGGCGTCCACCGAAAATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGTAACGGAGGGTAATGCGGTtaat            <  1:1501948/99‑1 (MQ=255)                                                                                     tACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaataa          <  1:2031584/100‑1 (MQ=255)                                                                                     tACTAAAAACCGCCCAGCGGCGTCCACCGTATAGAGGCACGATAAAGGAGTAGGGAACATGCAATAAGGCGCCAGAAACGGAGGGTAATGCGGTtaataa          <  1:891373/100‑1 (MQ=255)                                                                                      cgatctaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaa         <  1:1916505/94‑1 (MQ=255)                                                                                      aCTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaa         <  1:921917/100‑1 (MQ=255)                                                                                       cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa        <  1:1007186/100‑1 (MQ=255)                                                                                       cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa        >  2:1521059/1‑100 (MQ=255)                                                                                         aaaatCCGTCGATCGGCGTCCACCGAATATAAGCAACACAAACTAGCGGGTACCACGCACTAATGCGCATGTAAGGGAGCGTAATGCGCCTAATGAagag      >  1:2506598/1‑100 (MQ=255)                                                                                         aaaaaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag      <  1:1080920/100‑1 (MQ=255)                                                                                         aaaaaCCGCCCATCGGCGTCCACCGACTATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag      >  2:1896504/1‑100 (MQ=255)                                                                                         aaaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGCAGGGAACACGCAATAATGCGCCAGACACGGAGGGTAATGCGGTTAATAAAAag      >  1:920951/1‑100 (MQ=255)                                                                                          aaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAaga     <  2:603141/100‑1 (MQ=255)                                                                                          aaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAaga     >  2:257360/1‑100 (MQ=255)                                                                                           aaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAagag    <  1:638917/100‑1 (MQ=255)                                                                                            aaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGt   <  1:130268/100‑1 (MQ=255)                                                                                            aaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaagatcg   >  2:1916515/1‑95 (MQ=255)                                                                                             aCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGtt  >  2:1876656/1‑100 (MQ=255)                                                                                                 |                                                                                                CGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTT  >  NC_000913/1543658‑1543849

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: ‑