Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	344,961	T→C	N262N (AAT→AAC)	yahL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	344,961	0	T	C	95.2%	206.0 / 0.4	63	N262N (AAT→AAC)	yahL	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (35/25);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAAT  >  NC_000913/344865‑345054                                                                                                 |                                                                                              tGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAg                                                                                            <  1:159848/100‑1 (MQ=255)    tCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCAc                                                                                         <  1:1024104/100‑1 (MQ=255)       ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGtt                                                                                      >  1:624509/1‑100 (MQ=255)       ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGtt                                                                                      >  1:706570/1‑100 (MQ=255)       ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGtt                                                                                      <  1:1264543/100‑1 (MQ=255)       ggACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGtt                                                                                      >  2:830180/1‑100 (MQ=255)         aCGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTaa                                                                                    >  2:102430/1‑100 (MQ=255)          cGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTaaa                                                                                   >  2:413363/1‑100 (MQ=255)          cGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTaaa                                                                                   >  1:1074209/1‑100 (MQ=255)            tAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAAt                                                                                 >  2:499888/1‑100 (MQ=255)            tAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAAt                                                                                 >  2:151274/1‑100 (MQ=255)            tAAAAACACACAGAGAGCATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAAt                                                                                 >  2:843569/1‑100 (MQ=255)             naaaaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATc                                                                                >  1:1196306/2‑100 (MQ=255)              aaaaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATcc                                                                               >  2:904087/1‑100 (MQ=255)              aaaaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATcc                                                                               <  2:1201424/100‑1 (MQ=255)                aaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGc                                                                             >  2:1136166/1‑100 (MQ=255)                aaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGc                                                                             >  2:159529/1‑100 (MQ=255)                aaCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATAATATCCTGTAAACAAGCACGTTAAAATCCGc                                                                             >  1:1161008/1‑100 (MQ=255)                 aCATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCt                                                                            <  1:118952/100‑1 (MQ=255)                    tACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTaca                                                                         >  1:740807/1‑100 (MQ=255)                      cAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACAct                                                                       >  1:782926/1‑100 (MQ=255)                      cAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACAct                                                                       >  1:211176/1‑100 (MQ=255)                       gatcttgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACActc                                                                      <  1:814491/95‑1 (MQ=255)                            tgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACActcagatc                                                                 >  2:814485/1‑95 (MQ=255)                            tgatgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAg                                                                 <  2:159311/100‑1 (MQ=255)                              atgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGtt                                                               <  1:616529/100‑1 (MQ=255)                               tgCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTc                                                              <  2:1251722/100‑1 (MQ=255)                                    gATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGaaa                                                         <  2:296051/100‑1 (MQ=255)                                       tcaatgaatgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTTGTTCAGAAAggg                                                      <  1:843575/98‑1 (MQ=255)                                       tAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAgg                                                       <  2:988195/99‑1 (MQ=255)                                        aaaTGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAggg                                                      >  1:532464/1‑99 (MQ=255)                                        aaaTGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGt                                                     >  1:73800/1‑100 (MQ=255)                                        aaaTGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGt                                                     >  1:587291/1‑100 (MQ=255)                                         aatgaatgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGtt                                                    <  1:499892/100‑1 (MQ=255)                                         aatgaatgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGtt                                                    <  2:792307/100‑1 (MQ=255)                                           tgaatgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTcc                                                  >  1:816765/1‑100 (MQ=255)                                           tgaatgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTcc                                                  >  2:1033960/1‑100 (MQ=255)                                              atgGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCAcc                                               <  2:104689/100‑1 (MQ=255)                                                ggATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTg                                             <  1:1136176/100‑1 (MQ=255)                                                  atatCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCa                                           <  1:138723/100‑1 (MQ=255)                                                     tCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACACGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCCCCTGCAAg                                         >  1:577916/1‑99 (MQ=255)                                                       tAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAgttg                                      >  2:688898/1‑100 (MQ=255)                                                                 cAGGAATTAGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCAcc                            >  1:810096/1‑100 (MQ=255)                                                                         aGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                    <  2:874179/100‑1 (MQ=255)                                                                         aGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                    <  2:331806/100‑1 (MQ=255)                                                                         aGGTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGCTCCACCTGCAAGTTGTTATTCCACCTGAGAGTg                    >  2:1175786/1‑100 (MQ=255)                                                                          ggTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGATAGTGa                   <  2:1074199/100‑1 (MQ=255)                                                                          ggTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGa                   <  1:103870/100‑1 (MQ=255)                                                                          ggTGACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGa                   >  1:924446/1‑100 (MQ=255)                                                                             gACAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATc                <  1:110226/100‑1 (MQ=255)                                                                              aCAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCg               <  2:444630/100‑1 (MQ=255)                                                                              aCAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCg               >  1:1278620/1‑100 (MQ=255)                                                                              aCAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCg               <  2:1024528/100‑1 (MQ=255)                                                                               cAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGc              >  1:341583/1‑100 (MQ=255)                                                                               cAGATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGc              >  2:200790/1‑100 (MQ=255)                                                                                 gATATTATCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCaa            <  2:1219595/100‑1 (MQ=255)                                                                                   tattatCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATg          >  1:153496/1‑100 (MQ=255)                                                                                    attatCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGc         <  2:924438/100‑1 (MQ=255)                                                                                      tctCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATgn         <  2:1196295/96‑2 (MQ=255)                                                                                       atCCTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAgg      >  1:814835/1‑100 (MQ=255)                                                                                         ccTGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTa    >  2:1155244/1‑100 (MQ=255)                                                                                         ccTGTAAACAAGAACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAggca    >  2:891993/1‑98 (MQ=255)                                                                                           tGTAAACAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAAt  >  1:1022565/1‑100 (MQ=255)                                                                                                 |                                                                                              TGTTCTGGACGTAAAAACATACAGAGATGATGCGCGATTAAATGAATGGATATCTAATAGAGTGCAGGAATTAGGTGACAGATATTATCCTGTAAATAAGCACGTTAAAATCCGCTACACTCTTTAGTTCAGAAAGGGTTCCACCTGCAAGTTGTTATTCCACCTGAGAGTGAATCGCAATGCAGGTAAT  >  NC_000913/344865‑345054

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