Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,315	G→A	12.1%	Y79Y (TAC→TAT)	tpn ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,315	0	G	A	12.1%	152.1 / 14.3	73	Y79Y (TAC→TAT)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base A (3/6);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

AAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36180‑36455                                                                                                                                        |                                                                                                                                             aaaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGt                                                                                                                                         >  2:528858/1‑141 (MQ=32)  aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGGATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCCCCATTTTCCTTTTATTTTGATGTACGtc                                                                                                                                        >  1:496123/1‑141 (MQ=17)  aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtc                                                                                                                                        <  2:391854/141‑1 (MQ=32)  aaTTGGTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtc                                                                                                                                        <  2:60829/141‑1 (MQ=17)   aTTGTTTAATGAGACGATTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtcc                                                                                                                                       >  1:304844/1‑140 (MQ=17)         ttaTGAGAAGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCaa                                                                                                                                 <  1:28464/139‑1 (MQ=255)         tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                                                                 >  1:232590/1‑141 (MQ=32)          aaTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                                                                >  1:379554/1‑141 (MQ=32)           aTGAGACGTTTTATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                                                               <  2:24372/141‑1 (MQ=17)           aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGGCTCATCACTa                                                                                                                               >  1:494296/1‑141 (MQ=17)              agaCGTTTAATAAACGGATATGCTGAATGATTATCTCGTTGCCTACGCAACCAAAGATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  1:232474/141‑1 (MQ=11)                aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTCCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGTAATCAATGAAAGAGCTACATTTTTCCTGAATTTTGATATACGTCTCATCAATACGCCa                                                                                                                          >  1:536405/1‑141 (MQ=11)                   tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCTTCAATACGCCAttt                                                                                                                       >  1:70300/1‑141 (MQ=17)                    ttGATAAACGCATAAGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCAtttt                                                                                                                      >  1:336644/1‑140 (MQ=17)                      taTAAAGCCATATTCTTGATGTTTATGCCGTTCCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGta                                                                                                                    <  2:541692/140‑1 (MQ=11)                      gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtg                                                                                                                    >  2:113425/1‑140 (MQ=17)                                 atGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTCCGTTTCAACAATACGCCATTTGTAATAAGCttttt                                                                                                         <  1:145630/141‑1 (MQ=255)                                   gCTGAATGATTATCTCGTTGCTTACGCAATCCATTATCTTATGTATGACCCTCTGCATCTATGGCACCATATAAATACCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTa                                                                                                       <  2:205022/141‑1 (MQ=11)                                      gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  1:421508/141‑1 (MQ=40)                                      gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:302756/141‑1 (MQ=38)                                      gAATGATTATCGTGTTGCTTACGCAACCAAATATCTTATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:517391/141‑1 (MQ=255)                                          gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTGGTACGAAGCTTTTTTTTGCtttt                                                                                                >  1:464354/1‑141 (MQ=255)                                                     tGCTTACGCAACCAAATATCTAATGTATGTCCCTCTACCTCAATGGAACGATATAAATAGCTCCATTTTCCTTTTCTTTTGATGAACGTTTAATCACTACGCCATTTGTAATAGAATTTTTTGTGTTTTTTGTTCTaacgt                                                                                     >  1:269360/1‑138 (MQ=11)                                                        ttACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGa                                                                                  >  2:576262/1‑141 (MQ=40)                                                        ttACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTAAGTCTTATGAAGACGCCTTTTGTAATCAGCTTTTTTATGGTTTTTCTTCCAAATTTGa                                                                                  >  1:1681/1‑141 (MQ=11)                                                        gtACGCAACCAAAGATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGa                                                                                  <  1:171576/140‑1 (MQ=255)                                                         tACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGATTTTTTCTGCGTTTTCTTCCCACTTTGAt                                                                                 >  2:426853/1‑141 (MQ=14)                                                          aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  1:431712/1‑141 (MQ=255)                                                          aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  1:527987/1‑141 (MQ=40)                                                           cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAc                                                                               <  1:312544/141‑1 (MQ=40)                                                           agCAACCAAATATCTAATGTATGTCTCTCTGCATCAATTGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAc                                                                               <  1:355394/140‑1 (MQ=255)                                                             cAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACaa                                                                             <  2:246764/141‑1 (MQ=255)                                                             cAACCAAATATCTAATGTATGGCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACaa                                                                             <  2:435338/141‑1 (MQ=255)                                                               aCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAATCTTTTTTATGGTTTTTCTTCCAAATTTGATACaaaa                                                                           >  1:540058/1‑141 (MQ=255)                                                               aCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAACTTTTTTATGCTTTTTCCTCCAAATTTGATATaaaa                                                                           >  1:243191/1‑141 (MQ=255)                                                                 cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAAtt                                                                         >  2:74738/1‑141 (MQ=255)                                                                     tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTgggg                                                                     <  2:232590/141‑1 (MQ=255)                                                                     tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGATCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCTTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATACAAAATTgggg                                                                     >  2:336857/1‑141 (MQ=255)                                                                      atCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGc                                                                    >  2:174817/1‑141 (MQ=255)                                                                              tatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCtt                                                            >  2:486745/1‑141 (MQ=255)                                                                               ttgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                           <  1:61055/140‑1 (MQ=255)                                                                                 gtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGac                                                         >  2:332576/1‑140 (MQ=255)                                                                                  tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAc                                                        >  2:182778/1‑141 (MQ=21)                                                                                  tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAACTGGGGGCATCTTGtttcac                                                        >  1:380277/1‑137 (MQ=14)                                                                                  tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGTGGCATATTCTTGAAc                                                        >  1:142263/1‑141 (MQ=255)                                                                                  tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAc                                                        <  1:123951/141‑1 (MQ=40)                                                                                    ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAccc                                                      <  1:377296/141‑1 (MQ=40)                                                                                          cATCAATGGTACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATATGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCTTATCCTTGAACCCAACGGt                                                <  1:455784/141‑1 (MQ=255)                                                                                              aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACGAAATTGGGGCATATCCTTGAACCCaaagtgagcc                                            >  2:483970/1‑133 (MQ=14)                                                                                                  gCACGTTATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTCTACCGtt                                        <  1:216200/141‑1 (MQ=255)                                                                                                     cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTTaa                                     >  1:68961/1‑141 (MQ=255)                                                                                                     cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGaa                                     >  1:96042/1‑141 (MQ=40)                                                                                                       atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:379554/141‑1 (MQ=40)                                                                                                           aattAGCTCCATTTTCCTCTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAAtgatga                               <  1:483970/138‑1 (MQ=255)                                                                                                            aaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGATATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGaa                              >  1:541872/1‑141 (MQ=11)                                                                                                            aaTACCGCCATTTTCCTTTTATTTTGATGTACGTCGCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATCAGAAAAAACACTGGGGCATATTCTTGAACCCAACGGTAGACCGTTCAATGATGaa                              <  2:585779/141‑1 (MQ=11)                                                                                                               aGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCACCGGTAGACCGTTGGATGATGAACGt                           >  1:544553/1‑141 (MQ=255)                                                                                                               aGCTCCATTTTCCTGTGATTTTGATGTACGCCTCATCAATACGCCATTTGTAATAAGCGTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATGGAGGAACATTCTTCAACCCAACCGACAACCCTGGGATGATGAACGt                           >  2:559397/1‑141 (MQ=11)                                                                                                                 cTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGttt                         >  1:392754/1‑141 (MQ=255)                                                                                                                    cATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca                      >  1:481481/1‑141 (MQ=40)                                                                                                                     aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATACAACATTGGGGCATATTCTTGAACCCAACGGTAGACCGCTGCATGGTCACCGtttttac                     >  2:247622/1‑137 (MQ=14)                                                                                                                     aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                     >  1:525587/1‑141 (MQ=40)                                                                                                                      ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                    <  1:38464/141‑1 (MQ=40)                                                                                                                          ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGt                <  2:96042/141‑1 (MQ=40)                                                                                                                          ccTTTTATTTTGACAAACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGTACGTTTACACCACGt                <  2:467641/141‑1 (MQ=11)                                                                                                                            ttttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTc              <  1:455602/141‑1 (MQ=255)                                                                                                                             tttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTcc             >  1:32683/1‑141 (MQ=40)                                                                                                                              ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            <  1:174817/141‑1 (MQ=40)                                                                                                                              ttATTGTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            <  1:300952/141‑1 (MQ=255)                                                                                                                                tttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGGTGTGTCTGCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  1:17232/140‑1 (MQ=14)                                                                                                                                aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  2:461496/141‑1 (MQ=39)                                                                                                                                     gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGGAGACCGTTGAATGATGAACGTTTAAACAACGTGCCCTTGatat     >  1:507434/1‑141 (MQ=14)                                                                                                                                     gATGTACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAAATGGGAAATATAATTACACCCAACGGTAACACGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat     >  1:580722/1‑141 (MQ=14)                                                                                                                                        gTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAGATTTGATACAAAATTGGGGAATATTCGTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTACCTTAATATTTc  >  2:336583/1‑141 (MQ=255)                                                                                                                                        |                                                                                                                                             AAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36180‑36455

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