Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,653,550	A→G	4.2%	N143N (AAT→AAC)	aroK ←	shikimate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,653,550	0	A	G	4.2%	292.4 / 2.1	96	N143N (AAT→AAC)	aroK	shikimate kinase
Reads supporting (aligned to +/- strand):  ref base A (49/41);  new base G (2/2);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

ACTTGCACCTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTT  >  CP000730/1653413‑1653688                                                                                                                                          |                                                                                                                                           aCTTGCACCTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATta                                                                                                                                         >  2:1190407/1‑141 (MQ=255) aCTTGCACCTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATta                                                                                                                                         >  1:408860/1‑141 (MQ=255) aCTTGCACCTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATta                                                                                                                                         <  1:2756/141‑1 (MQ=255)   ttGCACCTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATtata                                                                                                                                       >  1:995072/1‑141 (MQ=255)         cTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTaa                                                                                                                                 <  1:433388/141‑1 (MQ=255)         cTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATGTCTTGAATTCGATTTCATTATATCTTaa                                                                                                                                 >  1:832169/1‑141 (MQ=255)            cGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAAtt                                                                                                                              >  2:1033566/1‑141 (MQ=255)             cgCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTAATATTTTTGTCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTc                                                                                                                             <  1:567656/140‑1 (MQ=255)                gATAGAGTCTTAACGTTCTATTTTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCcgc                                                                                                                          <  1:1152365/141‑1 (MQ=255)                  tAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAACTTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTGTCAAATTCCTTGAATGCGCTTTCCTTATATCTTAAATTCcgcgc                                                                                                                        >  2:74172/1‑140 (MQ=255)                  tAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGa                                                                                                                        >  2:1014688/1‑141 (MQ=255)                  tAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGa                                                                                                                        <  2:675475/141‑1 (MQ=255)                  tAGAGACTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATCTAGCAATACACAATATATTTCTAAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGa                                                                                                                        >  2:485269/1‑141 (MQ=255)                    gagTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGc                                                                                                                      <  2:1006148/141‑1 (MQ=255)                     cgTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGAGCTTTCACTATATCTTAAATTGCTCGAgga                                                                                                                     >  1:146300/2‑139 (MQ=255)                     agTCTTAACGCTTTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCa                                                                                                                     <  1:361587/141‑1 (MQ=255)                     agTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCa                                                                                                                     >  1:599705/1‑141 (MQ=255)                       tCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAAt                                                                                                                   >  2:750651/1‑141 (MQ=255)                                 ctATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCAtt                                                                                                         <  1:1121253/141‑1 (MQ=255)                                  tATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCAttt                                                                                                        >  2:1132248/1‑141 (MQ=255)                                  tATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCAttt                                                                                                        >  2:868925/1‑141 (MQ=255)                                   aTTCTCTCTAATACTGATCCCTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTa                                                                                                       <  1:5463/141‑1 (MQ=255)                                    ttCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTaa                                                                                                      >  2:1021698/1‑141 (MQ=255)                                         tctAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCt                                                                                                 <  2:995072/141‑1 (MQ=255)                                         tctAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTGGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCt                                                                                                 <  1:1180611/141‑1 (MQ=255)                                          ttAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTGGAATGCTATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCtt                                                                                                <  1:519691/140‑1 (MQ=255)                                           tAATCCTGATCACTCGATTTTAATATATTTAGCATTTCATAATATTTTTCCTACAGTTATATCAAATAACTAACATAATTCTACAATGCGATTTACTTATATCTTAAATTCCGCTAGCAATACAAGTCATTTAACTGCTTg                                                                                               >  2:235867/1‑141 (MQ=255)                                             aTACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGAt                                                                                             >  1:899885/1‑141 (MQ=255)                                                cTGATCCCTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGATCATTACAAGTCATTTAACTGCTTGATTGt                                                                                          <  2:355589/141‑1 (MQ=255)                                                  taTCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTGCACACACAATACAAGTCATTTAACTGCTTTAtgttct                                                                                        >  2:1100168/2‑136 (MQ=255)                                                  gATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATCCAAGTCATTTAACTGCTTGATTGTCt                                                                                        >  2:144726/1‑141 (MQ=255)                                                    tCACTCGTTTTTATTAAATTTAGCAATTCATAATATATTTCTTAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCtta                                                                                      <  2:769841/141‑1 (MQ=255)                                                     taCTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTTAAATTTATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCTATACAAGTCATTTAACTGCTTGATTGTCttat                                                                                     <  2:556335/140‑1 (MQ=255)                                                      aCTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCttatt                                                                                    >  2:174785/1‑141 (MQ=255)                                                        tCGCTTTTATTTAATTTAGCATTTCATAATATTTGTCTTAAATTGATATCTAATGACTATCAGAGTTATTGAATTCGGTTTAATTTTAACATCTATTCCGCCTGCTATACAAGTCATTTAACTGCTTGATTCTCttagttt                                                                                  <  2:546798/141‑5 (MQ=255)                                                        tCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTtactat                                                                                  <  2:1037988/141‑1 (MQ=255)                                                         cGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTCTATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTattatt                                                                                 >  1:882722/1‑141 (MQ=255)                                                          gCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTAttt                                                                                >  2:407496/1‑141 (MQ=255)                                                                 ttaAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAg                                                                         >  1:511232/1‑141 (MQ=255)                                                                   aaaTTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGt                                                                       <  1:41036/141‑1 (MQ=255)                                                                       ttAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGAt                                                                   >  1:672660/1‑141 (MQ=255)                                                                         gacAATTCATAATATAGTTCTAAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCTCGAGCAATACAAGTCAGTTAACTGCTTGATTGTTTTATTATTTGCATTAGGTCGAtgt                                                                 <  2:585450/139‑1 (MQ=255)                                                                         aGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGAtgt                                                                 <  2:756493/141‑1 (MQ=255)                                                                          gCAGTTCATAATGTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGAATTAGGTCGAtgtg                                                                >  2:700074/1‑141 (MQ=255)                                                                          gCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGAtgtg                                                                >  1:404091/1‑141 (MQ=255)                                                                            aaTTCATAATATATTTCTGAAATTAATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGAAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTggg                                                              >  2:532152/1‑141 (MQ=255)                                                                              ttCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGACTGTCTTATTATTTGCATTAGGTCGATGTGGGtc                                                            >  2:666892/1‑141 (MQ=255)                                                                               tCATAATATATTTCTGAAATTTATAGCAAATGGCTATCATATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtca                                                           <  2:484171/141‑1 (MQ=255)                                                                                cATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCGTTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtcat                                                          >  1:46339/1‑141 (MQ=255)                                                                                cATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCGTTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtcat                                                          >  1:46428/1‑141 (MQ=255)                                                                                cATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtcat                                                          <  1:685700/141‑1 (MQ=255)                                                                                  taataTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCGTTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtcatca                                                        <  2:46339/141‑1 (MQ=255)                                                                                  taataTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCGTTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtcatca                                                        <  2:46428/141‑1 (MQ=255)                                                                                  taataTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGtcatca                                                        <  2:96731/141‑1 (MQ=255)                                                                                  acatatatTTCTGAAATTGATATCAACTGACTATCAAAGTTCTTGAATGCGATTTCAGTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTGTTATTTGCATTAGGTCGATGTGGGtcatca                                                        <  2:338588/139‑1 (MQ=255)                                                                                   aataTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTGAAGGCGTTAGGTCGATGTGGGTcatcat                                                       <  1:1187971/141‑1 (MQ=255)                                                                                   aataTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTgatcat                                                       >  2:627212/1‑141 (MQ=255)                                                                                   aataTATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATCAGGTCGATGTGGGTcatcat                                                       <  1:318053/141‑1 (MQ=255)                                                                                    atatatTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCAtt                                                      >  2:46546/1‑141 (MQ=255)                                                                                    atatatTTCTGAAATTGATAGCAAACGACAATCAAATGTCGTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGGCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCAtt                                                      <  1:512486/141‑1 (MQ=255)                                                                                    agatatTTCTGAAATTGACAGCAAATGACCACCAAAACTCCGTAATGCCATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCAtt                                                      <  1:512451/139‑1 (MQ=255)                                                                                     tacatTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTg                                                     <  1:464017/138‑1 (MQ=255)                                                                                        atTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGAtt                                                  >  2:914464/1‑141 (MQ=255)                                                                                         tttCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTc                                                 >  1:701276/1‑141 (MQ=255)                                                                                         tttCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTc                                                 <  1:841700/141‑1 (MQ=255)                                                                                         tttCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTc                                                 >  1:899848/1‑141 (MQ=255)                                                                                          ttCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCg                                                >  2:1199925/1‑141 (MQ=255)                                                                                             tGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTGGATTTTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACt                                             >  2:1147790/1‑141 (MQ=255)                                                                                              gAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACta                                            >  2:615771/1‑141 (MQ=255)                                                                                              gAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACta                                            <  2:628988/141‑1 (MQ=255)                                                                                              gAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACta                                            >  1:670449/1‑141 (MQ=255)                                                                                                aaTTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCGTTGATTCGACtatg                                          >  1:42029/1‑140 (MQ=255)                                                                                                aaTTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtata                                          >  1:197546/1‑141 (MQ=255)                                                                                                  tttaTAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtatata                                        <  2:889744/138‑1 (MQ=255)                                                                                                  ttGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTCTTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtatata                                        >  1:922780/1‑141 (MQ=255)                                                                                                  ttGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGTCtatata                                        <  2:1036430/141‑1 (MQ=255)                                                                                                  ttGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtatata                                        >  1:472538/1‑141 (MQ=255)                                                                                                    gATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtatgtata                                      >  1:549322/1‑141 (MQ=255)                                                                                                    gATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtatatata                                      >  1:995154/1‑141 (MQ=255)                                                                                                    gATAGCAAATGACTATCAAATTTCTTGAATGCGATGTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACtatatata                                      >  1:1066100/1‑141 (MQ=255)                                                                                                      tAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATataat                                    >  1:405783/1‑141 (MQ=255)                                                                                                      tAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATataat                                    >  1:1107432/1‑141 (MQ=255)                                                                                                         cAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGAGCGCGATGCAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTACATATAATATc                                 <  1:868925/141‑1 (MQ=255)                                                                                                          aaaTGACTATCAAATTTCTTTAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCa                                <  2:580631/141‑1 (MQ=255)                                                                                                          aaaTGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCa                                <  2:686867/141‑1 (MQ=255)                                                                                                          aaaTGACTATCAAATGTCTTGAATGCGATTTCATTATATCTTTAATTCCGCGAGCAATACAAATCCTTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACGATATGTAATATCa                                >  2:1117583/1‑141 (MQ=255)                                                                                                             tGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAata                             <  1:237044/141‑1 (MQ=255)                                                                                                                    cAAATTTCTTGAATTCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAAt                      <  2:863498/141‑1 (MQ=255)                                                                                                                       aTTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCCTCATTGCTTCGACTATCCATAATATCAATATTCCAAtaca                   >  2:562612/1‑138 (MQ=255)                                                                                                                              gAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAaat            >  1:54838/1‑141 (MQ=255)                                                                                                                               aaTGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAaata           <  1:1132248/141‑1 (MQ=255)                                                                                                                                aTGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTAATCATTGATTCGACTATATATAATATCAAAATTCCAATCTAACCAaataa          >  1:254629/1‑141 (MQ=255)                                                                                                                                     aTTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTAGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGttt     >  1:159253/1‑141 (MQ=255)                                                                                                                                       ttCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGttttt   <  1:508895/141‑1 (MQ=255)                                                                                                                                       ttCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGttttt   <  1:396311/141‑1 (MQ=255)                                                                                                                                        tCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGtttttt  >  2:332521/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                           ACTTGCACCTTCGGCGATAGAGTCTTAACGCTCTATTCTCTCTAATACTGATCACTCGCTTTTATTAAATTTAGCAATTCATAATATATTTCTGAAATTGATAGCAAATGACTATCAAATTTCTTGAATGCGATTTCATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTT  >  CP000730/1653413‑1653688

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