Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	585,116	T→C	5.6%	Q47R (CAG→CGG) ‡	ansA ←	type 1 L‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	585,116	0	T	C	5.6%	274.4 / 3.6	91	Q47R (CAG→CGG) ‡	ansA	type 1 L‑asparaginase
Reads supporting (aligned to +/- strand):  ref base T (43/42);  new base C (2/3);  total (46/45)
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.20e-01

TGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAGCA  >  NC_002947/584970‑585263                                                                                                                                                    |                                                                                                                                                    tGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTgg                                                                                                                                                   >  6:25301/1‑149 (MQ=255)         ccGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCAtt                                                                                                                                           >  1:106869/1‑149 (MQ=255)          cGTGCAGCACCAACACGCCATCGTGGCGCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                          <  8:284917/149‑1 (MQ=255)          cGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTg                                                                                                                                          >  4:203549/1‑149 (MQ=255)          cGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCACCCCAGGTGCAATg                                                                                                                                          <  1:343345/149‑1 (MQ=255)           gTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTgg                                                                                                                                         <  6:171172/149‑1 (MQ=255)             gcagcaCCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGcc                                                                                                                                       <  3:210284/149‑1 (MQ=255)                gcaCCAACACGCCATCGTGGCCCGCGCCATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCTTATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCCCCTGGGGTGCATTGGCCagc                                                                                                                                    >  6:150676/1‑147 (MQ=255)                   ccAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTgg                                                                                                                                 >  1:173258/1‑149 (MQ=255)                    cAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTggg                                                                                                                                <  2:203539/149‑1 (MQ=255)                        acGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                            <  4:202889/149‑1 (MQ=255)                           ccATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAgtg                                                                                                                         >  1:229861/1‑149 (MQ=255)                             aTCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAgtgt                                                                                                                        >  2:172225/1‑148 (MQ=255)                              tCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTc                                                                                                                      <  5:25301/149‑1 (MQ=255)                               cGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCACCCCAGGTGCAATgga                                                                                                                                        >  1:218620/1‑129 (MQ=255)                               cGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCAATGCACCCCAGGTGCAATgga                                                                                                                                        <  2:218620/130‑2 (MQ=255)                                   gCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCTCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                            <  5:331744/138‑1 (MQ=255)                                   gCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                            >  6:331744/1‑138 (MQ=255)                                    cccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTCGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCAt                                                                                                                >  8:165940/1‑149 (MQ=255)                                     ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                            <  4:230013/136‑1 (MQ=255)                                     ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                            >  3:230013/1‑136 (MQ=255)                                     ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATc                                                                                                               >  6:388921/1‑149 (MQ=255)                                        cgACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGt                                                                                                            >  5:152207/1‑149 (MQ=255)                                         gACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGTAGCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGt                                                                                                            <  1:299089/149‑1 (MQ=255)                                              ccACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTc                                                                                                       >  2:321882/1‑148 (MQ=255)                                                aCCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaa                                                                                                     >  5:12288/1‑148 (MQ=255)                                                 ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                                   >  2:372082/1‑149 (MQ=255)                                                 ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                                   >  1:234766/1‑149 (MQ=255)                                                 ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                                   <  6:62341/149‑1 (MQ=255)                                                     cgcgACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAaccac                                                                                               <  8:118146/149‑1 (MQ=255)                                                     cgcgACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAaccac                                                                                               >  7:118146/1‑149 (MQ=255)                                                     cgcgACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCAAATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAaccac                                                                                               >  3:280514/1‑149 (MQ=255)                                                       cgacgaTGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCAccg                                                                                             >  4:271551/1‑149 (MQ=255)                                                             tGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGt                                                                                       <  1:328321/149‑1 (MQ=255)                                                                 atcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCa                                                                                   <  1:172225/149‑1 (MQ=255)                                                                 atcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCa                                                                                   <  2:229861/149‑1 (MQ=255)                                                                  tcgcccatagcCAGCCAGTTGTGCTGTTGCATATTGGCGATGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTTCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAg                                                                                  <  4:101680/149‑1 (MQ=255)                                                                    gcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAgg                                                                                 <  4:109791/148‑1 (MQ=255)                                                                      gcctcgcCAGCCAGGTGTGCTGTTGCATCTTTGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAg                                                                                  <  4:96344/142‑1 (MQ=255)                                                                         tcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTcc                                                                           >  2:183785/1‑149 (MQ=255)                                                                           gccagccaGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCgg                                                                         <  3:271551/149‑1 (MQ=255)                                                                           gccagccaGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCgg                                                                         <  2:173258/149‑1 (MQ=255)                                                                           gccagccaGTCGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCgg                                                                         <  7:165940/149‑1 (MQ=255)                                                                                ccagTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGttt                                                                    <  3:203549/149‑1 (MQ=255)                                                                                  agTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCg                                                                  <  6:12288/149‑1 (MQ=255)                                                                                   gTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGa                                                                 <  1:372082/149‑1 (MQ=255)                                                                                    ttGTGCTGTTGCAAATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAg                                                                <  4:280514/149‑1 (MQ=255)                                                                                      gtgCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                               <  6:152207/148‑1 (MQ=255)                                                                                       tgCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCAt                                                             >  3:184147/1‑149 (MQ=255)                                                                                       tgCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCAt                                                             <  5:388921/149‑1 (MQ=255)                                                                                              gCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCAACTTGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAACACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCAGGCCGATg                                                      >  5:264149/1‑149 (MQ=255)                                                                                               cATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTg                                                                                                           >  5:128263/1‑95 (MQ=255)                                                                                               cATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTg                                                                                                           <  6:128263/95‑1 (MQ=255)                                                                                                     ggCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCAc                                               >  4:166494/1‑149 (MQ=255)                                                                                                        gcTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                            <  2:234766/149‑1 (MQ=255)                                                                                                        gcTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCgg                                            >  3:235441/1‑149 (MQ=255)                                                                                                            tCGAGCAGCGGGTTCAGCTGGT‑TCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGta                                          <  6:264149/147‑1 (MQ=255)                                                                                                            tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                                      >  8:13533/1‑135 (MQ=255)                                                                                                            tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                                      <  7:13533/135‑1 (MQ=255)                                                                                                            tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGta                                          <  2:106869/147‑1 (MQ=255)                                                                                                            ctcgacagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTTTTGACGCCTCCGTTCTTCAAAACCACGGGCCGGTGCCAGGGCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGtata                                        >  1:253365/6‑149 (MQ=255)                                                                                                              gAGCAGCGGGTTCAGCTCCT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTAACTCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGtata                                        >  6:67270/1‑147 (MQ=255)                                                                                                                 cagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                               >  7:168327/1‑121 (MQ=255)                                                                                                                 cagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                               <  8:168327/121‑1 (MQ=255)                                                                                                                      ggTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACgaga                              <  1:183785/149‑1 (MQ=255)                                                                                                                        ttCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACgag                               >  3:350352/1‑146 (MQ=255)                                                                                                                         tCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACgag                               >  3:108624/1‑145 (MQ=255)                                                                                                                           aGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGt                                           <  5:386189/131‑1 (MQ=255)                                                                                                                           aGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGt                                           >  6:386189/1‑131 (MQ=255)                                                                                                                           aGCTCGT‑GCAGTGTCCATTGCACCATGGGTGCATTGGCCATCAGGGCAAAGTGTTCACGCATACGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATTCCGATGGTGCCACCGGTATAGAGAACGAGAAGAt                          >  7:51944/1‑148 (MQ=255)                                                                                                                            gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGc                                                                                >  7:308714/1‑93 (MQ=255)                                                                                                                            gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGc                                                                                <  8:308714/93‑1 (MQ=255)                                                                                                                            gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCTATGGTGCCACCGGTATagag                                     <  2:229512/136‑1 (MQ=255)                                                                                                                            gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATagag                                     >  1:229512/1‑136 (MQ=255)                                                                                                                                 t‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCGACCGGTATAGAGAACGAGAAGATTCTTTAc                   <  1:30746/149‑1 (MQ=255)                                                                                                                                       tgtCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgc              >  7:104502/1‑149 (MQ=255)                                                                                                                                         tCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgcgc            <  4:184147/149‑1 (MQ=255)                                                                                                                                          ccATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCa           >  7:88086/1‑149 (MQ=255)                                                                                                                                          ccAATGCACCCCAGGTGCAATGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGagaaga                           >  5:150011/1‑133 (MQ=255)                                                                                                                                          ccAATGCACCCCAGGTGCAATGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGagaaga                           <  6:150011/133‑1 (MQ=255)                                                                                                                                            aTTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgc                <  3:324706/142‑1 (MQ=255)                                                                                                                                            aTTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgc                >  4:324706/1‑142 (MQ=255)                                                                                                                                               gCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGa       >  7:157531/1‑148 (MQ=255)                                                                                                                                                cACCTGGTGTGCATTGGCCATCTGGTCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGGAt          <  3:106949/144‑1 (MQ=255)                                                                                                                                                cACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGACGATGCTCACACCGGTATAGAGAACGAGAAGATTCTTTACTCCGCGCAt          >  4:106949/1‑144 (MQ=255)                                                                                                                                                 aCCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgcg             <  5:6822/140‑1 (MQ=255)                                                                                                                                                 aCCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgcg             >  6:6822/1‑140 (MQ=255)                                                                                                                                                 aCCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGa       >  7:52639/1‑146 (MQ=255)                                                                                                                                                 aCCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAg    >  6:89668/1‑149 (MQ=255)                                                                                                                                                  ccTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAg    >  5:383645/1‑148 (MQ=255)                                                                                                                                                   cTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAGCa  <  6:383645/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGTCGGTACCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAGCA  >  NC_002947/584970‑585263

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