Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,559,215	C→A	5.8%	L104L (CTG→CTT)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,215	0	C	A	5.8%	193.7 / 2.9	69	L104L (CTG→CTT)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (39/26);  new base A (2/2);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

ATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAG  >  NC_000913/4559066‑4559364                                                                                                                                                      |                                                                                                                                                      aTCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAACGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGccc                                                                                                                                                       >  1:517995/1‑150 (MQ=255)  ccGCCACGTCTTTTACCACGGAACCCGTAATGGGGCGGGAAGGGACATGATAAGCGCCGGGCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGGCGGTGCCca                                                                                                                                                      <  1:887099/149‑1 (MQ=255)     ccccGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGGCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagca                                                                                                                                                   <  2:226068/147‑1 (MQ=255)     ccACGTCTTTTTCCACGGAACCCGAACTGGTGCGGGAAGGGACAGGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagca                                                                                                                                                   <  2:1340525/150‑1 (MQ=255)       aCGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGa                                                                                                                                                 <  2:242588/150‑1 (MQ=255)             ttttCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGaccaacc                                                                                                                                           >  2:230543/1‑150 (MQ=255)                 ccACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtg                                                                                                                                       >  2:518260/1‑150 (MQ=255)                         cccGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgtgacg                                                                                                                               >  2:854319/1‑150 (MQ=255)                             tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACAACTGACGTGACGCCCg                                                                                                                           >  1:658817/1‑150 (MQ=255)                                    gCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGCCGCCCGCTTCCGt                                                                                                                    >  1:191587/1‑150 (MQ=255)                                    gCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGt                                                                                                                    >  1:1104777/1‑150 (MQ=255)                                      gggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTcc                                                                                                                  >  2:129421/1‑149 (MQ=255)                                      gggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCa                                                                                                                  <  1:328441/150‑1 (MQ=255)                                      gggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCa                                                                                                                  <  1:416367/150‑1 (MQ=255)                                       ggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAg                                                                                                                 <  2:40190/150‑1 (MQ=255)                                       ggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCAAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAg                                                                                                                 >  1:625970/1‑150 (MQ=255)                                       ggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCAAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAg                                                                                                                 >  1:620226/1‑150 (MQ=255)                                        gAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAgg                                                                                                                <  2:455100/150‑1 (MQ=255)                                        gAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGGCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAgg                                                                                                                <  1:359141/150‑1 (MQ=255)                                             gACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACt                                                                                                           >  1:1043204/1‑150 (MQ=255)                                             gACATGATAAGCGCCGGTCAGCAACCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCACCAGACCAACAACTGACGAGACGCCCTCATCCGTCAGGCGACt                                                                                                           >  2:416768/1‑150 (MQ=255)                                                 tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCCGGCGACTTAgc                                                                                                       >  2:1296242/1‑150 (MQ=255)                                                 tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                       >  2:756047/1‑150 (MQ=255)                                                 tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                       >  2:667580/1‑150 (MQ=255)                                                 tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                       >  1:92503/1‑150 (MQ=255)                                                 tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                       >  2:608854/1‑150 (MQ=255)                                                 tGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCATCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                       <  2:1276574/150‑1 (MQ=255)                                                   aTAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGCCCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgcgc                                                                                                     >  2:110708/1‑150 (MQ=255)                                                   aTAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGCCGTGACGCCCGCTTCCGTAAGGCGACTTAgcgc                                                                                                     >  1:185748/1‑150 (MQ=255)                                                     aaGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGccc                                                                                                   >  1:190126/1‑149 (MQ=255)                                                       gcgcCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCCCt                                                                                                 >  1:951422/1‑150 (MQ=255)                                                            ggTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCgg                                                                                            >  1:1197577/1‑150 (MQ=255)                                                                     ccAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACAAACCACTTCCGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                                                                                   >  2:1259416/1‑150 (MQ=255)                                                                     ccAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGGGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCCCTTCCGGCGTGCGcgt                                                                                   >  1:144118/1‑150 (MQ=255)                                                                      cAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgtc                                                                                  >  1:721336/1‑150 (MQ=255)                                                                       aGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgtcg                                                                                 >  1:1352430/1‑150 (MQ=255)                                                                         gcgcTGATGCCTTCTTCTTTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTg                                                                               <  2:1258501/150‑1 (MQ=255)                                                                             tGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCACCCACTGACGTGACGCCCGCTTCCGTCCGGCGACTTAGCGCCCCTTCCGGCGTGCGCGTCGTGGGAc                                                                           >  2:646767/1‑150 (MQ=255)                                                                              gATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCACCCACTGACGTGACGCCCGCTTCCGTCCGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAcc                                                                          >  2:1113035/1‑150 (MQ=255)                                                                              gATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAcc                                                                          >  2:726827/1‑150 (MQ=255)                                                                              gATGCCTTCTTCATTGAGCGCACGGGGCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAcc                                                                          <  2:1174993/150‑1 (MQ=255)                                                                                  ccTTCTTCATTGCGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCt                                                                      <  1:292075/150‑1 (MQ=255)                                                                                  ccTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCt                                                                      >  2:53977/1‑150 (MQ=255)                                                                                   cttcttCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCtt                                                                     <  1:577891/150‑1 (MQ=255)                                                                                   cttcttCATTGAGCGCACGGGTCTTGGCGAGCAGGGAATCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCtt                                                                     <  1:663864/150‑1 (MQ=255)                                                                                    ttcttcATTGAGCACACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTc                                                                    <  2:331075/150‑1 (MQ=255)                                                                                          aTTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCCCTGACGTGACGCCCGCTTCCGTCCGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                                                              >  1:1213338/1‑150 (MQ=255)                                                                                          aTTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                                                              <  2:47477/150‑1 (MQ=255)                                                                                          aTTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTACGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                                                              >  1:1207937/1‑150 (MQ=255)                                                                                                cgcACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATATAGTCGGTGCCCAGCAGACCAACCACTGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGccgccaccgcc                                                        <  2:647756/150‑1 (MQ=255)                                                                                                      ggTCTTGGCGAGCAGCGATTCCCGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACAAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGCGGGACCTGCTTCGCCGCCACCGCCAATCaa                                                  <  1:1137169/150‑1 (MQ=255)                                                                                                        tCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAAt                                                >  1:749159/1‑150 (MQ=255)                                                                                                        tCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAAt                                                <  2:1352430/150‑1 (MQ=255)                                                                                                          ttGGCGAGCAGGGATTCCTGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCTCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATgg                                              <  2:1322365/150‑1 (MQ=255)                                                                                                            ggCGCGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGCCCTGCTTCGCCGCCCCCGCCAATCAAATGGAc                                            >  2:72149/1‑150 (MQ=255)                                                                                                            ggCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGGCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGAc                                            <  2:1116206/150‑1 (MQ=255)                                                                                                             gCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCAAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACg                                           <  2:625970/150‑1 (MQ=255)                                                                                                             gCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCAAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACg                                           <  2:620226/150‑1 (MQ=255)                                                                                                                 gCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGtt                                       >  1:777137/1‑150 (MQ=255)                                                                                                                  cAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCAGCTTCGCCGCCACCGCCAATCAAAAGGACGTGTTg                                      >  1:1128063/1‑150 (MQ=255)                                                                                                                         tCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATa                               <  2:115819/150‑1 (MQ=255)                                                                                                                         tCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATa                               <  2:622385/150‑1 (MQ=255)                                                                                                                           cGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGAAGTGTTGATCAATaaa                             >  2:1134057/1‑150 (MQ=255)                                                                                                                              ggtgctagagaTAGAGTCGGTGCCCAGCAGACCAAGCACTGACGTGTCGCCCGCTTCCGTCAGGCGACTTAGCGACACTTCCGGCGTGCGCTTCGTGGGACCTGCTTCGCCGCCACAGCCAATCAAATTTACGTGTTGATCAATAAAAcc                          <  1:75068/144‑1 (MQ=255)                                                                                                                              gTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAAcc                          >  2:397190/1‑150 (MQ=255)                                                                                                                              gTGGCGAGAGATAGAGTCGGGGCCCAGCAGACCAACCACTGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAAcc                          <  2:509321/150‑1 (MQ=255)                                                                                                                                ggCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTg                        <  1:756047/150‑1 (MQ=255)                                                                                                                                 gCGAGAGATATAGTCGGTGCCCAGCAGACCAACCTCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTgg                       <  2:228546/150‑1 (MQ=255)                                                                                                                                  cGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTggg                      >  1:907080/1‑150 (MQ=255)                                                                                                                                          agagTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGAt              >  2:312071/1‑150 (MQ=255)                                                                                                                                           gagTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATc             <  2:1197577/150‑1 (MQ=255)                                                                                                                                                ggTGCCCAGCAGACAAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTGGATCAATAAAACCTGGGCAGAGGATCTGccc        >  1:524183/1‑150 (MQ=255)                                                                                                                                                   gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCAcg     >  2:660430/1‑149 (MQ=255)                                                                                                                                                    cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACCTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTg    >  1:63883/1‑150 (MQ=255)                                                                                                                                                    cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGAATTAGCGCCACTTCCGGCGTGCGCGACGTGGGACCTGCTTCGCCGCCACCGCAAATCAAAAGGACAAGTTGATCAATAAAACCTCGGCAGAGGAACTCCCCACtt    >  2:1294345/1‑149 (MQ=255)                                                                                                                                                     ccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgt   >  1:52736/1‑149 (MQ=255)                                                                                                                                                      cagcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCACCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgag  >  2:536623/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      ATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAG  >  NC_000913/4559066‑4559364

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