Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,582,141	Δ1 bp	5.3%	coding (1169/1353 nt)	pnbA →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,582,141	0	T	.	5.3%	251.8 / 7.7	75	coding (1169/1353 nt)	pnbA	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base . (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CACATTGTTATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTAAAATGCCTTGGC  >  CP000730/2582006‑2582281                                                                                                                                        |                                                                                                                                              cacaTTGTTATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAG‑a                                                                                                                                        <  2:576289/141‑1 (MQ=255)   caTTGTTATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑at                                                                                                                                       <  1:340483/141‑1 (MQ=255)     ttGTTATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑atat                                                                                                                                     <  2:33875/141‑1 (MQ=255)       gTTATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAACATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAGAATATgg                                                                                                                                   <  1:423059/141‑1 (MQ=255)        ttATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGt                                                                                                                                  <  1:46412/141‑1 (MQ=255)         tATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTa                                                                                                                                 <  1:85079/141‑1 (MQ=255)          aTTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTAt                                                                                                                                >  1:508405/1‑141 (MQ=255)            ttCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTAttt                                                                                                                              <  1:82162/141‑1 (MQ=255)              cAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTg                                                                                                                            <  2:96662/141‑1 (MQ=255)                 aCAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGtt                                                                                                                         <  1:632957/141‑1 (MQ=255)                 aCAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGtt                                                                                                                         >  2:695413/1‑141 (MQ=255)                 aCAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGtt                                                                                                                         <  2:636392/141‑1 (MQ=255)                 aCAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGtt                                                                                                                         >  1:333261/1‑141 (MQ=255)                    gCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCgg                                                                                                                      >  1:97614/1‑141 (MQ=255)                    gCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCgg                                                                                                                      <  2:652933/141‑1 (MQ=255)                    gCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGAAAATTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCgg                                                                                                                      >  1:155577/1‑141 (MQ=255)                     ccGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTCGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGa                                                                                                                     >  2:594504/1‑141 (MQ=255)                             gCAATTTTTACAACAACTCAATATACAAGATTCAAACGCAAAACGACAGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACa                                                                                                             <  1:206968/141‑1 (MQ=255)                              cAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACaa                                                                                                            >  1:447482/1‑141 (MQ=255)                               aaTTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACaaa                                                                                                           <  1:295232/141‑1 (MQ=255)                               aaTTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACaaa                                                                                                           <  1:295213/141‑1 (MQ=255)                                aTTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAat                                                                                                          >  1:33614/1‑141 (MQ=255)                                 tttttACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAata                                                                                                         >  1:674395/1‑141 (MQ=255)                                  ttttACAACAACTCAATATACAAGATTCAAGCGCACATCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGGACACTAGCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACGAatag                                                                                                        <  2:137952/141‑2 (MQ=255)                                    ttACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATAGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGAAAACTTACAAATATTa                                                                                                      >  1:438344/1‑141 (MQ=255)                                     tACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                                                                     <  1:405275/141‑1 (MQ=255)                                        aacaacTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTagcag                                                                                                  <  2:566841/141‑1 (MQ=255)                                           aacTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcac                                                                                               <  2:101223/141‑1 (MQ=255)                                           aacTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcac                                                                                               >  2:427491/1‑141 (MQ=255)                                              tCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACGATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATc                                                                                            >  1:605535/1‑141 (MQ=255)                                              tCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATc                                                                                            <  2:102437/141‑1 (MQ=255)                                               cAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCa                                                                                           >  2:758604/1‑141 (MQ=255)                                                aaTATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCaa                                                                                          <  1:494750/141‑1 (MQ=255)                                                 atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAat                                                                                         >  1:435518/1‑141 (MQ=255)                                                  tataCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACATCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAata                                                                                        <  2:625393/141‑1 (MQ=255)                                                      cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                                                                    <  2:500953/141‑1 (MQ=255)                                                            tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACt                                                                              >  2:571659/1‑141 (MQ=255)                                                              aaaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGGATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGc                                                                            >  2:625094/1‑141 (MQ=255)                                                               aaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcattattttaa                                                                           >  2:487866/1‑132 (MQ=255)                                                                aCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATc                                                                                      <  2:669658/129‑1 (MQ=255)                                                                aCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATc                                                                                      >  1:669658/1‑129 (MQ=255)                                                                  gCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAc                                                                        >  1:112831/1‑141 (MQ=255)                                                                    acaACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                                                                      <  1:254003/141‑1 (MQ=255)                                                                      aaCTATGGCTTGCTGCATCTGCATGGCATGATACTTCTAGCGCACTCTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAACTTTAGCAGCACATCAATCTCCAACAACTGCTCACCTTa                                                                    <  2:213978/141‑1 (MQ=255)                                                                        cTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTaaa                                                                  <  2:318757/141‑1 (MQ=255)                                                                         tATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAt                                                                 <  2:330497/141‑1 (MQ=255)                                                                          aTGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtt                                                                <  1:54289/141‑1 (MQ=255)                                                                           tGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAttt                                                               >  1:602717/1‑141 (MQ=255)                                                                           tGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAttt                                                               >  1:696057/1‑141 (MQ=255)                                                                                tGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAg                                                          <  2:377567/141‑1 (MQ=255)                                                                                        ttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACaaa                                                  >  2:605146/1‑141 (MQ=255)                                                                                        ttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACaaa                                                  >  2:56195/1‑141 (MQ=255)                                                                                          gcatggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATg                                                <  1:147687/141‑1 (MQ=255)                                                                                           catggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGc                                               <  2:739462/141‑1 (MQ=255)                                                                                           catggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGc                                               <  2:711241/141‑1 (MQ=255)                                                                                            atggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCa                                              >  1:390187/1‑141 (MQ=255)                                                                                             tggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCaa                                             >  1:522878/1‑141 (MQ=255)                                                                                              ggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCaaa                                            <  2:41873/141‑1 (MQ=255)                                                                                                    aTACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGAtt                                      >  2:776335/1‑141 (MQ=255)                                                                                                      aCTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTa                                    >  1:287329/1‑141 (MQ=255)                                                                                                      aCTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTa                                    >  2:450730/1‑141 (MQ=255)                                                                                                            aGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAAt                              >  2:323557/1‑141 (MQ=255)                                                                                                             gcgcACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAAtt                             >  1:197467/1‑141 (MQ=255)                                                                                                                cacaCTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTg                          <  2:179611/141‑1 (MQ=255)                                                                                                                 acacTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGc                         <  1:303427/141‑1 (MQ=255)                                                                                                                     tccgaaccaaaatACCATA‑TCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTaaa                     >  1:628371/13‑141 (MQ=255)                                                                                                                      aTCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTaaaa                    >  1:309667/1‑141 (MQ=255)                                                                                                                       tCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAg                   >  2:533359/1‑141 (MQ=255)                                                                                                                         gCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCg                 <  1:336933/141‑1 (MQ=255)                                                                                                                           aGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGt               <  2:679753/141‑1 (MQ=255)                                                                                                                            gTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTa              >  2:674916/1‑141 (MQ=255)                                                                                                                            gTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTa              >  1:199684/1‑141 (MQ=255)                                                                                                                               aatACCATA‑TCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTaaaa           >  2:451652/3‑141 (MQ=255)                                                                                                                                   ccATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTAAAATGcc       <  1:59717/141‑1 (MQ=255)                                                                                                                                   ccATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTAAAATGcc       <  2:50759/141‑1 (MQ=255)                                                                                                                                       atTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTAAAATGCCTTgg   <  2:393363/141‑1 (MQ=255)                                                                                                                                        ttCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATGTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTAAAATGCCTTGGc  >  2:613229/1‑141 (MQ=255)                                                                                                                                        |                                                                                                                                              CACATTGTTATTTCAAACAGCCGATGAAGCAATTTTTACAACAACTCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAG‑ATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATGCAAAATGATTTAGCTAATTTTGCTAAAAGCGGTAAAATGCCTTGGC  >  CP000730/2582006‑2582281

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