Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,565,966	C→T	*427* (TAG→TAA)	yjiN ←	zinc‑type alcohol dehydrogenase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,565,966	0	C	T	100.0%	286.6 / NA	87	*427* (TAG→TAA)	yjiN	zinc‑type alcohol dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/40);  total (47/40)

AGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCCG  >  NC_000913/4565843‑4566100                                                                                                                            |                                                                                                                                       aGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         <  1:1039254/139‑1 (MQ=255)                cttccCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCg                                                                                                                            >  2:1364646/1‑121 (MQ=255)                cttccCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCg                                                                                                                            <  1:1364646/121‑1 (MQ=255)                             cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAAc                                                                                                                   <  1:1543874/117‑1 (MQ=255)                             cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAAc                                                                                                                   >  2:1543874/1‑117 (MQ=255)                               ggAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                                  >  2:20012/1‑116 (MQ=255)                               ggAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACa                                                                                                                  <  1:20012/116‑1 (MQ=255)                                  gaTAATGAAAACAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                        <  2:1230372/139‑1 (MQ=255)                                   aTAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                        >  2:222624/1‑138 (MQ=255)                                   aTAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                        <  1:222624/138‑1 (MQ=255)                                            agagGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGaa                                                                                                                    >  2:1345885/1‑101 (MQ=255)                                            agagGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGc                                                                                                       <  2:1399490/114‑1 (MQ=255)                                            agagGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGc                                                                                                       >  1:1399490/1‑114 (MQ=255)                                            agagGCCACATTATTACCATTCTTATAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGaa                                                                                                                    <  1:1345885/101‑1 (MQ=255)                                                 ccACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                        >  1:1682099/1‑124 (MQ=255)                                                 ccACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAatat                                                                                        <  2:1682099/124‑1 (MQ=255)                                                          taCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCa                                                                >  2:1320471/1‑139 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         <  2:1188141/77‑1 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         >  1:1188141/1‑77 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACaa                                                                      >  2:1318076/1‑128 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACaa                                                                      <  1:1318076/128‑1 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                                 <  2:358295/133‑1 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGcc                                                                 >  1:358295/1‑133 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                           >  2:1837936/1‑139 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                           >  1:89247/1‑139 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                           >  2:67545/1‑139 (MQ=255)                                                               ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                           >  1:629480/1‑139 (MQ=255)                                                                    ttAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTAc                                                      >  2:1607594/1‑139 (MQ=255)                                                                          tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCg                                                                                                                            <  1:1290817/63‑1 (MQ=255)                                                                          tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCg                                                                                                                            >  2:1290817/1‑63 (MQ=255)                                                                          tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCg                                                    >  2:727380/1‑135 (MQ=255)                                                                          tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCg                                                    <  1:727380/135‑1 (MQ=255)                                                                           ggTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGTTTTTTAAAAATTGCCGAGGGGGaa                                                                                                                    >  1:1440645/1‑70 (MQ=255)                                                                           ggTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGaa                                                                                                                    <  2:1440645/70‑1 (MQ=255)                                                                            gTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGaca                                              >  1:992979/1‑139 (MQ=255)                                                                             tttATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGacac                                             >  1:562033/1‑139 (MQ=255)                                                                              ttATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGcc                                                                              <  2:1224999/105‑1 (MQ=255)                                                                              ttATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGcc                                                                              >  1:1224999/1‑105 (MQ=255)                                                                               tATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGa                                           >  2:1552089/1‑139 (MQ=255)                                                                                 tttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAA‑CGCCAACCAGCGTACCGTTGACACGAATa                                        >  2:1692290/1‑139 (MQ=255)                                                                                     aTCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        <  1:292746/136‑1 (MQ=255)                                                                                     aTCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        >  2:292746/1‑136 (MQ=255)                                                                                      tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                         <  2:1528272/118‑1 (MQ=255)                                                                                      tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCg                                                         >  1:1528272/1‑118 (MQ=255)                                                                                         ctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        <  2:1083060/132‑1 (MQ=255)                                                                                         ctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        >  1:1083060/1‑132 (MQ=255)                                                                                           ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         >  1:627240/1‑49 (MQ=255)                                                                                           ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         <  2:627240/49‑1 (MQ=255)                                                                                           ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGa                               >  1:972775/1‑139 (MQ=255)                                                                                            tcGCGGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGAt                              >  2:1235688/1‑139 (MQ=255)                                                                                              gcgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAA‑CGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCtt                           <  1:1692290/139‑1 (MQ=255)                                                                                               cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                                   <  2:415243/67‑1 (MQ=255)                                                                                               cgGTTAGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCg                                                                                                   >  1:415243/1‑67 (MQ=255)                                                                                                 gTTAGCAACTTAGTTTTTAATGGATTTTTAAAAATGGCCGAGGGGAAACAAGGCCGGGAGCTGCGACGGCAAATATAAAATTAGCCCAATACACCCGCCAACCAGCGTACCGTTAACACGAATAAACTGCAGAT‑TTTGc                        <  1:1235688/139‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         >  1:245937/1‑40 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAgg                                                                                                                         <  2:245937/40‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAACTTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                      >  1:1343811/1‑139 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                   >  2:1297266/1‑78 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                   <  1:1297266/78‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGcc                                                                              >  1:440911/1‑83 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGcc                                                                              <  2:440911/83‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        >  1:209059/1‑121 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        <  1:714291/121‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        >  2:714291/1‑121 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        <  2:209059/121‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                      <  2:89247/139‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                      >  2:551668/1‑139 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                      <  2:562033/139‑1 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                      >  2:693646/1‑139 (MQ=255)                                                                                                    aGCCACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                      <  1:551668/139‑1 (MQ=255)                                                                                                      ccACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACCAGGCCGGGTGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCa                                                                >  2:257874/1‑95 (MQ=255)                                                                                                      ccACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCa                                                                <  1:257874/95‑1 (MQ=255)                                                                                                       cACTTAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGata                   <  2:1040636/139‑1 (MQ=255)                                                                                                          ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCa                                                                                             <  1:751045/62‑1 (MQ=255)                                                                                                          ttAGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCa                                                                                             >  2:751045/1‑62 (MQ=255)                                                                                                            aGTTTTTCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTaa              <  1:1320471/139‑1 (MQ=255)                                                                                                              tttttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        >  1:1165072/1‑111 (MQ=255)                                                                                                              tttttCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                        <  2:1165072/111‑1 (MQ=255)                                                                                                                  tCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                              <  1:496368/53‑1 (MQ=255)                                                                                                                  tCATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGc                                                                                              >  2:496368/1‑53 (MQ=255)                                                                                                                   cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAg                                                                                               >  2:201500/1‑51 (MQ=255)                                                                                                                   cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAg                                                                                               <  1:201500/51‑1 (MQ=255)                                                                                                                   cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                                >  1:1597961/1‑66 (MQ=255)                                                                                                                   cATGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAg                                                                                <  2:1597961/66‑1 (MQ=255)                                                                                                                     tGGATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTg     <  2:629480/139‑1 (MQ=255)                                                                                                                      ggATTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGc    >  1:1803594/1‑139 (MQ=255)                                                                                                                        aTTTTTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCcg  >  1:631304/1‑139 (MQ=255)                                                                                                                            |                                                                                                                                       AGTGACATGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCCG  >  NC_000913/4565843‑4566100

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