Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	364,024	C→T	R761H (CGC→CAC)	lacZ ←	beta‑D‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,024	0	C	T	94.7%	175.1 / 0.9	57	R761H (CGC→CAC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (34/20);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCG  >  NC_000913/363929‑364117                                                                                                |                                                                                              cGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGtt                                                                                           <  2:1015634/100‑1 (MQ=255)  gTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTa                                                                                          >  1:1097719/1‑100 (MQ=255)  gTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTa                                                                                          >  1:1399573/1‑100 (MQ=255)     aTCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTaa                                                                                         <  2:179623/98‑1 (MQ=255)     aTCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAAt                                                                                       >  1:312112/1‑100 (MQ=255)       ccAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTg                                                                                     <  2:1016550/100‑1 (MQ=255)       ccAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTg                                                                                     <  1:1197877/100‑1 (MQ=255)        cAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTgn                                                                                    <  2:636657/100‑2 (MQ=255)          gCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCa                                                                                  >  1:513045/1‑100 (MQ=255)             gTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACg                                                                               <  2:756556/100‑1 (MQ=255)                cACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCtt                                                                            >  2:1171394/1‑100 (MQ=255)                  cGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGGCAACGCttat                                                                          >  2:931882/1‑100 (MQ=255)                    ggTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatta                                                                        <  1:1519784/100‑1 (MQ=255)                     gTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAc                                                                       <  1:1426743/100‑1 (MQ=255)                        aaCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCa                                                                    >  2:475073/1‑100 (MQ=255)                         aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAg                                                                   <  1:587282/100‑1 (MQ=255)                          cagATCGCGCAGCAGCGTCAGCAGTTGTTTTGTATCGCCAATACACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGc                                                                  >  1:1066154/3‑100 (MQ=255)                               cgcgCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGAt                                                             >  1:1100969/1‑100 (MQ=255)                                 cgcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGCTGc                                                           >  2:689516/1‑100 (MQ=255)                                   cAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCaa                                                         >  2:1057724/1‑100 (MQ=255)                                       ggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAAt                                                     >  1:1146424/1‑100 (MQ=255)                                       ggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAAt                                                     >  1:1473519/1‑100 (MQ=255)                                          gTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCa                                                  >  1:1232284/1‑100 (MQ=255)                                            cagcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAtt                                                <  2:1232279/100‑1 (MQ=255)                                              gcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTc                                              >  2:1443489/1‑100 (MQ=255)                                                tcttgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGc                                            <  1:275152/98‑1 (MQ=255)                                                 gttgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCt                                           >  2:422068/1‑100 (MQ=255)                                                 gttgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCt                                           <  1:1201846/100‑1 (MQ=255)                                                 gttgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCt                                           >  1:199977/1‑100 (MQ=255)                                                   tgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGcaga                                         >  2:275152/1‑97 (MQ=255)                                                    gttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggt                                        <  1:128419/100‑1 (MQ=255)                                                     ttttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggtg                                       >  1:871774/1‑100 (MQ=255)                                                      tttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggtgg                                      >  1:688301/1‑100 (MQ=255)                                                      tttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggtgg                                      <  1:93846/100‑1 (MQ=255)                                                        tttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTc                                    >  2:721731/1‑100 (MQ=255)                                                          tATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGGTTATTACCCAGCTCGATGCAAAAATCCAGTTCGCTGGTGGTCa                                   <  2:1389211/99‑1 (MQ=255)                                                            tCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGAt                                >  2:783504/1‑100 (MQ=255)                                                                cAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCgg                            >  2:490933/1‑100 (MQ=255)                                                                  aTCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGa                          <  2:287750/100‑1 (MQ=255)                                                                    ccACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATg                        >  2:812677/1‑100 (MQ=255)                                                                     cacaTCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATgg                       >  1:636988/1‑100 (MQ=255)                                                                        aTCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGt                    >  1:550783/1‑100 (MQ=255)                                                                         tCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTg                   <  1:1174973/100‑1 (MQ=255)                                                                          cTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTgg                  <  1:783508/100‑1 (MQ=255)                                                                            gtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                 >  2:1463553/1‑99 (MQ=255)                                                                              gaaagaaaGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcg              >  2:319278/1‑100 (MQ=255)                                                                               aaagaaagCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcgc             >  1:512976/1‑100 (MQ=255)                                                                                aagaaagCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcgcg            >  1:1433549/1‑100 (MQ=255)                                                                                aagaaagCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcgcg            <  1:330576/100‑1 (MQ=255)                                                                                aagaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTGAGATGCGGGATGGCGTGGGAcgcg            >  1:1175504/1‑100 (MQ=255)                                                                                 agaaagCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACgcgg           >  1:892850/1‑100 (MQ=255)                                                                                     agCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCggg       >  1:1196570/1‑100 (MQ=255)                                                                                      gCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCgggg      >  2:1467329/1‑100 (MQ=255)                                                                                       ccTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGa     >  2:78118/1‑100 (MQ=255)                                                                                        cTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAg    <  1:1086021/100‑1 (MQ=255)                                                                                          gACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCg  >  2:932774/1‑100 (MQ=255)                                                                                          gACTGGTGGTTAAATTGCCAACGCTTATTAACCAGCTAGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGCCGCGGCGGGGAGc   <  1:316985/99‑1 (MQ=255)                                                                                                |                                                                                              CGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCG  >  NC_000913/363929‑364117

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