Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,516,856	G→A	A326V (GCG→GTG)	pdeA ←	putative c‑di‑GMP phosphodiesterase PdeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,516,856	0	G	A	100.0%	239.7 / NA	72	A326V (GCG→GTG)	pdeA	putative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/41);  total (31/41)

CATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAA  >  NC_000913/2516720‑2516989                                                                                                                                         |                                                                                                                                      catcatTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCt                                                                                                                                    >  2:545392/1‑140 (MQ=255) catcatTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCt                                                                                                                                    >  1:505205/1‑140 (MQ=255) catcatTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCt                                                                                                                                    >  2:206371/1‑140 (MQ=255)          gCCGTAATGACGACTCATAAACTCAAGACTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGcc                                                                                                                           >  1:515433/1‑140 (MQ=255)            cGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAg                                                                                                                         >  2:546720/1‑140 (MQ=255)              tAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcg                                                                                                                       <  1:755422/140‑1 (MQ=255)              tAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcg                                                                                                                       >  2:112254/1‑140 (MQ=255)              tAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcg                                                                                                                       >  2:730297/1‑140 (MQ=255)              tAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcg                                                                                                                       >  1:622093/1‑140 (MQ=255)              tAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcg                                                                                                                       >  1:429152/1‑140 (MQ=255)              tAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcg                                                                                                                       >  2:625492/1‑140 (MQ=255)                aTGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGTGTTAACGGATCAGTGAGCACCGGCAAATGCCAgcggc                                                                                                                     >  1:390248/1‑140 (MQ=255)                aTGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAgcggc                                                                                                                     >  2:429456/1‑140 (MQ=255)                  gacgacTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGa                                                                                                                   <  1:206371/140‑1 (MQ=255)                  aacgacTCTTAACCTCGAGATTGTCGATGCCCATGTCGCTAAACCTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCGCGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGa                                                                                                                   <  2:398933/139‑1 (MQ=255)                    cgacTCATAAACTCAAGATTATCAATGCGCAGGCAGCCAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGAtt                                                                                                                 <  1:232838/140‑1 (MQ=255)                    cgacTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGAtt                                                                                                                 <  2:531538/140‑1 (MQ=255)                     gacTCATAAACTCAAGATTATCAATGCGTAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTa                                                                                                                >  1:631177/1‑140 (MQ=255)                      acTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGAtt                                                                                                                 <  1:528073/138‑1 (MQ=255)                      acTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGAtt                                                                                                                 >  2:528073/1‑138 (MQ=255)                      acTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTaa                                                                                                               <  1:771345/140‑1 (MQ=255)                         cATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGcc                                                                                                            <  1:730297/140‑1 (MQ=255)                          aTAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCa                                                                                                           >  2:655883/1‑140 (MQ=255)                              aCTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCa                                                                                                       >  1:70183/1‑140 (MQ=255)                                  aaGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGc                                                                                                   >  1:746559/1‑140 (MQ=255)                                   aGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCg                                                                                                  <  1:625492/140‑1 (MQ=255)                                    gATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCGGCAGATGCCAGCGGCGACTAAGCCATTCACTGCga                                                                                                 >  2:339590/1‑140 (MQ=255)                                      ttATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCgatg                                                                                               >  1:166290/1‑140 (MQ=255)                                       tATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACt                                                                                                     >  2:713456/1‑133 (MQ=255)                                       tATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACt                                                                                                     <  1:713456/133‑1 (MQ=255)                                         tCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTTGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATa                                                                                            >  2:700696/1‑140 (MQ=255)                                            aTGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAAt                                                                                         <  1:207103/140‑1 (MQ=255)                                            aTGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAAt                                                                                         >  2:550114/1‑140 (MQ=255)                                              gcgcAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgc                                                                                       <  2:665200/140‑1 (MQ=255)                                              gcgcAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgc                                                                                       >  2:53965/1‑140 (MQ=255)                                              gcgcAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgc                                                                                       <  1:112254/140‑1 (MQ=255)                                              gcgcAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgc                                                                                       <  2:622093/140‑1 (MQ=255)                                              gcgcAGGCAGCAAAAACACTTGTCCGCCTCTTGCTCCGGCGCTTGTTCAAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgc                                                                                       <  2:281826/140‑1 (MQ=255)                                               cgcAGGCACCAAAAACTCTTGCCCGCCTCTTGCGCCGGCTCTTGTTCCTACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgcg                                                                                      <  2:390248/140‑1 (MQ=255)                                                gtaggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgcgc                                                                                     <  2:631177/138‑1 (MQ=255)                                                gcaggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgcgc                                                                                     <  1:545392/140‑1 (MQ=255)                                                gcaggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgcgc                                                                                     <  1:685460/140‑1 (MQ=255)                                                gcaggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATgcgc                                                                                     <  1:804592/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTTGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  1:700696/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  1:53965/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  1:546720/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  2:703324/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  2:505205/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  2:515433/140‑1 (MQ=255)                                                 caggcagCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCa                                                                                    <  2:545636/140‑1 (MQ=255)                                                    gcagcaAAAACTCTTGCCCGCCTCTTGCTCCTGCGCTTGGTCCAACGCTCGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCa                                                                                 <  2:252983/140‑1 (MQ=255)                                                                ctTGCCCGCCTCTTGATCCGGCTCTAGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGc                                                                     <  1:339590/140‑1 (MQ=255)                                                                 ttGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCa                                                                    <  2:824185/140‑1 (MQ=255)                                                                         ctctTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTaa                                                            <  2:543440/140‑1 (MQ=255)                                                                                tCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAAt                                                     <  2:746559/140‑1 (MQ=255)                                                                                    gcgcTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGc                                                 >  1:621864/1‑140 (MQ=255)                                                                                    gcgcTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGc                                                 <  2:166290/140‑1 (MQ=255)                                                                                       cTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCAcc                                              >  1:708682/1‑140 (MQ=255)                                                                                           ttCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCaa                                          <  2:516896/140‑1 (MQ=255)                                                                                                    aCGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATg                                 >  2:689886/1‑140 (MQ=255)                                                                                                     cGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGc                                <  2:621864/140‑1 (MQ=255)                                                                                                      gAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATaaa                                   >  2:172269/1‑136 (MQ=255)                                                                                                      gAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATaaa                                   <  1:172269/136‑1 (MQ=255)                                                                                                         aGTTGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGc                            <  2:429152/140‑1 (MQ=255)                                                                                                           ttGGGTAGGAGCGTTAACGGATCTGTCAGCAGCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGa                          <  2:271396/140‑1 (MQ=255)                                                                                                            tGGGTAGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAGTCACCATCGCAAGAATAAATGCCAGCGaa                         >  1:298644/1‑140 (MQ=255)                                                                                                                 aGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTc                    <  1:655883/140‑1 (MQ=255)                                                                                                                 aGGAGCGTTAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTc                    <  2:499767/140‑1 (MQ=255)                                                                                                                  ggAGCGTTAACGGATCTGTCAGCACCTGCTAATGCCATCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCg                   <  1:220053/140‑1 (MQ=255)                                                                                                                        ttAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTGAAGGAAATGAGCCCCGGAAGAATAAGTGCCGGCGAATATTCGGTTTcc             >  2:220053/1‑140 (MQ=255)                                                                                                                        ttAACGGATCTGTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTcc             >  2:549295/1‑140 (MQ=255)                                                                                                                                   gTCAGCACCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCaa  >  2:682694/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      CATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAA  >  NC_000913/2516720‑2516989

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