Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,644,444	T→G	6.6%	V230V (GTA→GTC)	nlpB ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,644,444	0	T	G	6.6%	168.2 / 6.9	61	V230V (GTA→GTC)	nlpB	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base G (2/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCAAACCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAGACCG  >  Exported/1644378‑1644513                                                                   |                                                                      ccAAACCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGt                                                                   >  1:1025283/1‑71 (MQ=255)  cAAACCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTc                                                                  <  1:121824/71‑1 (MQ=255)   aaaCCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTc                                                                  >  1:4098/1‑70 (MQ=255)    aaCCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCa                                                                <  1:1089862/71‑1 (MQ=255)    aaCCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCa                                                                >  1:1622030/1‑71 (MQ=255)       cacaTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTATGGACGTACAAAg                                                             <  1:2475901/71‑1 (MQ=255)        acaTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgt                                                            >  1:1096015/1‑71 (MQ=255)         caTTGAACGGCCCGCGTACGACCATCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtg                                                           <  1:215284/71‑1 (MQ=255)         caTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtg                                                           <  1:974574/71‑1 (MQ=255)         caTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtg                                                           <  1:282944/71‑1 (MQ=255)         caTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtg                                                           <  1:2276920/71‑1 (MQ=255)          aTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtgg                                                          >  1:1348224/1‑71 (MQ=255)          aTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtgg                                                          >  1:1307264/1‑71 (MQ=255)          aTTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtg                                                           >  1:972442/1‑70 (MQ=255)             gAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtggtg                                                        <  1:2561662/70‑1 (MQ=255)              aaCGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtggtg                                                        >  1:2880668/1‑69 (MQ=255)              aaCGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGtggtgg                                                       <  1:772668/70‑1 (MQ=255)              aaCGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGtggtgg                                                       <  1:873049/70‑1 (MQ=255)                 ggCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGc                                                   <  1:2122714/71‑1 (MQ=255)                    ccGCGTACGACCAGCATTGTTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACg                                                <  1:123862/71‑1 (MQ=255)                     cgcgTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGa                                               <  1:813563/71‑1 (MQ=255)                     cgcgTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGa                                               <  1:763239/71‑1 (MQ=255)                     cgcgTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGa                                               <  1:2490558/71‑1 (MQ=255)                     cgcgTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGa                                               <  1:1491880/71‑1 (MQ=255)                       cgTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGAtt                                             >  1:897294/1‑71 (MQ=255)                          aCGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTg                                          >  1:2070982/1‑71 (MQ=255)                           cGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTgc                                         >  1:2999412/1‑71 (MQ=255)                           cGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTgc                                         >  1:1050543/1‑71 (MQ=255)                              ccAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCa                                      <  1:26939/71‑1 (MQ=255)                                 gCATTGGTAAACCGGTGTCATCAGCTGCACTATGGACGTACAAAg                                                             <  1:2984836/45‑1 (MQ=25)                                  cATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtg                                                           <  1:3031863/46‑1 (MQ=255)                                   aTTGGTAAACCGGTGTCATCAGCTGCACTATGGACGTACAAAgtg                                                           >  1:418479/1‑45 (MQ=21)                                   aTTGGTAAACCGGTGTCATCAGCTGCACTATGGACGTACAAAgtg                                                           >  1:2844247/1‑45 (MQ=21)                                      ggTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtgg                                                          >  1:1092739/1‑43 (MQ=255)                                      ggTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtgg                                                          >  1:922105/1‑43 (MQ=255)                                      ggTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtgg                                                          >  1:576944/1‑43 (MQ=255)                                      ggTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTcg                               >  1:1415516/1‑70 (MQ=255)                                      ggTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGGGGAGGCACGATTTTGCGCAGCGTTcg                               <  1:697026/70‑1 (MQ=255)                                       gTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtggtg                                                        <  1:1300746/44‑1 (MQ=255)                                        tAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAgtggtg                                                        <  1:112498/43‑1 (MQ=255)                                         aaaCCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGa                                                      >  1:2216531/1‑44 (MQ=255)                                          aaCCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGa                                                      >  1:295652/1‑43 (MQ=255)                                                tgtCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTg                    <  1:2380876/71‑1 (MQ=255)                                                 gtCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTgg                   >  1:2946515/1‑71 (MQ=255)                                                    atcaGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTc                                >  1:561064/1‑55 (MQ=255)                                                    atcaGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAg                <  1:2787568/71‑1 (MQ=255)                                                         cTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGAGGCAGATTTa           <  1:985194/71‑1 (MQ=255)                                                           gCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATc         >  1:2568991/1‑71 (MQ=255)                                                            cACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATcc        <  1:1888968/71‑1 (MQ=255)                                                            cACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATcc        <  1:2946855/71‑1 (MQ=255)                                                            cACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATc         <  1:816481/70‑1 (MQ=255)                                                            cACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATc         <  1:1106688/70‑1 (MQ=255)                                                             aCTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCa       <  1:602418/71‑1 (MQ=255)                                                              cTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATc         >  1:989862/1‑68 (MQ=255)                                                              cTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAg      <  1:1655757/71‑1 (MQ=255)                                                               tttGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTc                        <  1:3003637/52‑1 (MQ=255)                                                                ttGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGc                                    >  1:2915630/1‑39 (MQ=255)                                                                ttGTACATCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAGAc    >  1:2474162/1‑71 (MQ=255)                                                                 tGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAGa     <  1:2697336/69‑1 (MQ=255)                                                                 tGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAGa     <  1:1972212/69‑1 (MQ=255)                                                                 tGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAGAcc   >  1:993025/1‑71 (MQ=255)                                                                   tACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGGCAGTGGCAGATTTATCCAGACcg  <  1:2693183/70‑1 (MQ=255)                                                                   |                                                                      CCAAACCACATTGAACGGCCCGCGTACGACCAGCATTGGTAAACCGGTGTCATCAGCTGCACTTTGTACGTCCATAGTGGTGGAGGCACGATTTTGCGCAGCGTTCGCGGCGTCAGTGGCAGATTTATCCAGACCG  >  Exported/1644378‑1644513

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