Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,145,108	G→A	19.6%	R148K (AGA→AAA)	yegJ →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,145,108	0	G	A	19.6%	136.3 / 7.8	59	R148K (AGA→AAA)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (9/3);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCT  >  CP009273/2144964‑2145257                                                                                                                                                  |                                                                                                                                                       cAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCa                                                                                                                                                    <  1:173207/149‑1 (MQ=255)  aaaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCaaa                                                                                                                                                  <  1:507162/150‑1 (MQ=255)      gTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTa                                                                                                                                              <  2:1023818/150‑1 (MQ=255)         gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGtt                                                                                                                                           <  2:396143/150‑1 (MQ=255)              tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCt                                                                                                                                       <  1:599379/149‑1 (MQ=255)                 ccAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATc                                                                                                                                   <  1:544975/150‑1 (MQ=255)                  cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑AAACTTTATTCTTTTTTTGGCTCCTTTCCCAAAAAAAAAACCCACGAGCTAAAAGAAAAAATAGATCTCTCAt                                                                                                                                    >  1:460056/1‑148 (MQ=255)                  cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCa                                                                                                                                  <  1:307650/150‑1 (MQ=255)                   aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ATCATTTATTCTTTTTTTGTCCCCCATCCAAAAAAAGAAGACGACGAGCCAAAAGCAAAAAAAATTTGCTAAt                                                                                                                                    >  1:731654/1‑147 (MQ=255)                     gAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGcc                                                                                                                               <  1:895305/150‑1 (MQ=255)                      aGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑g                                                                                                                             <  1:287248/150‑1 (MQ=255)                          gAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTta                                                                                                                          <  2:998193/149‑1 (MQ=255)                            aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTTATCCTTTTTTTGTTACCTATCCAAAAAAAAAAAGCGGCGAGCTAAAAGGAAAAATAATTTGCTAATCAAcccgtttat                                                                                                                        >  1:624299/1‑143 (MQ=255)                                 aagaaagaAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTTATCCTTTTTTTGTTACCTATCCAAAAAAAAAAACCGACGAGCTAAGAGCAAAAATAATTTGCTTATCAACC‑GGTTtatttttt                                                                                                                  >  1:239988/1‑150 (MQ=255)                                    aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑AACCTTGATCCCTTTTTTGTTTCCTATCCCAAAAAAGAAAACGACGAGCTAAGAGCAAAAATAATTTTCTAATCAGCC‑GGTATATtttttttt                                                                                                               >  1:693074/1‑150 (MQ=255)                                         agaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑ACCATTGATTCTTTTTTTGTCCCCTATCCCAAAAAAGAAACCGCCCAGCCAAAAGCAAAAATAGTTTGCTCATCAGCC‑GCTTTATTTTTTTTTTtaaa                                                                                                          >  1:694714/1‑149 (MQ=255)                                           atatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑AACATTGATTCCTTTTTTGTGTCCCTTCCAAAAAAAGAAAACGACCAGGTAAAAGGAAAAATAGTTTTCTAAACAACC‑GTTATATTTTTTTTTTtaaat                                                                                                         >  2:835866/1‑147 (MQ=255)                                             atTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑ACCATTTATTCTTTTTTTGTCACCTATTCCAAAAAAAAAAACGACGAGCTAAAAGCAAAAATAATTCGCTAATCAACC‑GGTATATTTTTTTTTTtaaattgc                                                                                                      >  1:164570/1‑145 (MQ=255)                                              tttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑TTTATTTATTATTTTTTGGTGGCCTATCCAAAAAAAAAAGCCGACGAGGTAAGAGCAAAAATAGTTCGCTAATAAACC‑GGTATATTATTTTTTTAAAATgcgg                                                                                                     >  2:125402/1‑147 (MQ=255)                                                tCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCgg                                                                                                    <  2:357146/149‑1 (MQ=255)                                                 cAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGt                                                                                                   <  1:1037296/149‑1 (MQ=255)                                                   gACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGa                                                                                                <  1:337333/150‑1 (MQ=255)                                                    aCTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACCATT‑TTTCTTTTTTTGTTTCCCATCCCAAAAAAAAAACCGACCAGCGAAAAAGAAAAATAAATTTCTCATCAACC‑GGTATTTTTTTTTTTTTAAATGGCggggtta                                                                                               >  2:702849/1‑145 (MQ=255)                                                    aCTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑ACCACTGATTCTTTTTTTGTTACCTCTCCCAAAAAAAAAACCGACGAGCTAAAAGCAAAAATAGTTTGCTAATCAGCCGGTTTTTTTTTTTTTTTAAATTGCGGGTTGa                                                                                                >  1:774071/1‑150 (MQ=255)                                                     cTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTTATCCTTTTTTTGTTACCTATCCCAAAAAAAAAACCGACGAGCCAAAAGCAAAAATAATTTGCTAATTAACC‑GTTATATTTTTTTTTTTAAATGGCGGGTggaa                                                                                              >  1:90214/1‑147 (MQ=255)                                                        gATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAAttg                                                                                            <  2:789076/149‑1 (MQ=255)                                                        gATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAAttg                                                                                            <  2:599732/149‑1 (MQ=255)                                                            tataCGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTc                                                                                       <  2:143667/150‑1 (MQ=255)                                                                            aaaTTTACGGCGGTTAT‑ACAATTGATCTTTTATTGGTCACCTATCAAAAAGAAGAAGCCGACGAGAAAAGAGAAAAATTAGTTCGCTAACCAGCC‑GTTATTTTTATTGTTAAACCGGGGGGTTGAATTGTTGTCTTCAACCGCCATgaga                                                                       <  2:140136/150‑1 (MQ=255)                                                                                    ggcggTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTGTCTTTCACCGCCCTGAGAAaatcact                                                               >  1:1025148/1‑150 (MQ=255)                                                                                      cggTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTTTTCAACCGCCATGAGAAACTCActct                                                             >  2:1140426/1‑150 (MQ=255)                                                                                      cggTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTCAATTGTTGGTTTTATCCCCCCTGAGAAACTCActct                                                             >  2:834388/1‑150 (MQ=255)                                                                                       ggTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTg                                                            >  1:476866/1‑150 (MQ=255)                                                                                        gTTAT‑ACCATTGATCTTTTATTGGTCACATATCAAAAAGAAAAAGCCAACGAGAAAAAAGAAAAATTAGTTCGAAAACAAGCA‑GTTATATTAATTTTTAAAACGGGCGGTTGATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGc                                                           <  1:203899/150‑1 (MQ=255)                                                                                        gTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAAATCACCCTGc                                                           >  1:856586/1‑150 (MQ=255)                                                                                              aCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTTTTCAACCGCCATGAGAAAATCACTCTGCCcac                                                       >  2:839504/1‑149 (MQ=255)                                                                                                cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgc                                                    >  2:173378/1‑150 (MQ=255)                                                                                                cATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCTACCCCCCTGAGAAAATCACTCTGccccccgc                                                    >  2:346803/1‑145 (MQ=255)                                                                                                   tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTGTCCTCCACCACCATGTGAAACACACTCTGCCCACAgcga                                                  >  2:234404/1‑148 (MQ=255)                                                                                                    gATCTTTTTTTGCTCCCCTATCAAAAAAAAGAAGCCGACGAGCAAAGAAAAAAATTAGTTCTCTAATCAGCC‑TTTATTTTTATTTTTAAAACGGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc                                                 <  2:624299/149‑1 (MQ=255)                                                                                                            aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCCTTAGAAAATCACTCTGCCCCCAGCGCCCGTGCCAAc                                        >  2:1044760/1‑150 (MQ=255)                                                                                                            aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGGGCCAAc                                        >  1:971076/1‑150 (MQ=255)                                                                                                             ttgtcccccTAACAAAAAAAAGAGGCCGACGAGCTAAGAGAAAAATTATTTCGCTAATCGCCC‑TTTATTTTTATTGTTAAACCTGGCGGTTGAATTGTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACg                                       <  2:731654/143‑1 (MQ=255)                                                                                                                   aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat                                 >  1:1004128/1‑150 (MQ=255)                                                                                                                    ccTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat                                 >  1:606401/1‑149 (MQ=255)                                                                                                                    ccTATCAAAAAGAAGAAGCAGACGAGATAAGAGAAAAATTATTTCGATAATCGGCT‑TTTATATTTATTGTTAAACCGGGGGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatc                                <  2:693074/150‑1 (MQ=255)                                                                                                                              agaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAAATCCCTCTGCCCACAGCGCCAGTGCCAACGCTTTTTCATCCTCCGt                       >  1:872388/1‑149 (MQ=255)                                                                                                                              agaagaagCCGACGAGCAAAAAAAAAAATTAGTTCGAAAATCAGCT‑TTTATATTTATTGTTAAAACGGGGGGTTAAATTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTa                      <  2:774071/150‑1 (MQ=255)                                                                                                                               gaggaggCCGACGAAAAAAGAGAAAAATTATTTCACTATTCAGCC‑TTTATATTTTTTGTTTAAACGGGGGGTTGATTTGTTGTTTTCAGCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                     <  2:33163/144‑1 (MQ=255)                                                                                                                                aagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGGTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACg                    >  2:14357/1‑150 (MQ=255)                                                                                                                                 agaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAAATCACTCTGCCCACAGCGCCCGTGCCAACGGTTTATCATCATCCGTACGc                   >  1:446662/1‑150 (MQ=255)                                                                                                                                   aagCCGAGGAGCAAAAAAAAAAATTATTTCGTTAATCACCG‑TTTATTTTTTTTGTTAAAACGGGCGTTTAAATTTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTc                 <  2:881234/150‑1 (MQ=255)                                                                                                                                   aagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAACGGTTTATCATCATCCGTACGCTc                 >  1:876550/1‑150 (MQ=255)                                                                                                                                      ccGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGtt              >  2:213716/1‑150 (MQ=255)                                                                                                                                      ccGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTCTCCTCAACCGCCATGAGAAAATCCCTCTCCCCACCGCCCCAATGCCAACGCCTTATCATTTTCCGTACCCTCgat              >  1:57248/1‑148 (MQ=255)                                                                                                                                          cgaGCTAAGAGAAAATTTATTCCGTTAATCACCC‑TTTATTTTTATTTTTTAAACGGGGGTTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAcc          <  1:170678/150‑1 (MQ=255)                                                                                                                                           gaGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCACGTTTTTTCATCATCCCTACCCCCGTTTAccg         >  2:425972/1‑150 (MQ=255)                                                                                                                                             gCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTTAACCGCCATTAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGTTTAccgcc       >  1:399957/1‑150 (MQ=255)                                                                                                                                               tAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCgg     >  2:400146/1‑150 (MQ=255)                                                                                                                                                 agagCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCCTACGCCCGTTTACCGCCGGGc   >  2:431716/1‑150 (MQ=255)                                                                                                                                                  gagCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCt  >  1:376525/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCC‑GTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCT  >  CP009273/2144964‑2145257

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