Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,840,197	T→G	L244R (CTT→CGT)	ynjE →	molybdopterin synthase sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,840,197	0	T	G	100.0%	286.1 / NA	73	L244R (CTT→CGT)	ynjE	molybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/39);  total (34/39)

AACGAAGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCTTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGT  >  NC_000913/1840049‑1840340                                                                                                                                                     |                                                                                                                                                aaCGAAGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGAttatgcgtt                                                                                                                                               <  2:642456/151‑1 (MQ=255) aaCGAAGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGAttatgcgtt                                                                                                                                               <  2:92451/151‑1 (MQ=255) aaCGAAGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGAttatgcgtt                                                                                                                                               >  1:833622/1‑151 (MQ=255)      aGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGAttatgcgttatgct                                                                                                                                          >  2:863147/1‑151 (MQ=255)      aGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGAttatgcgttatgc                                                                                                                                           >  1:223460/1‑150 (MQ=255)         ggAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGc                                                                                                                                       >  2:1286429/1‑151 (MQ=255)          gAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCg                                                                                                                                      <  2:580775/151‑1 (MQ=255)            aaGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTg                                                                                                                                    >  2:537842/1‑151 (MQ=255)                gAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGaa                                                                                                                                  >  2:139511/1‑149 (MQ=255)                  aCCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAGGAt                                                                                                                              <  2:1036709/151‑1 (MQ=255)                  aCCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGAt                                                                                                                              <  2:240866/151‑1 (MQ=255)                  aCCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGCATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGAt                                                                                                                              <  1:403492/151‑1 (MQ=255)                   ccGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGAtg                                                                                                                             >  1:1169290/1‑151 (MQ=255)                        gtgGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTgcgc                                                                                                                        <  1:99249/151‑1 (MQ=255)                          ggAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCct                                                                                                                      <  1:526972/151‑1 (MQ=255)                           gAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCctg                                                                                                                     <  2:1237741/151‑1 (MQ=255)                             aCAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCctgct                                                                                                                   >  2:1094036/1‑151 (MQ=255)                              cAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCctgctg                                                                                                                  >  1:1012776/1‑151 (MQ=255)                                aaGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGa                                                                                                                <  2:1161742/151‑1 (MQ=255)                                 aGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGAt                                                                                                               <  2:488563/151‑1 (MQ=255)                                  gTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTgg                                                                                                                 >  1:686013/1‑148 (MQ=255)                                   tttCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATgg                                                                                                             >  2:764154/1‑151 (MQ=255)                                        gatgaACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATggcg                                                                                                           <  1:440375/148‑1 (MQ=255)                                         atgaACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTg                                                                                                       >  2:1200846/1‑151 (MQ=255)                                               aacTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAAc                                                                                                 >  1:318064/1‑151 (MQ=255)                                                acTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAAcc                                                                                                >  2:411892/1‑151 (MQ=255)                                                acTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAAcc                                                                                                <  1:869267/151‑1 (MQ=255)                                                  tAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTg                                                                                              >  1:825941/1‑151 (MQ=255)                                                   aaaaGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTgg                                                                                             >  2:535875/1‑151 (MQ=255)                                                      aGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTcc                                                                                          <  1:1183234/151‑1 (MQ=255)                                                      aGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTcc                                                                                          >  1:237528/1‑151 (MQ=255)                                                      aGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGCATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTACGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTCGTcc                                                                                          >  2:669891/1‑151 (MQ=255)                                                       gCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCg                                                                                         >  2:88562/1‑151 (MQ=255)                                                        cGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCg                                                                                         >  2:199880/1‑150 (MQ=255)                                                        cGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCg                                                                                         >  2:100801/1‑150 (MQ=255)                                                        cGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCg                                                                                         >  1:552120/1‑150 (MQ=255)                                                        cGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCg                                                                                         >  2:432404/1‑150 (MQ=255)                                                         gATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGAc                                                                                       >  1:818805/1‑151 (MQ=255)                                                         gATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGAc                                                                                       <  1:759884/151‑1 (MQ=255)                                                          aTGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGAcg                                                                                      >  1:50510/1‑151 (MQ=255)                                                           tGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGAcg                                                                                      >  1:866269/1‑150 (MQ=255)                                                             cTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCgg                                                                                   >  2:1189210/1‑151 (MQ=255)                                                               ggCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGa                                                                                 <  2:621325/151‑1 (MQ=255)                                                                 cAAAACGCGGCATTCGCCATGACACCGCGGTCATTCTGTATGGGCGGGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTGATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACt                                                                               <  2:1137223/151‑1 (MQ=255)                                                                 cAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACt                                                                               <  1:262672/151‑1 (MQ=255)                                                                  aaaaCACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTg                                                                              >  2:237818/1‑151 (MQ=255)                                                                     acacGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCt                                                                           <  2:70834/151‑1 (MQ=255)                                                                     acacGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCt                                                                           <  1:912495/151‑1 (MQ=255)                                                                     acacGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTCATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCt                                                                           <  1:672063/151‑1 (MQ=255)                                                                       acGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGt                                                                         <  2:1229718/151‑1 (MQ=255)                                                                        cGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGtt                                                                        >  1:736682/1‑151 (MQ=255)                                                                        cGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGtt                                                                        <  1:259895/151‑1 (MQ=255)                                                                            aTTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAgc                                                                    <  2:935992/151‑1 (MQ=255)                                                                             ttCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATACGTTATGCTGGCGTTAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACATGGTCCGACGCGGGACTGCCTGTTGAg                                                                     <  1:351990/149‑1 (MQ=255)                                                                                gCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTTAGCGCgg                                                                >  1:808700/1‑151 (MQ=255)                                                                                   aTGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGTACTGCCTGTTGAGCGCGGAAc                                                             >  2:1150227/1‑151 (MQ=255)                                                                                   aTGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAAc                                                             <  1:825943/151‑1 (MQ=255)                                                                                      acacCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGcc                                                          <  2:218866/151‑1 (MQ=255)                                                                                       caccacGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCa                                                         <  2:678398/151‑1 (MQ=255)                                                                                             ggTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGaaa                                                   <  2:221022/151‑1 (MQ=255)                                                                                             ggTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCCGGTCCGACGTGGGACTGCCTGTTGAGCGCGGAACGCCTCCGaaa                                                   <  1:1292877/151‑1 (MQ=255)                                                                                              gTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGaa                                                    >  2:176405/1‑149 (MQ=255)                                                                                              gTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAg                                                  <  2:833655/151‑1 (MQ=255)                                                                                                 aTTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGa                                               <  2:306836/151‑1 (MQ=255)                                                                                                      gTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGaaa                                             <  2:576580/148‑1 (MQ=255)                                                                                                         tGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCgg                                         >  1:1092914/1‑149 (MQ=255)                                                                                                           ggCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAAcc                                     <  1:760912/151‑1 (MQ=255)                                                                                                                gACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGAtt                                <  2:142764/151‑1 (MQ=255)                                                                                                                   gTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGCTTTCg                             >  2:811071/1‑151 (MQ=255)                                                                                                                   gTATACGCTGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCg                             <  1:912410/151‑1 (MQ=255)                                                                                                                           tGCAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGa                       >  1:668078/1‑149 (MQ=255)                                                                                                                             cAGCGCGTGTGGCGCAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGAt                   <  2:551300/151‑1 (MQ=255)                                                                                                                                      tgGCGCAGATTATGCGTTATGATGGCGTGAAAGATGTGAGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGca            >  1:661669/1‑149 (MQ=255)                                                                                                                                         cgcAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAAcc       <  1:1301256/151‑1 (MQ=255)                                                                                                                                         cgcAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAAcc       <  1:1286350/151‑1 (MQ=255)                                                                                                                                         cgcAGATTATGAGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAAcc       >  1:61345/1‑151 (MQ=255)                                                                                                                                          gcAGATTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCg      <  2:529573/151‑1 (MQ=255)                                                                                                                                              aTTATGCGTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGt  <  2:1231369/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                AACGAAGTGGAAAGTGAACCGCTGTGGAACAAAGTTTCTGATGAACAACTAAAAGCGATGCTGGCAAAACACGGCATTCGCCATGACACCACGGTCATTCTGTATGGGCGTGACGTATACGCTGCAGCGCGTGTGGCGCAGATTATGCTTTATGCTGGCGTGAAAGATGTGCGCCTGCTGGATGGCGGCTGGCAAACCTGGTCCGACGCGGGACTGCCTGTTGAGCGCGGAACGCCACCGAAAGTGAAAGCGGAACCGGATTTCGGCGTGAAGATCCCGGCACAACCGCAGT  >  NC_000913/1840049‑1840340

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