Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,321,045	A→C	L68* (TTA→TGA)	phnH ←	ribophosphonate triphosphate synthase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,321,045	0	A	C	100.0%	313.3 / NA	78	L68* (TTA→TGA)	phnH	ribophosphonate triphosphate synthase subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (32/46);  total (32/46)

CGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCG  >  NC_000913/4320897‑4321190                                                                                                                                                     |                                                                                                                                                  cGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCAt                                                                                                                                                 <  2:570972/151‑1 (MQ=255)        aCGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAg                                                                                                                                          <  2:65022/151‑1 (MQ=255)         cggcggTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGa                                                                                                                                         <  1:504388/151‑1 (MQ=255)          ggcggTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGaa                                                                                                                                        <  1:1057621/151‑1 (MQ=255)           gcggTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGAAGATATCGTTATTTCATGGGGTAGaaa                                                                                                                                       <  1:596773/151‑1 (MQ=255)           gcggTGCCGGTGGAAAGGGCGTTGAGCTGTTCACTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGaaa                                                                                                                                       <  1:643529/151‑1 (MQ=255)            cggTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGaa                                                                                                                                        >  2:365592/1‑149 (MQ=255)              gTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGCTGACGATATCGTTATTTCATGGGGTAGAAAg                                                                                                                                      >  2:1082097/1‑149 (MQ=255)                   ggTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACAcc                                                                                                                               >  1:516995/1‑151 (MQ=255)                    gTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCg                                                                                                                              <  2:276001/151‑1 (MQ=255)                    gTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCg                                                                                                                              <  1:972607/151‑1 (MQ=255)                      ggAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGc                                                                                                                            >  2:1167270/1‑151 (MQ=255)                        aaaGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCgt                                                                                                                          <  2:1008196/151‑1 (MQ=255)                              cGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGtt                                                                                                                    >  1:408357/1‑151 (MQ=255)                              cGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGtt                                                                                                                    >  2:825956/1‑151 (MQ=255)                               gTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGtt                                                                                                                    >  2:169400/1‑150 (MQ=255)                                 tGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGTTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTAt                                                                                                                  <  1:592053/150‑1 (MQ=255)                                     cTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCgg                                                                                                               >  1:856852/1‑149 (MQ=255)                                      tGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGcc                                                                                                             >  1:1123310/1‑150 (MQ=255)                                       gTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGcc                                                                                                             <  1:1005763/149‑1 (MQ=255)                                       gTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGcc                                                                                                             <  2:859627/149‑1 (MQ=255)                                       gTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAg                                                                                                           <  1:670565/151‑1 (MQ=255)                                         tCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCg                                                                                                         >  1:794667/1‑151 (MQ=255)                                          cGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCg                                                                                                         <  1:1180829/150‑1 (MQ=255)                                                    cGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCAc                                                                                               <  2:672021/150‑1 (MQ=255)                                                     gCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGc                                                                                             <  1:1000237/151‑1 (MQ=255)                                                            ccGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtgg                                                                                       <  2:225837/150‑1 (MQ=255)                                                             cGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtggtg                                                                                     >  2:793061/1‑151 (MQ=255)                                                              gTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtggtgg                                                                                    <  1:1065080/151‑1 (MQ=255)                                                                cACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGc                                                                                   <  1:719879/150‑1 (MQ=255)                                                                 aCCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGa                                                                                 <  1:277432/151‑1 (MQ=255)                                                                  ccGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGat                                                                                <  1:324428/151‑1 (MQ=255)                                                                   cgcgAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGata                                                                               <  2:446147/151‑1 (MQ=255)                                                                     cgAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATAtt                                                                             <  1:146368/151‑1 (MQ=255)                                                                      gAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTc                                                                            <  2:1133759/151‑1 (MQ=255)                                                                        aGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTc                                                                            >  1:713951/1‑149 (MQ=255)                                                                        aGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAg                                                                          <  2:999130/151‑1 (MQ=255)                                                                            cGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGGCGATATCGTTACTTCCTGGGGTAGAACGCCACACCGGCGTGTCGGTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCg                                                                        >  2:487946/1‑149 (MQ=255)                                                                            cGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGt                                                                      <  2:747924/151‑1 (MQ=255)                                                                              ccTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTg                                                                    <  1:531445/151‑1 (MQ=255)                                                                              ccTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTg                                                                    >  1:630533/1‑151 (MQ=255)                                                                                tGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGGGGTGGCGATATTCAGCGGTTGcc                                                                  >  1:917837/1‑151 (MQ=255)                                                                                  ttCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCa                                                                 <  1:920312/150‑1 (MQ=255)                                                                                    ccGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGcc                                                              <  1:243781/151‑1 (MQ=255)                                                                                     cGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGcc                                                              >  2:1028927/1‑150 (MQ=255)                                                                                       gctggctgACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCcgcg                                                           >  1:936416/1‑151 (MQ=255)                                                                                           gctgACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTg                                                       <  1:476314/151‑1 (MQ=255)                                                                                             tgACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGtt                                                         <  2:844424/147‑1 (MQ=255)                                                                                             tgACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTg                                                       >  2:1161338/1‑149 (MQ=255)                                                                                                  aGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGAt                                                >  2:808355/1‑151 (MQ=255)                                                                                                  aGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGAt                                                <  1:1146600/151‑1 (MQ=255)                                                                                                   gcggcgCGTTGGTATGAAAACGCAGGCTCTGGTTGCCGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGAt                                                <  2:986248/150‑1 (MQ=255)                                                                                                    cggcgCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGc                                              <  2:863976/151‑1 (MQ=255)                                                                                                    cggcgCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGc                                              >  2:64682/1‑151 (MQ=255)                                                                                                     ggcgCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGc                                              >  2:835201/1‑150 (MQ=255)                                                                                                            tGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGAc                                      >  1:711996/1‑151 (MQ=255)                                                                                                             ggTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACa                                     <  1:998200/151‑1 (MQ=255)                                                                                                              gTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACaa                                    <  1:1017105/151‑1 (MQ=255)                                                                                                                  gAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACTCCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCAc                                <  2:334568/151‑1 (MQ=255)                                                                                                                  gAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCa                                 <  1:323826/150‑1 (MQ=255)                                                                                                                   aaaaCGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCAc                                >  2:381199/1‑150 (MQ=255)                                                                                                                   aaaaCGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACg                               <  1:157773/151‑1 (MQ=255)                                                                                                                     aaCGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGcc                             <  1:607104/151‑1 (MQ=255)                                                                                                                       cGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCg                           <  2:71399/151‑1 (MQ=255)                                                                                                                        gCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCgg                          >  1:591431/1‑151 (MQ=255)                                                                                                                            gCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcg                      <  1:365608/151‑1 (MQ=255)                                                                                                                            gCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcg                      <  2:292994/151‑1 (MQ=255)                                                                                                                                tGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcgctc                   >  2:4124/1‑150 (MQ=255)                                                                                                                                 ggTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCAt                 <  1:681530/151‑1 (MQ=255)                                                                                                                                  gTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATg                >  1:129497/1‑151 (MQ=255)                                                                                                                                    tGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCCCGCTCATgg               >  2:282338/1‑150 (MQ=255)                                                                                                                                     gACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGcc             >  2:540050/1‑151 (MQ=255)                                                                                                                                        gATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCttt          <  2:278462/151‑1 (MQ=255)                                                                                                                                           atCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAAc       <  2:408340/151‑1 (MQ=255)                                                                                                                                            tCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAAc       >  1:101879/1‑150 (MQ=255)                                                                                                                                              gTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGCTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAgg    >  2:1135041/1‑151 (MQ=255)                                                                                                                                               ttATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAgg    >  1:1106839/1‑150 (MQ=255)                                                                                                                                                tATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGcg  >  2:224678/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCG  >  NC_000913/4320897‑4321190

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