Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,085,461	A→G	S151P (TCC→CCC)	speA ←	biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,085,461	0	A	G	100.0%	289.4 / NA	81	S151P (TCC→CCC)	speA	biosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (45/36);  total (45/36)

CCAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAACG  >  NC_000913/3085326‑3085595                                                                                                                                        |                                                                                                                                       ccAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                     >  2:202133/1‑139 (MQ=255)          ccTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAg                                                                                                                            >  2:620984/1‑139 (MQ=255)              gtgCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGAccc                                                                                                                        >  2:936975/1‑139 (MQ=255)                          ctcGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGccc                                                                                                            >  2:1233040/1‑139 (MQ=255)                           tcGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCg                                                                                                           >  2:513450/1‑139 (MQ=255)                                aaTTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCgaatga                                                                                                      >  1:248553/1‑139 (MQ=255)                                      tGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAgg                                                                                                <  2:248553/139‑1 (MQ=255)                                          aGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACt                                                                                            <  2:650051/139‑1 (MQ=255)                                            gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          <  1:1299053/123‑1 (MQ=255)                                            gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          >  2:1299053/1‑123 (MQ=255)                                            gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          <  2:1130767/123‑1 (MQ=255)                                            gCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          >  1:1130767/1‑123 (MQ=255)                                                atatatTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                     >  2:508928/1‑92 (MQ=255)                                                atatatTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGAcc                                                                                                                                     <  1:508928/92‑1 (MQ=255)                                                     ttCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                 <  1:202133/139‑1 (MQ=255)                                                          ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGcggcgg                                                                              <  1:739693/137‑1 (MQ=255)                                                          ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGcggcgg                                                                              >  2:739693/1‑137 (MQ=255)                                                          ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            >  2:1369478/1‑139 (MQ=255)                                                          ggTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            >  2:700024/1‑139 (MQ=255)                                                              tttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGa                                                                                                                           >  2:950822/1‑88 (MQ=255)                                                              tttATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGa                                                                                                                           <  1:950822/88‑1 (MQ=255)                                                                      cGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACttt                                                                >  2:618963/1‑139 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGGCCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                              >  1:944297/1‑139 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCg                                                                                                                                    >  1:800643/1‑69 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCg                                                                                                                                    <  2:800643/69‑1 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          >  1:1377209/1‑95 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          <  2:1377209/95‑1 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                              >  2:955064/1‑139 (MQ=255)                                                                        ttGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCACGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                              >  2:960461/1‑139 (MQ=255)                                                                          gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                            >  1:151773/1‑139 (MQ=255)                                                                          gCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATc                                                            >  1:436324/1‑139 (MQ=255)                                                                              aCGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGAt                                                        <  2:289373/139‑1 (MQ=255)                                                                                gATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATaa                                                      >  2:1171068/1‑139 (MQ=255)                                                                                   gACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCGAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGGTCGGATCAACa                                                   <  1:936975/139‑1 (MQ=255)                                                                                   gACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACa                                                   >  1:409597/1‑139 (MQ=255)                                                                                        tACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            <  2:1423292/109‑1 (MQ=255)                                                                                        tACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            >  1:1423292/1‑109 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGa                                                                                                                           <  2:243021/60‑1 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGa                                                                                                                           >  1:243021/1‑60 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          <  2:1317020/77‑1 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCga                                                                                                          >  1:1317020/1‑77 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            <  2:1342654/107‑1 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            >  1:1342654/1‑107 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCgg                                                          >  2:101700/1‑125 (MQ=255)                                                                                          cGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCgg                                                          <  1:101700/125‑1 (MQ=255)                                                                                             gTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAAt                                         <  1:620984/139‑1 (MQ=255)                                                                                              tCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATc                                        <  1:1233040/139‑1 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                 <  1:340973/96‑1 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                 >  2:340973/1‑96 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            >  2:398619/1‑101 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            <  1:398619/101‑1 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                              >  1:1380208/1‑115 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGa                                                              <  2:1380208/115‑1 (MQ=255)                                                                                                aTGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGc                                      >  2:1107946/1‑139 (MQ=255)                                                                                                  gCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCg                                    >  1:325132/1‑139 (MQ=255)                                                                                                    cAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGtt                                  <  1:700024/139‑1 (MQ=255)                                                                                                     aGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTta                                 <  2:151773/139‑1 (MQ=255)                                                                                                           tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCg                           <  1:1107946/139‑1 (MQ=255)                                                                                                           tgCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCg                           <  1:513450/139‑1 (MQ=255)                                                                                                            gCCAGTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGt                          >  2:213335/1‑139 (MQ=255)                                                                                                                gTACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGa                      >  1:1091072/1‑139 (MQ=255)                                                                                                                 tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATa                                                       <  1:105008/105‑1 (MQ=255)                                                                                                                 tACTGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATa                                                       >  2:105008/1‑105 (MQ=255)                                                                                                                    tGCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTcc                  >  1:1272220/1‑139 (MQ=255)                                                                                                                     gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                 <  1:715632/75‑1 (MQ=255)                                                                                                                     gCCATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACgcgg                                                                                 >  2:715632/1‑75 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGc                                                                                                              <  2:297119/45‑1 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGc                                                                                                              >  1:297119/1‑45 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCa                                                                                      <  2:1574596/69‑1 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCa                                                                                      >  1:1574596/1‑69 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                             >  2:915349/1‑94 (MQ=255)                                                                                                                      ccATCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGAt                                                             <  1:915349/94‑1 (MQ=255)                                                                                                                        aTCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCg              <  2:325132/139‑1 (MQ=255)                                                                                                                         tCAACTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGc             >  2:737706/1‑139 (MQ=255)                                                                                                                           aaCTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCg                                       >  1:1557121/1‑111 (MQ=255)                                                                                                                           aaCTCGGCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCg                                       <  2:1557121/111‑1 (MQ=255)                                                                                                                                ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            >  2:620123/1‑69 (MQ=255)                                                                                                                                ggCTTTGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTg                                                                            <  1:620123/69‑1 (MQ=255)                                                                                                                                    ttGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAAcg  <  2:436324/139‑1 (MQ=255)                                                                                                                                     tGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                       >  1:1386935/1‑69 (MQ=255)                                                                                                                                     tGGGACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGt                                                                       <  2:1386935/69‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       CCAGATAGACCTTGTGCCCCATCTTCTCGCCAATTAATGCCAGGCGGATATATTCGCGGTCTTTATAACCGTTGCAGACGATGACGCTACGGGTCATGCCAGCATGTGCCAGTACTGCCATCAACTCGGCTTTGGAACCGGCTTCCAGACCCAGCGGTTCGCCCGAATGAATCAGGGACTCAATCACGCGGCGGTGCTGGTTAACTTTGATCGGATAAACAAGGAAGTAATCGCCGTTATAGCCGTAGGATTCCCTCGCACGTTTGAACG  >  NC_000913/3085326‑3085595

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