Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,785	C→T	*239* (TAG→TAA)	nanC ←	N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,785	0	C	T	100.0%	258.9 / NA	80	*239* (TAG→TAA)	nanC	N‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/42);  total (38/42)

GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGT  >  NC_000913/4538649‑4538914                                                                                                                                         |                                                                                                                                  gtCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                 >  1:1475767/1‑139 (MQ=255)  tCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTG‑TTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               >  1:1494352/1‑139 (MQ=255)   cGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               <  2:1719213/139‑1 (MQ=255)     ttATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGt                                                                                                                             >  2:775093/1‑139 (MQ=255)        tCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                          <  2:1475767/139‑1 (MQ=255)          aaTTGCTCCGGTAACACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                        <  1:154118/139‑1 (MQ=255)                  cGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGacacac                                                                                                                >  2:1289090/1‑139 (MQ=255)                   ggTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACAcc                                                                                                               <  2:1495870/139‑1 (MQ=255)                   ggTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTG‑TTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCa                                                                                                              <  2:1494352/139‑1 (MQ=255)                    gTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCa                                                                                                              >  2:1175435/1‑139 (MQ=255)                     tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                             <  1:775093/139‑1 (MQ=255)                     tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                             <  1:1289090/139‑1 (MQ=255)                               aaaTGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               >  2:655794/1‑111 (MQ=255)                               aaaTGGCACAGGAGTTTCCGGTAATACAGACGCCTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACa                                                                                                                               <  1:655794/111‑1 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACAcc                                                                                                               <  2:174640/121‑1 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACAcc                                                                                                               >  1:174640/1‑121 (MQ=255)                                     cacaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTc                                                                                             >  1:363256/1‑139 (MQ=255)                                       caGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                 <  1:1678023/101‑1 (MQ=255)                                                     aaTACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             <  1:881811/139‑1 (MQ=255)                                                           gACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGCCACACCACTGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTg                                                                       >  2:806383/1‑139 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                     >  2:1198467/1‑106 (MQ=255)                                                              gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTACTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                     <  1:1198467/106‑1 (MQ=255)                                                                 tttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGTGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAaca                                                                 >  1:824700/1‑139 (MQ=255)                                                                  ttGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAacac                                                                >  2:1621590/1‑139 (MQ=255)                                                                   tGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACAcc                                                               >  1:562072/1‑139 (MQ=255)                                                                     cGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACACTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCt                                                             <  1:1249495/139‑1 (MQ=255)                                                                      ggCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGTGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTg                                                            <  2:824700/139‑1 (MQ=255)                                                                       gCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGa                                                           >  2:1750792/1‑139 (MQ=255)                                                                        cAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGAc                                                          <  1:806383/139‑1 (MQ=255)                                                                          aaTGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAAt                                                                                     >  1:1325504/1‑110 (MQ=255)                                                                          aaTGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAAt                                                                                     <  2:1325504/110‑1 (MQ=255)                                                                                gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                  <  1:1621590/139‑1 (MQ=255)                                                                                       tGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           <  1:1219569/139‑1 (MQ=255)                                                                                          tAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATAtt                                        <  1:1175435/139‑1 (MQ=255)                                                                                           aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCg                                                                          <  1:1002730/104‑1 (MQ=255)                                                                                           aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCg                                                                          >  2:1002730/1‑104 (MQ=255)                                                                                           aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                              <  1:1825755/132‑1 (MQ=255)                                                                                           aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                              >  2:1825755/1‑132 (MQ=255)                                                                                             cAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCta                                     >  1:257992/1‑139 (MQ=255)                                                                                               aGCTCCGTTATTGTTTTATTCGTTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCGCTg                                                                             <  1:1179958/97‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTcc                                                                                            <  2:908264/82‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTcc                                                                                            >  1:908264/1‑82 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATACATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAATCATATTCtatg                                   >  2:572445/1‑139 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             >  2:1179958/1‑97 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                  >  2:1426123/1‑124 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                  <  1:1426123/124‑1 (MQ=255)                                                                                               aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatg                                   <  2:562072/139‑1 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                           <  2:650212/97‑1 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                           >  1:650212/1‑97 (MQ=255)                                                                                                 cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                 <  1:1750792/139‑1 (MQ=255)                                                                                                    cGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATgg                              >  1:372867/1‑139 (MQ=255)                                                                                                     gTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCaa                                                    >  1:1639219/1‑116 (MQ=255)                                                                                                     gTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCaa                                                    <  2:1639219/116‑1 (MQ=255)                                                                                                      ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           <  1:1678065/124‑1 (MQ=255)                                                                                                      ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                           >  2:1678065/1‑124 (MQ=255)                                                                                                      ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCg                            <  2:372867/139‑1 (MQ=255)                                                                                                               tATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCgg                   <  2:363256/139‑1 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             >  1:861477/1‑79 (MQ=255)                                                                                                                 ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                             <  2:861477/79‑1 (MQ=255)                                                                                                                    aTTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                     >  1:522592/1‑52 (MQ=255)                                                                                                                    aTTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                     <  2:522592/52‑1 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                <  2:621124/135‑1 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                >  1:621124/1‑135 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCat            >  1:818204/1‑139 (MQ=255)                                                                                                                      tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCat            >  2:1289703/1‑139 (MQ=255)                                                                                                                        tgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata           >  2:693770/1‑138 (MQ=255)                                                                                                                        tgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata           <  1:693770/138‑1 (MQ=255)                                                                                                                         gACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATg         <  2:257992/139‑1 (MQ=255)                                                                                                                           cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                <  2:331470/130‑1 (MQ=255)                                                                                                                           cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                >  1:331470/1‑130 (MQ=255)                                                                                                                            tGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTa      <  1:1817079/139‑1 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCa             >  2:847213/1‑130 (MQ=255)                                                                                                                              tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCa             <  1:847213/130‑1 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACAGCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATgg        >  1:1315623/1‑133 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGt                                                       <  1:1248106/86‑1 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATgg        <  2:1315623/133‑1 (MQ=255)                                                                                                                                tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGt  >  1:1283778/1‑139 (MQ=255)                                                                                                                                tcCGGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGt                                                       >  2:1248106/1‑86 (MQ=255)                                                                                                                                   tGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTaaa                                                                   <  2:1794969/71‑1 (MQ=255)                                                                                                                                   tGTCTATTACAGTTTAAATGACAAACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTaaa                                                                   >  1:1794969/1‑71 (MQ=255)                                                                                                                                     tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCt                                      >  2:908172/1‑98 (MQ=255)                                                                                                                                     tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCt                                      <  1:908172/98‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GTCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGT  >  NC_000913/4538649‑4538914

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