Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	4,586,635	G→A	14.3%	P365L (CCG→CTG)	hsdM ←	DNA methylase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	4,586,635	0	G	A	14.3%	190.1 / 19.8	70	P365L (CCG→CTG)	hsdM	DNA methylase M
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base A (4/6);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

GCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCG  >  AP009048/4586565‑4586703                                                                       |                                                                     gCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGAATCg                                                                      >  1:3274249/1‑71 (MQ=255)  ccGAATCTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCgg                                                                     <  1:4525880/71‑1 (MQ=255)  ccGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCg                                                                      <  1:1045263/70‑1 (MQ=255)   cGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCgg                                                                     >  1:3612849/1‑70 (MQ=255)    gAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGTTTCGGGt                                                                   <  1:839386/71‑1 (MQ=255)     aaaCTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGtt                                                                  <  1:1587825/71‑1 (MQ=255)     aaaCTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGtt                                                                  >  1:1551108/1‑71 (MQ=255)     aaaCTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGtt                                                                  >  1:4349585/1‑71 (MQ=255)           ggCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGcc                                                              <  1:1074727/69‑1 (MQ=255)             cATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCgg                                                                     >  1:2051794/1‑60 (MQ=255)                atTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTcc                                                       >  1:1824533/1‑71 (MQ=255)                  tGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCggg                                                                    >  1:2983241/1‑56 (MQ=255)                  tGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCggg                                                                    >  1:3961565/1‑56 (MQ=255)                  tGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCAGCGTccc                                                      <  1:4267150/70‑1 (MQ=255)                  tGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGt                                                         <  1:4566061/67‑1 (MQ=255)                  tGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCt                                                     >  1:1045606/1‑71 (MQ=255)                   ggTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCt                                                     <  1:1345941/70‑1 (MQ=255)                   ggTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCt                                                     <  1:1354987/70‑1 (MQ=255)                    gTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCttt                                                   >  1:2844762/1‑71 (MQ=255)                    gTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCtt                                                    <  1:4255832/70‑1 (MQ=255)                    gTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCtt                                                    <  1:2741607/70‑1 (MQ=255)                      aCGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCggg                                                                    <  1:1373643/52‑1 (MQ=255)                       cGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTgg                                                 >  1:74933/1‑70 (MQ=255)                          aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTaaa                                             <  1:4101916/71‑1 (MQ=255)                            gTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGtt                                                                  <  1:64637/48‑1 (MQ=255)                             tCATACACCCACACATCATCGGTACAGTTCTTATCCTTATTCGGGTTCGCCACCGTCCCTTTGGTAAAGaa                                          <  1:3667588/71‑1 (MQ=255)                             tCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGaa                                          <  1:2965523/71‑1 (MQ=255)                               aTACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCAc                                                            >  1:717490/1‑51 (MQ=255)                               aTACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGaa                                          >  1:1708185/1‑69 (MQ=255)                                tACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGc                                                               <  1:753303/47‑1 (MQ=255)                                 acacCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCtt                                                    >  1:4172121/1‑57 (MQ=255)                                  cacCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCggc                                                              <  1:3787781/46‑3 (MQ=255)                                  cacCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGaa                                          >  1:3235056/1‑66 (MQ=255)                                   acCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTc                                                                 >  1:935022/1‑42 (MQ=255)                                   acCCACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCAcc                                                           >  1:4096890/1‑48 (MQ=255)                                   acCCACACATCATCGGTACAGTTCTTAACCTGATTCGGGTTCGCCa                                                             >  1:2248124/1‑46 (MQ=255)                                     ccACACATCATCGGTACAGTTCTTATCCTGATTCgc                                                                     >  1:2467801/1‑35 (MQ=255)                                     ccACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGt                                  <  1:2885498/71‑1 (MQ=255)                                     ccACACATCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGt                                  <  1:746677/71‑1 (MQ=255)                                     ccACACATCATCGGTACAGTTCTTATACTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACa                                        >  1:3756641/1‑65 (MQ=255)                                      cacacaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCa                                                             >  1:4535044/1‑43 (MQ=255)                                       acacaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCa                                                             >  1:3581979/1‑42 (MQ=255)                                       acacaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTaa                                              <  1:570976/57‑1 (MQ=255)                                       acacaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACa                                        >  1:3335965/1‑63 (MQ=255)                                       acacaTCATCGGTACAGTTCTTATCCTGATTCGGGTACGCCACCGTCCCTTTGGTAAAGAACAGCACg                                   >  1:3884243/1‑68 (MQ=255)                                       acacaTCATCGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTg                                >  1:1903466/1‑71 (MQ=255)                                        cacaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTgg                               <  1:161101/71‑1 (MQ=255)                                         acaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGt                              <  1:3336889/71‑1 (MQ=255)                                          catcatCGGTACAGTTCTTATCCTGATTCGGGTTc                                                                 >  1:3217258/1‑35 (MQ=255)                                          catcatCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTc                             <  1:2656420/71‑1 (MQ=255)                                           atcatcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCt                            <  1:1386693/71‑1 (MQ=255)                                           atcatcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCt                            <  1:1582270/71‑1 (MQ=255)                                           atcatcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCt                            <  1:1587617/71‑1 (MQ=255)                                           atcatcGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTg                                <  1:2243010/67‑1 (MQ=255)                                            tcatcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTaaa                                             >  1:2079748/1‑53 (MQ=255)                                            tcatcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCtt                           >  1:1545450/1‑71 (MQ=255)                                             catcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTaa                                              >  1:45779/1‑51 (MQ=255)                                              atcGGTACAGTTCTTATCCTGATTCGGGTTCGCCAc                                                            >  1:4262767/1‑36 (MQ=255)                                              atcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCa                         >  1:1284452/1‑71 (MQ=255)                                              atcGGTACAGTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCa                         >  1:299202/1‑71 (MQ=255)                                               tcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCAc                        <  1:4500749/71‑1 (MQ=255)                                                cGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAg                                            >  1:397909/1‑50 (MQ=255)                                                 ggTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACg                       <  1:2054502/70‑1 (MQ=255)                                                  gTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                                       >  1:2113645/1‑53 (MQ=255)                                                   tACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACg                       <  1:3550033/68‑1 (MQ=255)                                                     cAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCt                   <  1:3317148/70‑1 (MQ=255)                                                     cAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCt                   <  1:3680226/70‑1 (MQ=255)                                                     cAGTTCTTATCCCGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCt                   >  1:1169085/1‑70 (MQ=255)                                                       gTTCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGTTAAAGAACAGCACGt                                  <  1:1321738/53‑1 (MQ=255)                                                        ttCTTATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCt                   <  1:423335/67‑1 (MQ=255)                                                          cTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCa                         >  1:2992385/1‑59 (MQ=255)                                                           ttATCCTGATTCAGGTTCGCCACCGTCCCTTTGGTAAAg                                            <  1:3850954/39‑1 (MQ=255)                                                             aTCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                                       <  1:1575834/42‑1 (MQ=255)                                                             aTCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTg                                <  1:92804/49‑1 (MQ=255)                                                                  gATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATa     >  1:1345569/1‑71 (MQ=255)                                                                   aTTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATa     <  1:1786139/70‑1 (MQ=255)                                                                   aTTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGAAAAAAAt      <  1:451709/69‑1 (MQ=255)                                                                   aTTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGGCTTCACGCCCTGAGCGTAAAAAAt      <  1:693740/69‑1 (MQ=255)                                                                      cGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCg  <  1:2672324/70‑1 (MQ=255)                                                                       |                                                                     GCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCG  >  AP009048/4586565‑4586703

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