Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	681,881	T→A	I181F (ATC→TTC)	nagD ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	681,881	0	T	A	97.6%	158.0 / ‑4.2	42	I181F (ATC→TTC)	nagD	UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (21/20);  total (22/20)
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.82e-01

ACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGATCCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGCACACAACGCGCC  >  REL606/681810‑681953                                                                        |                                                                         aCGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCAc                                                                      >  2:606380/1‑76 (MQ=255) aCGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCAc                                                                      <  1:875520/76‑1 (MQ=255)     aGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGAc                                                                  <  1:868403/76‑1 (MQ=255)      ggTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACt                                                                 <  1:1351746/76‑1 (MQ=255)        ttATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCg                                                               <  2:744726/76‑1 (MQ=255)          aTCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGATCCACGGACTCGGt                                                             >  1:380233/1‑76 (MQ=255)          aTCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGt                                                             >  2:389250/1‑76 (MQ=255)             gCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTa                                                          <  1:1059549/76‑1 (MQ=255)             gCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTa                                                          <  2:576958/76‑1 (MQ=255)                gACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCa                                                       <  1:1332227/76‑1 (MQ=255)                   aaTCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACa                                                    >  2:898083/1‑76 (MQ=255)                         cGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAAc                                              >  1:737186/1‑76 (MQ=255)                           tttCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACgg                                            >  1:976026/1‑76 (MQ=255)                             tcttcCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGtt                                          >  1:1521350/1‑76 (MQ=255)                              cttcCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGttt                                         <  1:418351/76‑1 (MQ=255)                               ttcCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTa                                        >  1:1086502/1‑76 (MQ=255)                                 ccGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACg                                      <  1:579298/76‑1 (MQ=255)                                  cGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGc                                     <  1:534100/76‑1 (MQ=255)                                    aaTGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGccc                                   <  2:1456904/76‑1 (MQ=255)                                      tGCGCCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCgg                                 <  1:1390307/76‑1 (MQ=255)                                        cgcCTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGgag                               >  2:1313532/1‑76 (MQ=255)                                           cTGCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGAtt                            <  2:974091/76‑1 (MQ=255)                                             gCATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGAtttt                          >  1:1470795/1‑76 (MQ=255)                                              cATTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGAttttt                         >  2:876339/1‑76 (MQ=255)                                               aTTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGAtttttt                        >  1:844134/1‑76 (MQ=255)                                               aTTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGAtttttt                        >  2:461150/1‑76 (MQ=255)                                               aTTTTGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGAtttttt                        <  2:211542/76‑1 (MQ=255)                                                 tttGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCg                      <  2:995273/76‑1 (MQ=255)                                                  ttGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGa                     >  1:995395/1‑76 (MQ=255)                                                  ttGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGa                     >  1:120079/1‑76 (MQ=255)                                                   tGTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGAt                    >  2:406320/1‑76 (MQ=255)                                                    gTTTAATGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATg                   >  2:74306/1‑76 (MQ=255)                                                        aaTGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTg               >  1:1348688/1‑76 (MQ=255)                                                        aaTGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTg               <  1:846298/76‑1 (MQ=255)                                                        aaTGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTg               >  1:1040068/1‑76 (MQ=255)                                                        aaTGCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTg               >  1:486774/1‑76 (MQ=255)                                                           gCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGcac            >  2:1029556/1‑76 (MQ=255)                                                           gCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGcac            >  2:741243/1‑76 (MQ=255)                                                           gCCGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGcac            <  1:1410910/76‑1 (MQ=255)                                                           gCCGCGCGGATGAACCACGGACTCGGTTTACCAACAGAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGcac            <  2:1536990/76‑1 (MQ=255)                                                            ccGCGCGGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGcaca           <  1:172841/76‑1 (MQ=255)                                                             cgcgcgGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGcacac          <  1:620899/76‑1 (MQ=255)                                                                gcgGATGAACCACGGACTCGGTTTACCAACATAGAACGGTTTACGGCCGGAGATTTTTTCGATGCCTGCACACAAc       <  1:71483/76‑1 (MQ=255)                                                                     tgaaCCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGCACACAACGCGcc  >  2:1355643/1‑76 (MQ=255)                                                                        |                                                                         ACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCCGCGCGGATGATCCACGGACTCGGTTTACCAACATAGAACGGTTTACGCCCGGAGATTTTTTCGATGCCTGCACACAACGCGCC  >  REL606/681810‑681953

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