Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,876,944	G→A	A122T (GCC→ACC)	iap →	aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,876,944	0	G	A	100.0%	246.0 / NA	75	A122T (GCC→ACC)	iap	aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/35);  total (40/35)

TAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTA  >  NC_000913/2876811‑2877074                                                                                                                                      |                                                                                                                                   tAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCg                                                                                                                               >  2:1353519/1‑139 (MQ=255) tAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCg                                                                                                                               <  2:259762/139‑1 (MQ=255)       tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGc                                                                                                                         >  1:1051842/1‑139 (MQ=255)              atTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTg                                                                                                                  >  2:255189/1‑139 (MQ=255)                 tataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACg                                                                                                               >  1:1026349/1‑139 (MQ=255)                  ataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                              <  2:679241/139‑1 (MQ=255)                   taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                             >  2:697803/1‑139 (MQ=255)                    aCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCg                                                                                                            >  2:986048/1‑139 (MQ=255)                     ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                        <  1:289018/126‑1 (MQ=255)                     ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                        >  2:289018/1‑126 (MQ=255)                     ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                           >  1:805739/1‑139 (MQ=255)                       gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                         >  2:158155/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGTAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  1:95022/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:278223/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:65799/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:253037/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:188724/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:682742/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:1299735/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:874838/1‑139 (MQ=255)                        cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  1:1030443/1‑139 (MQ=255)                            cgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                                    <  1:65799/139‑1 (MQ=255)                                aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                                <  1:1299735/139‑1 (MQ=255)                                aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                                <  2:661440/139‑1 (MQ=255)                                        gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                        <  2:1030443/139‑1 (MQ=255)                                        gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGACTg                                                                                        <  1:278223/139‑1 (MQ=255)                                          gCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                                                                      >  2:124847/1‑139 (MQ=255)                                           cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                     <  1:95969/139‑1 (MQ=255)                                           cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                     <  1:253037/139‑1 (MQ=255)                                           cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                     >  2:763921/1‑139 (MQ=255)                                               cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                                 >  2:547249/1‑139 (MQ=255)                                                aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                                <  1:874838/139‑1 (MQ=255)                                                aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                                <  1:986048/139‑1 (MQ=255)                                                aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                                <  1:255189/139‑1 (MQ=255)                                                   aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                                                                              <  1:225443/122‑1 (MQ=255)                                                   aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                                                                              >  2:225443/1‑122 (MQ=255)                                                    cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            >  2:772367/1‑139 (MQ=255)                                                    cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            <  2:1026349/139‑1 (MQ=255)                                                    cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            <  2:1051842/139‑1 (MQ=255)                                                        aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                                  >  2:783154/1‑129 (MQ=255)                                                        aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                                  <  1:783154/129‑1 (MQ=255)                                                        aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                                        >  2:622511/1‑139 (MQ=255)                                                            gAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGc                                                                                          <  2:1136943/117‑1 (MQ=255)                                                            gAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGc                                                                                          >  1:1136943/1‑117 (MQ=255)                                                              aGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGTCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATa                                                                  >  2:971548/1‑139 (MQ=255)                                                              aGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATa                                                                  >  2:462571/1‑139 (MQ=255)                                                                    gTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                        <  2:776223/111‑1 (MQ=255)                                                                    gTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                        >  1:776223/1‑111 (MQ=255)                                                                    gTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                            <  1:124847/139‑1 (MQ=255)                                                                    gTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                            >  2:1199638/1‑139 (MQ=255)                                                                     tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                           <  1:158155/139‑1 (MQ=255)                                                                          gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                           <  2:1437541/118‑1 (MQ=255)                                                                          gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                           >  1:1437541/1‑118 (MQ=255)                                                                          gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt                                                      >  2:604190/1‑139 (MQ=255)                                                                                cATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                                <  1:697803/139‑1 (MQ=255)                                                                                 aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtgt                                               <  1:763921/139‑1 (MQ=255)                                                                                     aGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                           <  1:622511/139‑1 (MQ=255)                                                                                                  cagcagATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAg                              <  1:772367/139‑1 (MQ=255)                                                                                                    gcagATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                            <  1:1199638/139‑1 (MQ=255)                                                                                                                  tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAc              <  1:164280/139‑1 (MQ=255)                                                                                                                  tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAc              >  2:1408874/1‑139 (MQ=255)                                                                                                                    gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                           <  2:320493/76‑1 (MQ=255)                                                                                                                    gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                           >  1:320493/1‑76 (MQ=255)                                                                                                                    gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc            >  1:1037679/1‑139 (MQ=255)                                                                                                                     cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGcc           >  2:347994/1‑139 (MQ=255)                                                                                                                       cATCTGGATACTTACACGCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         >  2:920498/1‑139 (MQ=255)                                                                                                                       cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                           >  2:1251556/1‑73 (MQ=255)                                                                                                                       cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                           <  1:1251556/73‑1 (MQ=255)                                                                                                                       cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         <  2:1320326/139‑1 (MQ=255)                                                                                                                       cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         <  2:164849/139‑1 (MQ=255)                                                                                                                       cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         >  1:164849/1‑139 (MQ=255)                                                                                                                       cATCTGGATACTTACACACCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         <  2:95022/139‑1 (MQ=255)                                                                                                                        aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa               >  1:805717/1‑132 (MQ=255)                                                                                                                        aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa               <  2:805717/132‑1 (MQ=255)                                                                                                                              aTACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTa  <  1:1408874/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   TAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTA  >  NC_000913/2876811‑2877074

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