Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,685,709	G→A	intergenic (‑202/‑126)	glyA ← / → hmp	serine hydroxymethyltransferase/nitric oxide dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,685,709	0	G	A	100.0%	180.3 / NA	50	intergenic (‑202/‑126)	glyA/hmp	serine hydroxymethyltransferase/nitric oxide dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/24);  total (26/24)

GCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTGATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTT  >  NC_000913/2685573‑2685840                                                                                                                                         |                                                                                                                                    gCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATg                                                                                                                                  >  2:268851/1‑140 (MQ=255)    aaCAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCaa                                                                                                                               >  1:839306/1‑140 (MQ=255)         tataACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAAtttt                                                                                                                          >  1:484733/1‑140 (MQ=255)          ataACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAAttttt                                                                                                                         <  1:1184962/140‑1 (MQ=255)           taACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTc                                                                                                                        >  2:703127/1‑140 (MQ=255)                   tCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTtcatc                                                                                                                >  1:261222/1‑140 (MQ=255)                   tCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTtcatc                                                                                                                >  2:382902/1‑140 (MQ=255)                   tCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTtcatc                                                                                                                >  2:667502/1‑140 (MQ=255)                   tCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTtcatc                                                                                                                >  2:613625/1‑140 (MQ=255)                        gtgCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACAtt                                                                                                           >  2:71485/1‑140 (MQ=255)                          gCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTtc                                                                                                                   <  2:704772/130‑1 (MQ=255)                          gCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTtc                                                                                                                   >  1:704772/1‑130 (MQ=255)                             aTAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACAttctttc                                                                                                      >  2:791192/1‑140 (MQ=255)                               aaCAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTTCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTg                                                                                                    <  2:1085903/140‑1 (MQ=255)                                    gTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGaaaaa                                                                                               <  1:613625/140‑1 (MQ=255)                                    gTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGaaaaa                                                                                               <  1:791192/140‑1 (MQ=255)                                        tGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACAcc                                                                                           >  2:689967/1‑140 (MQ=255)                                            aaaGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAg                                                                                       <  2:261222/140‑1 (MQ=255)                                              aGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGaa                                                                                                  >  2:466249/1‑127 (MQ=255)                                              aGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGaa                                                                                                  <  1:466249/127‑1 (MQ=255)                                              aGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGaa                                                                                     >  2:533676/1‑140 (MQ=255)                                                gAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAAcc                                                                                   <  1:667502/140‑1 (MQ=255)                                                            aCGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCa                                                                       <  2:484733/140‑1 (MQ=255)                                                               aTTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAggg                                                                    <  1:533676/140‑1 (MQ=255)                                                                          gCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAggg                                                                    <  1:721628/129‑1 (MQ=255)                                                                          gCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAggg                                                                    >  2:721628/1‑129 (MQ=255)                                                                            tACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATaa                                                       >  2:1000004/1‑140 (MQ=255)                                                                              cGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGa                                                     >  1:931586/1‑140 (MQ=255)                                                                              cGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGa                                                     >  2:377754/1‑140 (MQ=255)                                                                              cGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGa                                                     <  1:818114/140‑1 (MQ=255)                                                                               gCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGAt                                                    <  2:839306/140‑1 (MQ=255)                                                                                        tGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTgag                                           >  1:129070/1‑140 (MQ=255)                                                                                        tGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTgag                                           >  1:521843/1‑140 (MQ=255)                                                                                         ggAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTgaga                                          <  2:931586/140‑1 (MQ=255)                                                                                           agaATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATa                                        <  1:1217556/140‑1 (MQ=255)                                                                                            gaATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATAc                                       <  1:377754/140‑1 (MQ=255)                                                                                               tAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCAt                                                >  1:811966/1‑128 (MQ=255)                                                                                               tAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCAt                                                <  2:811966/128‑1 (MQ=255)                                                                                               tAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATc                                    >  1:871584/1‑140 (MQ=255)                                                                                               tAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATc                                    >  2:546340/1‑140 (MQ=255)                                                                                                  agagCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAAt                                 <  1:546340/140‑1 (MQ=255)                                                                                                    agCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATAc                                       >  2:193057/1‑132 (MQ=255)                                                                                                    agCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATAc                                       <  1:193057/132‑1 (MQ=255)                                                                                                       ttGCAACCGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGa                            >  2:655492/1‑140 (MQ=255)                                                                                                             ccGGAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaa                      <  1:296947/140‑1 (MQ=255)                                                                                                                gAAACGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCaaaaaa                   >  2:455984/1‑140 (MQ=255)                                                                                                                  aaCGGATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAg                 <  2:871584/140‑1 (MQ=255)                                                                                                                     ggATTTCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGaa              <  1:455984/140‑1 (MQ=255)                                                                                                                         ttCTTTTCAGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGAcc          <  1:268851/140‑1 (MQ=255)                                                                                                                                 aGGTTTGTAATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCtt  <  1:1094769/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    GCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTCTTTTCAGGTTTGTGATGCAAATTTTTCACTTCATCACATTCTTTCTGAAAAACACCAAAGAACCATTTACATTGCAGGGCTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTT  >  NC_000913/2685573‑2685840

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