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%	259.4 / NA	79	A122T (GCC→ACC)	iap	aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/36);  total (43/36)

GGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACC  >  NC_000913/2876819‑2877069                                                                                                                              |                                                                                                                              ggTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGa                                                                                                                  >  1:1724156/1‑139 (MQ=255) ggTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGa                                                                                                                  <  1:1386530/139‑1 (MQ=255) ggTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGa                                                                                                                  >  1:648636/1‑139 (MQ=255) ggTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGa                                                                                                                  >  2:1386530/1‑139 (MQ=255)          ataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                         >  1:1914782/1‑139 (MQ=255)             ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                      >  1:1704902/1‑139 (MQ=255)             ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                      >  1:64779/1‑139 (MQ=255)             ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                      >  2:178284/1‑139 (MQ=255)             ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCCCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                      >  1:33638/1‑139 (MQ=255)                cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                   >  2:1976880/1‑139 (MQ=255)                 ccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGcc                                                                                                  >  1:424729/1‑139 (MQ=255)                  cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCg                                                                                                                          <  2:1477048/114‑1 (MQ=255)                  cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCg                                                                                                                          >  1:1477048/1‑114 (MQ=255)                        aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                           <  2:1724156/139‑1 (MQ=255)                                gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                   <  2:1704902/139‑1 (MQ=255)                                   cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                        <  2:1056210/115‑1 (MQ=255)                                   cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                        >  1:1056210/1‑115 (MQ=255)                                   cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                          >  1:1098262/1‑129 (MQ=255)                                   cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                          <  2:1098262/129‑1 (MQ=255)                                   cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                >  1:996198/1‑139 (MQ=255)                                   cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                >  1:895035/1‑139 (MQ=255)                                    tGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGt                                                                               >  1:624595/1‑139 (MQ=255)                                     ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                              >  2:210017/1‑139 (MQ=255)                                      gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                             <  2:64779/139‑1 (MQ=255)                                        aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                           <  1:178284/139‑1 (MQ=255)                                          aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                         <  2:624595/139‑1 (MQ=255)                                           aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                        >  2:966604/1‑139 (MQ=255)                                           aCGTGACAGGAAGTACGGTGATTGACGCTCATGAAGGCAAGGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                        >  2:1256139/1‑139 (MQ=255)                                             gTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                      <  1:1256139/139‑1 (MQ=255)                                                aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                    <  2:1452323/122‑1 (MQ=255)                                                aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                    <  2:248474/122‑1 (MQ=255)                                                aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                    >  1:1452323/1‑122 (MQ=255)                                                aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                    >  1:248474/1‑122 (MQ=255)                                                     aaGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                              <  2:648636/139‑1 (MQ=255)                                                       gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATaa                                                            >  2:1523199/1‑139 (MQ=255)                                                           ggTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                        <  1:966604/139‑1 (MQ=255)                                                            gTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                       <  2:686805/139‑1 (MQ=255)                                                             tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                      <  2:424729/139‑1 (MQ=255)                                                             tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                      <  1:1523199/139‑1 (MQ=255)                                                                 tGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACAcc                                                                                                                             >  2:1785109/1‑64 (MQ=255)                                                                 tGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACAcc                                                                                                                             <  1:1785109/64‑1 (MQ=255)                                                                 tGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                           >  2:1042275/1‑114 (MQ=255)                                                                 tGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                           <  1:1042275/114‑1 (MQ=255)                                                                     gCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                        >  1:289559/1‑129 (MQ=255)                                                                     gCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                        <  2:289559/129‑1 (MQ=255)                                                                     gCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAg                                              <  2:1914782/139‑1 (MQ=255)                                                                        cATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                           <  2:1649918/139‑1 (MQ=255)                                                                         aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtgt                                          <  1:210017/139‑1 (MQ=255)                                                                                   aGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                                >  1:1123911/1‑139 (MQ=255)                                                                                   aGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                                >  2:1341774/1‑139 (MQ=255)                                                                                     cgccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCGATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAtt                              >  2:623164/1‑139 (MQ=255)                                                                                        cgcAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGa                                                                                  >  1:1574669/1‑84 (MQ=255)                                                                                        cgcAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGa                                                                                  <  2:1574669/84‑1 (MQ=255)                                                                                         gcagcaGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                   >  1:1340267/1‑66 (MQ=255)                                                                                         gcagcaGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                   <  2:1340267/66‑1 (MQ=255)                                                                                         gcagcaGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAt                               <  1:1390205/134‑1 (MQ=255)                                                                                         gcagcaGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAt                               >  2:1390205/1‑134 (MQ=255)                                                                                          cagcagATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAg                         <  2:33638/139‑1 (MQ=255)                                                                                                  catcatTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCt                                                                                                              <  2:1811931/46‑1 (MQ=255)                                                                                                  catcatTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCt                                                                                                              >  1:1811931/1‑46 (MQ=255)                                                                                                     catTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAg                                              >  2:650530/1‑107 (MQ=255)                                                                                                     catTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAg                                              <  1:650530/107‑1 (MQ=255)                                                                                                          tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                    <  2:1225422/64‑1 (MQ=255)                                                                                                          tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                    >  1:1225422/1‑64 (MQ=255)                                                                                                            gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc       >  1:1632233/1‑139 (MQ=255)                                                                                                             cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                            <  1:1349015/69‑1 (MQ=255)                                                                                                             cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                            >  2:1349015/1‑69 (MQ=255)                                                                                                             cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGcc      >  2:9522/1‑139 (MQ=255)                                                                                                               cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa    >  2:972523/1‑139 (MQ=255)                                                                                                               cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa    >  1:1521232/1‑139 (MQ=255)                                                                                                               cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa    >  1:153867/1‑139 (MQ=255)                                                                                                                aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGttt                                                <  2:305762/94‑1 (MQ=255)                                                                                                                aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGttt                                                >  1:305762/1‑94 (MQ=255)                                                                                                                aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTAc   <  2:1123911/139‑1 (MQ=255)                                                                                                                 tCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTTGGCGGGCTGACGTTACAAGGa                                                                       >  2:1352106/1‑70 (MQ=255)                                                                                                                 tCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTTGGCGGGCTGACGTTACAAGGa                                                                       <  1:1352106/70‑1 (MQ=255)                                                                                                                 tCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTAcc  <  2:1061507/139‑1 (MQ=255)                                                                                                                  cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                           <  1:97822/97‑1 (MQ=255)                                                                                                                  cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                           >  2:97822/1‑97 (MQ=255)                                                                                                                              |                                                                                                                              GGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACC  >  NC_000913/2876819‑2877069

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