Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,810,941	T→C	G76G (GGA→GGG)	b1730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,810,941	0	T	C	98.7%	230.2 / ‑4.2	76	G76G (GGA→GGG)	b1730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (50/24);  total (51/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 = 2.98e-01

CGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCTCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTCACAT  >  NC_000913/1810855‑1811028                                                                                       |                                                                                        cGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGc                                                                                   >  2:2327506/1‑93 (MQ=255)    aTTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCttt                                                                                <  1:5300428/93‑1 (MQ=255)     tttCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTg                                                                               >  1:7871017/1‑93 (MQ=255)      ttCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGa                                                                              <  2:354412/93‑1 (MQ=255)      ttCACTACCATAAAGATTATTAATTAATACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGa                                                                              >  2:1554809/1‑93 (MQ=255)         aCTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGc                                                                           <  1:5723259/93‑1 (MQ=255)          cTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGca                                                                          >  2:3585600/1‑93 (MQ=255)             ccATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGcatca                                                                       <  2:6922917/93‑1 (MQ=255)             ccATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGcatca                                                                       <  2:6922918/93‑1 (MQ=255)             acATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGcatca                                                                       <  2:802161/92‑1 (MQ=255)               aTAAAGATTATTAATTATTACTCCATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCACCTGGCTTTGAGGCATCATa                                                                     <  2:5415132/93‑1 (MQ=255)                tAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAg                                                                    >  2:4153793/1‑93 (MQ=255)                 aaaGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGt                                                                   <  1:743548/93‑1 (MQ=255)                  aaGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGtt                                                                  >  2:3871100/1‑93 (MQ=255)                    gATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAACACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTaa                                                                >  1:3856372/1‑93 (MQ=255)                     attattaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAAt                                                               <  2:1234135/93‑1 (MQ=255)                     attattaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAAt                                                               >  2:997077/1‑93 (MQ=255)                       tattaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGt                                                             >  2:3990344/1‑93 (MQ=255)                        attaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGtt                                                            >  2:1113231/1‑93 (MQ=255)                         ttaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTa                                                           >  2:973479/1‑93 (MQ=255)                               attaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCTCCTGGCTTTGAGGCATCATAGTTAATGTTAGACata                                                     >  2:4121435/1‑93 (MQ=255)                                 taCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACAtata                                                   <  2:3070704/93‑1 (MQ=255)                                 taCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACAtata                                                   <  1:5813388/93‑1 (MQ=255)                                 taCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACAtata                                                   >  2:4776054/1‑93 (MQ=255)                                   ctctATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAgg                                                 <  1:7582841/93‑1 (MQ=255)                                    tctATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGc                                                <  1:903381/93‑1 (MQ=255)                                      tataTTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCaa                                              >  1:255529/1‑93 (MQ=255)                                       atatTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAAt                                             <  1:1512542/93‑1 (MQ=255)                                       ataaTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTTCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAAt                                             <  1:1383929/93‑1 (MQ=255)                                        tatTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATc                                            >  1:7524281/1‑93 (MQ=255)                                        tatTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATc                                            >  2:1262348/1‑93 (MQ=255)                                         atTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCa                                           >  2:5364434/1‑93 (MQ=255)                                         atTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCa                                           >  1:4059004/1‑93 (MQ=255)                                          tttAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCat                                          >  2:6849040/1‑93 (MQ=255)                                           ttAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCata                                         >  1:4350273/1‑93 (MQ=255)                                           ttAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCata                                         >  2:462434/1‑93 (MQ=255)                                           ttAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCata                                         >  2:2781040/1‑93 (MQ=255)                                           ttAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGATAATGTTAGACATATAGGCAATCata                                         >  2:462433/1‑93 (MQ=255)                                             aaCATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCAtata                                       >  2:2810135/1‑93 (MQ=255)                                               cATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAAt                                     <  2:4412094/93‑1 (MQ=255)                                               cATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAAt                                     <  2:1220517/93‑1 (MQ=255)                                                aTTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAAtt                                    >  2:2761244/1‑93 (MQ=255)                                                aTTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAAtt                                    >  2:7255704/1‑93 (MQ=255)                                                 ttGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTa                                   >  2:4673802/1‑93 (MQ=255)                                                 ttGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTa                                   >  1:4091645/1‑93 (MQ=255)                                                  tGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTaa                                  <  2:175750/93‑1 (MQ=255)                                                   gTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTaaa                                 >  2:3036979/1‑93 (MQ=255)                                                      cTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGtt                              >  2:2972221/1‑93 (MQ=255)                                                      cTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGtt                              >  1:4043594/1‑93 (MQ=255)                                                       tGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGttt                             >  2:2322902/1‑93 (MQ=255)                                                        gCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTg                            >  1:4757235/1‑93 (MQ=255)                                                         cAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGa                           <  1:1426098/93‑1 (MQ=255)                                                           gATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGAtt                         >  2:5752303/1‑93 (MQ=255)                                                            atatGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTc                        >  1:1348147/1‑93 (MQ=255)                                                             tatGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCt                       >  2:2485446/1‑93 (MQ=255)                                                               tGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGc                     >  1:5334576/1‑93 (MQ=255)                                                               tGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGc                     >  1:5567549/1‑93 (MQ=255)                                                               tGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGc                     >  2:1111676/1‑93 (MQ=255)                                                               tGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGc                     >  2:1204401/1‑93 (MQ=255)                                                                gTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCt                    >  2:678050/1‑93 (MQ=255)                                                                 tATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTa                   >  2:3324066/1‑93 (MQ=255)                                                                 tATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTa                   <  1:5626963/93‑1 (MQ=255)                                                                 tATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTa                   <  2:5376194/93‑1 (MQ=255)                                                                  aTCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTAt                  >  2:1088931/1‑93 (MQ=255)                                                                  aTCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTAt                  >  2:5545946/1‑93 (MQ=255)                                                                  aTCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTAt                  <  1:4334469/93‑1 (MQ=255)                                                                   tCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTAtt                 >  1:3346663/1‑93 (MQ=255)                                                                   tCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTAtt                 >  2:6863698/1‑93 (MQ=255)                                                                    cAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTAttt                <  2:6170235/93‑1 (MQ=255)                                                                       aaCTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCAt             >  1:7621637/1‑93 (MQ=255)                                                                       aaCTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCAt             >  2:1125181/1‑93 (MQ=255)                                                                        aCTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATc            >  1:84029/1‑93 (MQ=255)                                                                          tAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGc          >  1:380949/1‑93 (MQ=255)                                                                           aaGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGcg         >  1:6814590/1‑92 (MQ=255)                                                                                 tGGCGCCCCTGGCGTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTcaca   <  1:4282351/93‑1 (MQ=255)                                                                                  ggCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTCACAt  <  2:2404865/93‑1 (MQ=255)                                                                                       |                                                                                        CGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCTCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTCACAT  >  NC_000913/1810855‑1811028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑