Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	117,798	A→G	L283P (CTA→CCA)	nadC ←	quinolinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	117,798	0	A	G	100.0%	269.9 / NA	74	L283P (CTA→CCA)	nadC	quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/35);  total (39/35)

CAAGAAATTTGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTAGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCATCTGTTCTGTTT  >  NC_000913/117669‑117930                                                                                                                                  |                                                                                                                                     cAAGAAATTTGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGa                                                                                                                             <  2:1294516/139‑1 (MQ=255)    gAAATTTGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACgg                                                                                                                          >  2:1466282/1‑139 (MQ=255)         ttGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     <  1:417519/139‑1 (MQ=255)         ttGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     <  1:1466282/139‑1 (MQ=255)          tGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATaa                                                                                                                    >  2:886841/1‑139 (MQ=255)          tGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATaa                                                                                                                    >  2:589696/1‑139 (MQ=255)                gCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCa                                                                                                              >  2:770333/1‑139 (MQ=255)                 cGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCAc                                                                                                             >  2:1495795/1‑139 (MQ=255)                  gTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACg                                                                                                            <  1:510371/139‑1 (MQ=255)                       gaagaaAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGt                                                                                                       <  1:770333/139‑1 (MQ=255)                           aaaGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGt                                                                                                                 >  1:1355274/1‑125 (MQ=255)                           aaaGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGt                                                                                                                 <  2:1355274/125‑1 (MQ=255)                           aaaGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCg                                                                                                   >  2:1173816/1‑139 (MQ=255)                              gTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGTGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCa                                                                                                <  1:302279/139‑1 (MQ=255)                                       ggggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       >  2:1150939/1‑139 (MQ=255)                                       ggggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:1495795/139‑1 (MQ=255)                                       ggggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:1090294/139‑1 (MQ=255)                                       ggggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:589696/139‑1 (MQ=255)                                       ggggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:565275/139‑1 (MQ=255)                                       ggggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       >  2:1090294/1‑139 (MQ=255)                                         ggAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAg                                                                                     >  1:105415/1‑139 (MQ=255)                                             aTTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACg                                                                                        >  2:328372/1‑132 (MQ=255)                                             aTTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACg                                                                                        <  1:328372/132‑1 (MQ=255)                                                  ccgcccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTAAGTTGGCGCACCGACGGAGATAAAGTCTACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTc                                                                            >  2:971665/1‑139 (MQ=255)                                                   cgcccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCa                                                                           <  2:3451/139‑1 (MQ=255)                                                   cgcccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCa                                                                           >  2:852914/1‑139 (MQ=255)                                                    gcccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAg                                                                                     >  1:345251/1‑128 (MQ=255)                                                    gcccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAg                                                                                     <  2:345251/128‑1 (MQ=255)                                                     cccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:417064/125‑1 (MQ=255)                                                     cccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       >  2:417064/1‑125 (MQ=255)                                                     cccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGt                                                                         >  1:440717/1‑139 (MQ=255)                                                     cccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGt                                                                         >  1:712555/1‑139 (MQ=255)                                                      ccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCAc                                                                                         >  2:213112/1‑122 (MQ=255)                                                      ccgcGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCAc                                                                                         <  1:213112/122‑1 (MQ=255)                                                       cgcgATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACg                                                                                                                           <  2:281337/87‑1 (MQ=255)                                                       cgcgATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACg                                                                                                                           >  1:281337/1‑87 (MQ=255)                                                         cgATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACg                                                                     >  1:628281/1‑139 (MQ=255)                                                         cgACGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAATCCACGCCCGTTTCGACAAATTCACGCAGTGTTTTGTCAGTGAct                                                                     >  1:1038391/1‑138 (MQ=255)                                                             ggAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGc                                                                 <  2:440717/139‑1 (MQ=255)                                                               aGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCa                                                               <  2:712555/139‑1 (MQ=255)                                                                gCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAg                                                              <  1:1150939/139‑1 (MQ=255)                                                                gCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAg                                                              <  2:628281/139‑1 (MQ=255)                                                                  ggATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     >  1:1442358/1‑82 (MQ=255)                                                                  ggATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     <  2:1442358/82‑1 (MQ=255)                                                                    aTAAATCTGTCCACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACGGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGacac                                                          >  1:936348/1‑139 (MQ=255)                                                                    aTAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     <  1:568331/80‑1 (MQ=255)                                                                    aTAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     >  2:568331/1‑80 (MQ=255)                                                                           tGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGgcg                                                                                                                                   <  1:1448057/59‑1 (MQ=255)                                                                           tGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGgcg                                                                                                                                   >  2:1448057/1‑59 (MQ=255)                                                                             tCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       >  2:1274322/1‑101 (MQ=255)                                                                             tCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:1274322/101‑1 (MQ=255)                                                                             tCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGtt                                                                                 >  1:1257397/1‑107 (MQ=255)                                                                             tCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGtt                                                                                 <  2:1257397/107‑1 (MQ=255)                                                                             tCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGAGGTTGCCAGACACTTCCAGTAg                                                 >  2:947572/1‑139 (MQ=255)                                                                             tCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAg                                                 >  2:1038306/1‑139 (MQ=255)                                                                                  tATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     <  1:1021841/66‑1 (MQ=255)                                                                                  tATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATa                                                                                                                     >  2:1021841/1‑66 (MQ=255)                                                                                  tATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCt                                            >  1:438446/1‑139 (MQ=255)                                                                                            aaCGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACg                                                                                                                           >  2:971129/1‑50 (MQ=255)                                                                                            aaCGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACg                                                                                                                           <  1:971129/50‑1 (MQ=255)                                                                                              cGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGcc                                                                                                          <  2:368182/65‑1 (MQ=255)                                                                                              cGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGcc                                                                                                          >  1:368182/1‑65 (MQ=255)                                                                                              cGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACg                                                                                        >  1:762333/1‑83 (MQ=255)                                                                                              cGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACg                                                                                        <  2:762333/83‑1 (MQ=255)                                                                                                          tCGCGTGCTTGTACGTGTTTAGTTGTCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCtt                    <  2:1038391/139‑1 (MQ=255)                                                                                                             aGTGCTTGTACGTATTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTcgc                 >  2:944698/1‑139 (MQ=255)                                                                                                              gTGCTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTcgcg                <  2:438446/139‑1 (MQ=255)                                                                                                                 cTTGTACGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCAt             >  2:1276060/1‑139 (MQ=255)                                                                                                                      aCGTGTTTAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCAtctgtt        <  2:1418325/139‑1 (MQ=255)                                                                                                                           tttAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       >  2:628739/1‑55 (MQ=255)                                                                                                                           tttAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  1:628739/55‑1 (MQ=255)                                                                                                                           tttAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCa              >  2:165626/1‑128 (MQ=255)                                                                                                                           tttAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCa              <  1:165626/128‑1 (MQ=255)                                                                                                                           tttAGTTGGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCATctgttctgtt   >  2:796451/1‑139 (MQ=255)                                                                                                                            ttagttggCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCATCTGTTCTGttt  <  1:1173816/139‑1 (MQ=255)                                                                                                                             tagttggCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       >  1:1251468/1‑53 (MQ=255)                                                                                                                             tagttggCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGc                                                                                       <  2:1251468/53‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     CAAGAAATTTGCGAGGCGTTACGAAGAAAGTTGGGGAAGGGGAGATTATCCGCCCGCGATGGAGCGGATAAATCTGTCAACTATTAGCGAAAACGCATTGAAAGGTCGAGTGCTTGTACGTGTTTAGTTAGCGCACCGACGGAGATAAAGTCCACGCCCGTTTCGGCAAATTCACGCAGTGTTTTGTCAGTGACGTTGCCAGACACTTCCAGTAGCGCCTTGCCGTTGGTGCGTTTGACGGCTTCGCGCATCTGTTCTGTTT  >  NC_000913/117669‑117930

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: ‑