Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	585,125	T→A	5.6%	N44I (AAT→ATT)	ansA ←	type 1 L‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	585,125	0	T	A	5.6%	277.8 / 3.2	89	N44I (AAT→ATT)	ansA	type 1 L‑asparaginase
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base A (2/3);  total (44/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.53e-01

CCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAGCATCCGTAGCGG  >  NC_002947/584978‑585273                                                                                                                                                     |                                                                                                                                                     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                                                                                                                                      >  3:230013/1‑136 (MQ=255)                             ccGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCaaa                                                                                                                                      <  4:230013/136‑1 (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                                                                                                             >  1:234766/1‑149 (MQ=255)                                         ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                                             <  6:62341/149‑1 (MQ=255)                                         ccgccgCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCaaaa                                                                                                             >  2:372082/1‑149 (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                                                                                             <  2:229861/149‑1 (MQ=255)                                                         atcgcgcatcgcCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCa                                                                                             <  1:172225/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                                                                <  7:13533/135‑1 (MQ=255)                                                                                                    tCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATg                                                                >  8:13533/1‑135 (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                                                                         <  8:168327/121‑1 (MQ=255)                                                                                                         cagcGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGc                                                                         >  7:168327/1‑121 (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                                                                                          <  8:308714/93‑1 (MQ=255)                                                                                                                    gCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGc                                                                                          >  7:308714/1‑93 (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                                     <  6:150011/133‑1 (MQ=255)                                                                                                                                  ccAATGCACCCCAGGTGCAATGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGagaaga                                     >  5:150011/1‑133 (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                       >  6:6822/1‑140 (MQ=255)                                                                                                                                         aCCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTgcgcg                       <  5:6822/140‑1 (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)                                                                                                                                                     ttGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAGCATCCGTAGCgg  >  5:236918/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGTGCAGCACCAACACGCCATCGTGGCCCGCGACATCCACCGCCGCGACGATGGCATCGCGCATCGCCAGCCAGTTGTGCTGTTGCATATTGGCGCTGTCGAGCAGCGGGTTCAGCTCGT‑GCAGTGTCCATTGCACCTGGGGTGCATTGGCCATCTGGGCAAAGTGTTCACGCATCCGTGCTTCAAAACCACCGGCCGGTGCCAGGCCTTCCGGGGTTTCGAGCATGCCGATGGTGCCACCGGTATAGAGAACGAGAAGATTCTTTACTGCGCGCATGGAAAGCATCCGTAGCGG  >  NC_002947/584978‑585273

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: ‑