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

TCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTA  >  NC_000913/2876817‑2877074                                                                                                                                 |                                                                                                                                   tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGc                                                                                                                         <  1:1242192/139‑1 (MQ=255)              aCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCg                                                                                                            >  2:72418/1‑139 (MQ=255)                 gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGAt                                                                                                          >  2:311839/1‑139 (MQ=255)                  cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                        >  2:1342136/1‑123 (MQ=255)                  cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                        <  1:1342136/123‑1 (MQ=255)                  cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  1:737440/1‑139 (MQ=255)                  cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  2:815190/1‑139 (MQ=255)                  cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                        >  1:306030/1‑139 (MQ=255)                   ccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGcc                                                                                                       >  1:1014214/1‑139 (MQ=255)                           aTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                              >  1:451958/1‑124 (MQ=255)                           aTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                              <  2:451958/124‑1 (MQ=255)                            tCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                               <  1:311839/139‑1 (MQ=255)                            tCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                                                                              >  1:52238/1‑139 (MQ=255)                                 agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                         >  2:427837/1‑139 (MQ=255)                                  gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                        <  1:815190/139‑1 (MQ=255)                                     cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                     >  2:753618/1‑139 (MQ=255)                                       ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                                   >  2:1137422/1‑139 (MQ=255)                                        gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                                  <  1:427837/139‑1 (MQ=255)                                          aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                                <  1:72418/139‑1 (MQ=255)                                              cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAtctc                                                                                                 >  1:781792/1‑118 (MQ=255)                                              cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAtctc                                                                                                 <  2:781792/118‑1 (MQ=255)                                                tGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTg                                                                                                                  >  2:613368/1‑99 (MQ=255)                                                tGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTg                                                                                                                  <  1:613368/99‑1 (MQ=255)                                                  aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                             <  1:831708/102‑1 (MQ=255)                                                  aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                             >  2:831708/1‑102 (MQ=255)                                                  aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            >  1:24585/1‑135 (MQ=255)                                                  aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            <  2:24585/135‑1 (MQ=255)                                                  aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                                        >  2:735342/1‑139 (MQ=255)                                                       aaGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                                   <  1:735342/139‑1 (MQ=255)                                                        aGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGc                                                                                          <  1:519792/115‑1 (MQ=255)                                                        aGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGc                                                                                          >  2:519792/1‑115 (MQ=255)                                                          tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                             >  1:492256/1‑94 (MQ=255)                                                          tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                             <  2:492256/94‑1 (MQ=255)                                                                gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            >  1:1426882/1‑121 (MQ=255)                                                                gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                            <  2:1426882/121‑1 (MQ=255)                                                                 aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                               >  2:629704/1‑101 (MQ=255)                                                                 aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                               <  1:629704/101‑1 (MQ=255)                                                                  ttGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                                   <  1:203374/112‑1 (MQ=255)                                                                  ttGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                                   >  2:203374/1‑112 (MQ=255)                                                                    gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt                                                      >  2:683014/1‑139 (MQ=255)                                                                           aTGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtgt                                               <  2:737440/139‑1 (MQ=255)                                                                           aTGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtgt                                               >  2:680332/1‑139 (MQ=255)                                                                                 gCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                                   >  1:219456/1‑113 (MQ=255)                                                                                 gCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                                   <  2:219456/113‑1 (MQ=255)                                                                                     aGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAg                                                                                                                          >  1:246318/1‑54 (MQ=255)                                                                                     aGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAg                                                                                                                          <  2:246318/54‑1 (MQ=255)                                                                                     aGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                                     >  1:603720/1‑139 (MQ=255)                                                                                     aGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                                     >  2:748680/1‑139 (MQ=255)                                                                                          cgcAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                              <  1:22480/93‑1 (MQ=255)                                                                                          cgcAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                              >  2:22480/1‑93 (MQ=255)                                                                                          cgcAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                             <  2:1152271/110‑1 (MQ=255)                                                                                          cgcAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                             >  1:1152271/1‑110 (MQ=255)                                                                                           gcagcaGATCATCATTAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAGTTGGCa                               >  1:871193/1‑139 (MQ=255)                                                                                                         ttAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                            >  2:367024/1‑128 (MQ=255)                                                                                                         ttAT‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                            <  1:367024/128‑1 (MQ=255)                                                                                                            t‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                                <  2:390346/57‑1 (MQ=255)                                                                                                            t‑GGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                                >  1:390346/1‑57 (MQ=255)                                                                                                               gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGTTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc            <  2:306030/139‑1 (MQ=255)                                                                                                               gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGTTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc            <  2:1014214/139‑1 (MQ=255)                                                                                                               gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc            >  2:194089/1‑139 (MQ=255)                                                                                                               gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGACAGaa                            <  1:1087351/123‑1 (MQ=255)                                                                                                               gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGACAGaa                            >  2:1087351/1‑123 (MQ=255)                                                                                                                cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGcc           >  2:1221698/1‑139 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                              >  2:843815/1‑70 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                              <  1:843815/70‑1 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt                                                      <  2:466667/94‑1 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt                                                      >  1:466667/1‑94 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaa                   >  2:224392/1‑129 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaa                   <  1:224392/129‑1 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         <  1:680332/139‑1 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         <  1:387719/139‑1 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         >  2:387719/1‑139 (MQ=255)                                                                                                                  cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa         <  1:683014/139‑1 (MQ=255)                                                                                                                     cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                               <  2:1157163/50‑1 (MQ=255)                                                                                                                     cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                               >  1:1157163/1‑50 (MQ=255)                                                                                                                     cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTg                      <  1:372310/123‑1 (MQ=255)                                                                                                                     cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTg                      >  2:372310/1‑123 (MQ=255)                                                                                                                     cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCg      <  1:1221698/139‑1 (MQ=255)                                                                                                                         aTACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTa  <  1:1137422/139‑1 (MQ=255)                                                                                                                             ttACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                            >  1:759868/1‑77 (MQ=255)                                                                                                                             ttACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                            <  2:759868/77‑1 (MQ=255)                                                                                                                                 |                                                                                                                                   TCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTAT‑GGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTA  >  NC_000913/2876817‑2877074

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