Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	376,673	T→C	intergenic (+3/‑99)	mhpT → / → yaiL	3‑hydroxyphenylpropionic transporter/DUF2058 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	376,673	0	T	C	100.0%	213.3 / NA	60	intergenic (+3/‑99)	mhpT/yaiL	3‑hydroxyphenylpropionic transporter/DUF2058 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/30);  total (30/30)

CGTTACTGGCCGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGTGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCAAAGCAGGATTAGTC  >  NC_000913/376525‑376819                                                                                                                                                     |                                                                                                                                                   cGTTACTGGCCGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACgcgc                                                                                                                                                  >  2:389877/1‑151 (MQ=255)    tACTGGCCGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCtt                                                                                                                                               >  1:516924/1‑151 (MQ=255)     aCTGGCCGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCtt                                                                                                                                               >  1:1164963/1‑150 (MQ=255)     aCTGGCCGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCtt                                                                                                                                               <  1:941747/150‑1 (MQ=255)          ccGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAgg                                                                                                                                         >  2:185232/1‑151 (MQ=255)          ccGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAgg                                                                                                                                         <  1:792039/151‑1 (MQ=255)                    gCTGGCATTAGGCACTGGCACGTTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGa                                                                                                                                <  2:639309/150‑1 (MQ=255)                        gCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAggg                                                                                                                           >  1:216415/1‑151 (MQ=255)                        gCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAggg                                                                                                                           <  1:495500/151‑1 (MQ=255)                               gCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCAc                                                                                                                    <  1:637165/151‑1 (MQ=255)                                    ggCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCaa                                                                                                               <  1:747191/151‑1 (MQ=255)                                      cACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATa                                                                                                             <  1:932427/151‑1 (MQ=255)                                       aCGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATa                                                                                                             >  2:1205838/1‑150 (MQ=255)                                             ggCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTAc                                                                                                      <  1:967659/151‑1 (MQ=255)                                             ggCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTAc                                                                                                      <  2:593725/151‑1 (MQ=255)                                                gTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACgg                                                                                                    <  1:309765/150‑1 (MQ=255)                                                gTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACgg                                                                                                    <  2:938360/150‑1 (MQ=255)                                                 tAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCg                                                                                                  >  1:291749/1‑151 (MQ=255)                                                    tGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGttt                                                                                               >  2:149372/1‑151 (MQ=255)                                                         gCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAgg                                                                                          >  2:628594/1‑151 (MQ=255)                                                         gCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAgg                                                                                          >  2:580933/1‑151 (MQ=255)                                                                 aCCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCg                                                                                  >  2:3612/1‑151 (MQ=255)                                                                 aCCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCg                                                                                  <  1:466545/151‑1 (MQ=255)                                                                  ccGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGc                                                                                 >  1:335810/1‑151 (MQ=255)                                                                   cGGGTATTCTTGTTGCTGGGTTGGTGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGc                                                                                 <  1:482424/150‑1 (MQ=255)                                                                     ggTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTaa                                                                              <  1:509068/151‑1 (MQ=255)                                                                            ttgttgCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAattat                                                                       >  2:1187585/1‑151 (MQ=255)                                                                               ttgCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTc                                                                     >  2:615815/1‑150 (MQ=255)                                                                                 gCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAg                                                                   <  1:545626/150‑1 (MQ=255)                                                                                  cTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCg                                                                 >  1:544092/1‑151 (MQ=255)                                                                                    gggTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCggg                                                               <  2:1036576/151‑1 (MQ=255)                                                                                    gggTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCggg                                                               <  1:577093/151‑1 (MQ=255)                                                                                       ttGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGttt                                                            >  2:544620/1‑151 (MQ=255)                                                                                              tgtTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAAc                                                      >  2:393835/1‑150 (MQ=255)                                                                                                tttATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACaa                                                    >  1:529103/1‑150 (MQ=255)                                                                                                   aTTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGa                                                 <  2:825646/150‑1 (MQ=255)                                                                                                   aTTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGa                                                 <  2:462359/150‑1 (MQ=255)                                                                                                   aTTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGAt                                                <  1:1205824/151‑1 (MQ=255)                                                                                                    ttttGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATg                                               >  2:230249/1‑151 (MQ=255)                                                                                                     tttGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATg                                               <  2:449405/150‑1 (MQ=255)                                                                                                     tttGACGAGCCGGAGATCACGAATACAGCCGTGAGCCGATGCCGGCCGCGCCTTATTAAGCAGGGGCTAAAGCGACTTAACCCCAATACAACTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACaaata                                                 <  2:111898/148‑3 (MQ=255)                                                                                                        gatgaGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCaa                                           <  1:1066176/151‑1 (MQ=255)                                                                                                           gaGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCaaaa                                         <  1:100094/150‑1 (MQ=255)                                                                                                           gaGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAAc                                        <  2:506899/151‑1 (MQ=255)                                                                                                              ccGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAActt                                      >  2:110280/1‑150 (MQ=255)                                                                                                              ccGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAActt                                      >  2:723345/1‑150 (MQ=255)                                                                                                                tgagaTCACAAATACCGCCGTGCGCCGACGCCTGACGCGCCCTATTAGGCAGGAGCGAACGGGACTTCACCCCAATACACCTCCGGCGATTACAGGTATACTCGCTAAAAATTATTCAGGGGGTTTGGAAACAAAGACGGCAAAAatta                                     >  2:132071/2‑149 (MQ=255)                                                                                                                   gaTCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAActtacctt                                 >  1:1187345/1‑150 (MQ=255)                                                                                                                   gaTCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAActtacctt                                 >  1:943544/1‑150 (MQ=255)                                                                                                                     tCACGAATACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAActtaccttac                               >  2:560943/1‑150 (MQ=255)                                                                                                                           aTACAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCa                        <  1:158140/151‑1 (MQ=255)                                                                                                                              cAGCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGtt                     >  2:362068/1‑151 (MQ=255)                                                                                                                                gCCGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGc                   <  2:746616/151‑1 (MQ=255)                                                                                                                                 ccGTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCt                  <  1:87645/151‑1 (MQ=255)                                                                                                                                   gTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTc                 >  1:1043639/1‑150 (MQ=255)                                                                                                                                   gTGCGCCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTc                 >  2:668499/1‑150 (MQ=255)                                                                                                                                     gcgcCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCaaa              <  2:1164980/151‑1 (MQ=255)                                                                                                                                     gcgcCGATGCCTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCaaa              <  1:976093/151‑1 (MQ=255)                                                                                                                                              ccTGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCAAAGCAGGAtt      >  1:529057/1‑150 (MQ=255)                                                                                                                                                tGACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCAAAGCAGGATTAGt   >  2:783569/1‑151 (MQ=255)                                                                                                                                                 gACGCGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCAAAGCAGGATTAGTc  >  1:790185/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGTTACTGGCCGGGAAAATGCTGGCATTAGGCACTGGCACGGTCGGCGTAATGGCCGCTTCTGCACCGGGTATTCTTGTTGCTGGGTTGGCGGTGTTTATTTTGATGAGCCGGAGATCACGAATACAGCCGTGCGCCGATGCCTGACGTGCCTTATTAGGCAGGGGCGAAAGGGACTTCACCCCAATACACCTACGGCGTTTACAGGTATACTCGCTAAAAATTATTCAGCGGGTTTGGAAACAAAGATGGCAAAACTTACCTTACAAGAGCAGTTGCTCAAAGCAGGATTAGTC  >  NC_000913/376525‑376819

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