Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	699,862	A→G	A155A (GCT→GCC)	umpH ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	699,862	0	A	G	100.0%	227.1 / NA	63	A155A (GCT→GCC)	umpH	UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/31);  total (32/31)

GGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCG  >  NC_000913/699720‑700006                                                                                                                                               |                                                                                                                                                 ggTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGAt                                                                                                                                           <  2:631902/150‑1 (MQ=255)    tATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATaaa                                                                                                                                        >  1:540374/1‑150 (MQ=255)     aTCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATaaaaa                                                                                                                                      <  1:214172/151‑1 (MQ=255)       cGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAAc                                                                                                                                     <  1:662272/150‑1 (MQ=255)        gCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGCTaaaa                                                                                                                                       >  2:635954/1‑147 (MQ=255)            aCAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGcc                                                                                                                               >  1:27939/1‑151 (MQ=255)             cAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGccc                                                                                                                              <  1:665256/151‑1 (MQ=255)              aaTCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGGTAAAAACCCCGccc                                                                                                                              >  2:516212/1‑150 (MQ=255)                 cACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTgg                                                                                                                          <  1:874368/151‑1 (MQ=255)                  aCCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTggg                                                                                                                         >  1:585453/1‑151 (MQ=255)                   ccGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGgt                                                                                                                        <  1:277907/151‑1 (MQ=255)                      tttCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGgtg                                                                                                                       >  2:451403/1‑149 (MQ=255)                        tcttcCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTcc                                                                                                                    <  2:100072/150‑1 (MQ=255)                           tcCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGa                                                                                                                 <  1:708677/150‑1 (MQ=255)                             cGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATtg                                                                                                              <  2:497015/151‑1 (MQ=255)                                aTGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGGCAAAAAACGCCCCCGGGGGTGTCCgtggtggtg                                                                                                           >  1:804344/1‑143 (MQ=255)                                aTGCGCCGGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACGCGACGCGCGACAGGCGGGATAAAAACCGCGCCCGTGGGTGGCCGGATtggtg                                                                                                           <  2:696799/151‑1 (MQ=255)                                     ccTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGGTAAAAAACGCGCCCGTGGGTGTCCGGATTGGTGGCGa                                                                                                       >  2:91448/1‑150 (MQ=255)                                     ccTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATtggtgg                                                                                                          >  1:152587/1‑147 (MQ=255)                                     ccTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGGGTCCGGATTGGTGGCg                                                                                                        >  1:259760/1‑149 (MQ=255)                                        gCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATa                                                                                                     <  1:853322/149‑1 (MQ=255)                                              tGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAAACGCGCCCGTGGGGGTCCGGATTGGTGGCGATAAAACATGc                                                                                             >  2:578686/1‑151 (MQ=255)                                                tttAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCAc                                                                                           >  1:168883/1‑151 (MQ=255)                                                tttAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGGGTACGGATTGGGGGCGATAAAACGTGCa                                                                                            >  2:457758/1‑150 (MQ=255)                                                      gctgcGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTa                                                                                      <  2:513258/150‑1 (MQ=255)                                                         gcgcGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcg                                                                                   <  2:657593/150‑1 (MQ=255)                                                         gcgcGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGc                                                                                    <  1:267118/149‑1 (MQ=255)                                                          cgcgGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCg                                                                                                        <  2:356743/128‑1 (MQ=255)                                                          cgcgGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgac                                                                                 <  1:290063/151‑1 (MQ=255)                                                          cgcgGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGGGGCg                                                                                                        >  1:356769/1‑128 (MQ=255)                                                           gcgGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgac                                                                                 <  1:162808/150‑1 (MQ=255)                                                             ggATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacg                                                                                <  2:550500/149‑1 (MQ=255)                                                                tgaaCCGCGGGCTGGGCTTACCAACTTAGCACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGGGGGTGTCCCGATTGGTGGGGGTAAAACGTGCATCGGTTGcgccaaa                                                                              >  2:655213/1‑143 (MQ=255)                                                                  atCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGGCCGGGTTGGGGGCGATAAAACGTGCACCGTTAGCGACGAAATAgg                                                                         >  1:406215/1‑151 (MQ=255)                                                                  atCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCCTGGGGGTCCGGATTTGGGGCGATAAAACCGGCACCGGTAGCGACCGAATa                                                                           >  2:559419/1‑149 (MQ=255)                                                                  atCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCCCGCCCCCGGGTGGCCGGATTGGTGGCGATAAAACGTGCCCCCTTTGCGACGAAATAgg                                                                         >  1:568389/1‑151 (MQ=255)                                                                    ccACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGgc                                                                        <  2:232252/150‑1 (MQ=255)                                                                     cACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGAATTGGTGGCGATAAAACGTGCACCGCTAGCGACGAAAAAGgccg                                                                      >  2:447371/1‑149 (MQ=255)                                                                       cGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGgctgc                                                                     <  2:457971/150‑1 (MQ=255)                                                                            tgggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCCTGGGTGGCCGGATTGGGGGGGATAAAACGGGCACCGTTAGCGACGAAATAGGCTGCTTTAt                                                                >  2:662569/1‑150 (MQ=255)                                                                               gctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat                                                            <  2:135326/151‑1 (MQ=255)                                                                                ctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat                                                            >  1:148887/1‑150 (MQ=255)                                                                                ctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGGGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat                                                            >  1:771173/1‑150 (MQ=255)                                                                                 ttACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGGGTCCGGCTTGGGGGGGATAAAACGTGCACCGTTAGCGACGAAATAAGCTGCTTTATGcat                                                            >  2:529203/1‑149 (MQ=255)                                                                                  tACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCCTGGGGGTCCGGATTGGTGGCGATAAAACATGCACCGGTAGCGACGAAAAAGGCTGCTTTATGcatcat                                                         >  2:152987/1‑151 (MQ=255)                                                                                  tACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGAAAAAAACCGCGCCCGTTGGTGTCCGGATTAGGGGCGATAAAACGATCACCGTTAACGACGAAAAATGCTGGTTTATgga                                                             >  2:613011/1‑145 (MQ=255)                                                                                   aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat                                                         <  2:503957/150‑1 (MQ=255)                                                                                   aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat                                                         <  1:570720/150‑1 (MQ=255)                                                                                            gAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTg                                               <  2:795589/151‑1 (MQ=255)                                                                                                  tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGGTAAAAAACGCGCCCGTGGGTGTCCGGATTGGTGGGGGTAAAACGTGCACCGGTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCCGTTGTa                                             >  2:471061/1‑147 (MQ=255)                                                                                                    tGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAAc                                        <  2:520702/150‑1 (MQ=255)                                                                                                         ccGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGttt                                  >  2:439456/1‑151 (MQ=255)                                                                                                              gaTTTTCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCca                             <  1:356805/151‑1 (MQ=255)                                                                                                                ttttCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCcaa                            >  1:796214/1‑150 (MQ=255)                                                                                                                 tttCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCcaac                           <  1:158018/150‑1 (MQ=255)                                                                                                                 tttCTCAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCCCGCCCGTGGGGGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAAAAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGGACGCGTTTCGCcaa                            >  1:617131/1‑149 (MQ=255)                                                                                                                   tctcAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCcaacaa                         <  2:609513/150‑1 (MQ=255)                                                                                                                   tctcAATCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGGGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCcaacaa                         >  1:364123/1‑150 (MQ=255)                                                                                                                       aaTCCCTGCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACg                    <  1:423488/151‑1 (MQ=255)                                                                                                                              gCACACAACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCa              <  1:526997/150‑1 (MQ=255)                                                                                                                                  acaACGCGCCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGGTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAggggt         >  1:274821/1‑149 (MQ=255)                                                                                                                                        cgcCACAGGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATc   <  1:783495/151‑1 (MQ=255)                                                                                                                                          ccACAGGCGGGATAAAAACCGCGCCCGTGGGTGTGCTGGTTTGGGGGGCTAAAACGTGCACCGTTAGCGCCGAAATAGGCTGCTTTATGCATCATGTCCCAGCTGTAGGAACACGTTTCGCCAACAATCACGAAATCAGGGTTCACATCg  >  2:419380/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                 GGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCG  >  NC_000913/699720‑700006

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑