Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,613,804	C→T	S6S (AGC→AGT)	b1527 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,613,804	0	C	T	98.8%	283.7 / NA	82	S6S (AGC→AGT)	b1527	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/44);  minor base G (1/0);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGCATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAT  >  NC_000913/1613716‑1613891                                                                                         |                                                                                        tttGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCa                                                                                     <  2:3964541/93‑1 (MQ=255)    gATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCa                                                                                  >  2:945406/1‑93 (MQ=255)     aTAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCaa                                                                                 <  2:4403444/93‑1 (MQ=255)     aTAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCaa                                                                                 <  2:3189712/93‑1 (MQ=255)     aTAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCaa                                                                                 <  1:4578624/93‑1 (MQ=255)      tAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAAt                                                                                >  1:3730137/1‑93 (MQ=255)       aGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAAtt                                                                               <  1:1793439/93‑1 (MQ=255)        gTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTa                                                                              >  1:1236035/1‑93 (MQ=255)        gTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTa                                                                              >  1:994955/1‑93 (MQ=255)         ttCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGCCTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTTTCAGCACTTAc                                                                             >  1:2148464/1‑93 (MQ=255)         ttCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTAc                                                                             <  1:1649631/93‑1 (MQ=255)         ttCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTAc                                                                             <  2:506742/93‑1 (MQ=255)         ttCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTAc                                                                             >  2:4625543/1‑93 (MQ=255)          tCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACa                                                                            <  2:4202271/93‑1 (MQ=255)             tGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACatat                                                                         >  2:3684647/1‑93 (MQ=255)              tcgtttAATATGCTCTATATAGTGATGTGCCGATGACTTATGACTATATGGGGCAAATATGGTTGCGCCAGTAAGTATCAGCAATTACATAtc                                                                        <  2:2653857/89‑1 (MQ=255)              gCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtc                                                                        >  2:3691681/1‑93 (MQ=255)              gCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtc                                                                        <  2:2741400/93‑1 (MQ=255)                ggTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtctc                                                                      >  2:706410/1‑93 (MQ=255)                ggTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtctc                                                                      <  1:2832840/93‑1 (MQ=255)                ggTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtctc                                                                      >  1:1650593/1‑93 (MQ=255)                 gTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATctct                                                                     >  1:1915221/1‑93 (MQ=255)                 gTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATctct                                                                     >  1:826996/1‑93 (MQ=255)                 gTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATctct                                                                     <  2:4673962/93‑1 (MQ=255)                 gTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATctct                                                                     <  2:3135807/93‑1 (MQ=255)                  ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                    >  2:4429567/1‑93 (MQ=255)                  ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                    >  2:1306717/1‑93 (MQ=255)                  ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                    >  2:885462/1‑93 (MQ=255)                   tAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTc                                                                   <  1:3206184/93‑1 (MQ=255)                   tAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTc                                                                   >  1:90842/1‑93 (MQ=255)                    aaTATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTcc                                                                  <  1:2468560/93‑1 (MQ=255)                     atatGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCt                                                                 <  2:1151944/93‑1 (MQ=255)                        tGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCtgat                                                              >  1:1767585/1‑93 (MQ=255)                        tGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCtgat                                                              <  1:993223/93‑1 (MQ=255)                          ctctATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCtgatga                                                            <  1:2960685/93‑1 (MQ=255)                            ctATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAtt                                                          <  2:264494/93‑1 (MQ=255)                             tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                         >  1:4504524/1‑93 (MQ=255)                             tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                         >  1:4130850/1‑93 (MQ=255)                             tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                         <  1:4096749/93‑1 (MQ=255)                             tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                         >  2:2080990/1‑93 (MQ=255)                               tataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTc                                                       >  2:381288/1‑93 (MQ=255)                               tataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTc                                                       >  1:2737993/1‑93 (MQ=255)                                 taGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTccc                                                     <  1:1421135/93‑1 (MQ=255)                                    tGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgta                                                  >  2:649488/1‑93 (MQ=255)                                      aTGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgtacg                                                >  1:429704/1‑93 (MQ=255)                                      aTGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgtacg                                                >  2:3634995/1‑93 (MQ=255)                                         ttCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTaa                                             >  1:2857859/1‑93 (MQ=255)                                           ccGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAata                                           <  2:4239753/93‑1 (MQ=255)                                              aTGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTg                                        >  1:2790008/1‑93 (MQ=255)                                                 aCTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCAc                                     <  1:175993/93‑1 (MQ=255)                                                 aCTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGAATCAGCAATTACATATCTCTTCCTGATGATTTTTCCGTACGTAATATAGCAc                                     >  2:3010060/1‑93 (MQ=255)                                                  cTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCAcc                                    >  2:2380508/1‑93 (MQ=255)                                                  cTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGGATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCAcc                                    >  2:1805053/1‑93 (MQ=255)                                                     aTGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCa                                 <  2:3577233/93‑1 (MQ=255)                                                      tGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAg                                <  2:2123462/93‑1 (MQ=255)                                                          tataTGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaa                            >  2:139687/1‑93 (MQ=255)                                                            tatGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaaaa                          >  1:3943483/1‑93 (MQ=255)                                                               ggggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaa                       <  1:54941/93‑1 (MQ=255)                                                                gggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaag                      >  2:609206/1‑93 (MQ=255)                                                                 ggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGt                     <  1:102253/93‑1 (MQ=255)                                                                   cAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGttt                   <  2:1462472/93‑1 (MQ=255)                                                                    aaaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGtttt                  <  2:4647586/93‑1 (MQ=255)                                                                     aaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTa                 >  2:3945315/1‑93 (MQ=255)                                                                     aaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTa                 <  1:1619531/93‑1 (MQ=255)                                                                     aaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTa                 <  1:3059536/93‑1 (MQ=255)                                                                      atatGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTaa                >  2:2598048/1‑93 (MQ=255)                                                                       tatGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTaaa               <  2:1427484/93‑1 (MQ=255)                                                                       tatGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTaaa               <  1:4089240/93‑1 (MQ=255)                                                                       tatGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTaaa               >  1:2051920/1‑93 (MQ=255)                                                                        atGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAg              <  2:3745239/93‑1 (MQ=255)                                                                        atGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAg              <  2:3147822/93‑1 (MQ=255)                                                                        atGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAg              <  2:75979/93‑1 (MQ=255)                                                                          ggTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGa            <  1:750676/93‑1 (MQ=255)                                                                           gTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGAt           <  1:910673/93‑1 (MQ=255)                                                                           gTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGAt           <  1:3852699/93‑1 (MQ=255)                                                                           gTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATAATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGAt           >  1:2517253/1‑93 (MQ=255)                                                                               cGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTa       <  2:3132993/93‑1 (MQ=255)                                                                                gCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTAt      <  1:3992694/93‑1 (MQ=255)                                                                                 ccAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTAtt     >  1:531982/1‑93 (MQ=255)                                                                                  cAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTg    <  2:647597/93‑1 (MQ=255)                                                                                  cAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTg    >  2:497738/1‑93 (MQ=255)                                                                                   agtaagtaTCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGa   >  1:3033636/1‑93 (MQ=255)                                                                                    gTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAt  <  1:1862429/93‑1 (MQ=255)                                                                                    gTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAt  <  2:2111600/93‑1 (MQ=255)                                                                                         |                                                                                        TTTGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGCATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAT  >  NC_000913/1613716‑1613891

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