Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,258,531	C→G	7.0%	A65A (GCC→GCG)	asd →	aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,258,531	0	C	G	7.0%	211.3 / 3.1	68	A65A (GCC→GCG)	asd	aspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base G (3/2);  total (35/33)
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.67e-01

CCTTGGTGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAACGCCGAACG  >  NC_002947/2258392‑2258677                                                                                                                                              |                                                                                                                                                    ccTTGGTGCAGATCCTCGAAGAACTGGCATTCGCCGGTTGCCACTCTGCACCTGATGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTc                                                                                                                                               >  1:201593/1‑147 (MQ=255) ccTTGGTGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAg                                                                                                                                            >  5:172385/1‑149 (MQ=255)     ggTGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTgg                                                                                                                                        >  7:243385/1‑149 (MQ=255)     ggTGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTg                                                                                                                                         <  1:203868/148‑1 (MQ=255)      gTGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACATGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCTCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                                       >  8:83523/1‑149 (MQ=255)       tGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGcc                                                                                                                                      <  2:80485/149‑1 (MQ=255)       tGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGc                                                                                                                                       <  6:172385/148‑1 (MQ=255)                     gaaCTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCgg                                                                                                                        >  8:45050/1‑149 (MQ=255)                         tGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCct                                                                                                                                     >  4:50060/1‑132 (MQ=255)                         tGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCct                                                                                                                                     <  3:50060/132‑1 (MQ=255)                            cATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                                 <  8:23895/149‑1 (MQ=255)                            cATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                                 <  4:295086/149‑1 (MQ=255)                              ttC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGt                                                                                                               >  8:102812/1‑149 (MQ=255)                                         ccACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGc                                                                                                     >  7:259103/1‑149 (MQ=255)                                          cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                                    >  7:182718/1‑149 (MQ=255)                                           aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                                    >  1:116935/1‑148 (MQ=255)                                           aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                                    <  2:116935/148‑1 (MQ=255)                                            ctctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCttt                                                                                                  >  5:320006/1‑149 (MQ=255)                                               tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGcc                                                                                               >  5:234361/1‑149 (MQ=255)                                                      gctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAgg                                                                                        >  4:165406/1‑149 (MQ=255)                                                       ctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGgc                                                                                       <  7:83523/149‑1 (MQ=255)                                                          gCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCt                                                                                    <  6:183731/149‑1 (MQ=255)                                                            cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                                  >  6:9936/1‑149 (MQ=255)                                                            cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                                  <  6:17764/149‑1 (MQ=255)                                                            cAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                                  <  8:243385/149‑1 (MQ=255)                                                              gCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAg                                                                                >  5:171529/1‑149 (MQ=255)                                                               cATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAgg                                                                               <  6:320006/149‑1 (MQ=255)                                                               cATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAgg                                                                               >  2:292773/1‑149 (MQ=255)                                                                 tGGAATAGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑g                                                                            <  3:240352/149‑1 (MQ=255)                                                                  ggAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCGAGGCGCTGCAGGC‑gg                                                                           >  5:273703/1‑149 (MQ=255)                                                                   gAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑ggg                                                                          >  2:84999/1‑149 (MQ=255)                                                                              cagcagCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCt                                                                        <  4:245387/140‑1 (MQ=255)                                                                              cagcagCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCt                                                                        >  3:245387/1‑140 (MQ=255)                                                                                 cagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCa                                                                     >  6:157668/1‑140 (MQ=255)                                                                                 cagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCa                                                                     <  5:157668/140‑1 (MQ=255)                                                                                     gTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCt                                                        >  3:102811/1‑149 (MQ=255)                                                                                       gATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGt                                                      >  3:308914/1‑149 (MQ=255)                                                                                           ttCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCcg                                                   <  5:9936/148‑1 (MQ=255)                                                                                               cTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCggcggc                                              >  2:103019/1‑149 (MQ=255)                                                                                                  gCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTg                                           <  1:292773/149‑1 (MQ=255)                                                                                                    aaGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGga                                         <  7:45050/149‑1 (MQ=255)                                                                                                      gAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgacg                                       >  8:268873/1‑149 (MQ=255)                                                                                                         aCTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGacgacg                                    >  3:252240/1‑149 (MQ=255)                                                                                                             aaGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCGAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccct                                <  6:273703/149‑1 (MQ=255)                                                                                                             aaGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCa                                                                     <  1:314273/112‑1 (MQ=255)                                                                                                             aaGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCa                                                                     >  2:314273/1‑112 (MQ=255)                                                                                                               ggTGCGCGAAGTCGACAGCT‑TCGCCTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctgg                              >  2:179188/1‑149 (MQ=255)                                                                                                                ctgggcgAAGTCGA‑AGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggc                             >  6:218466/5‑149 (MQ=255)                                                                                                                 tgggcgAAGTCGA‑AGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGAc                                                          <  6:318327/116‑1 (MQ=38)                                                                                                                 tgggcgAAGTCGA‑AGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggcc                            >  4:306256/4‑149 (MQ=255)                                                                                                                 tgggcgAAGTCGA‑AGCTGTCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACCGTCATCGAc                                                          >  5:318327/4‑119 (MQ=38)                                                                                                                    gcgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg                         <  6:234361/149‑1 (MQ=255)                                                                                                                     cgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTg                                                       <  6:248187/118‑1 (MQ=255)                                                                                                                     cgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTg                                                       >  5:248187/1‑118 (MQ=255)                                                                                                                     cgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGACGGCCTGGACGACGccctggccctgg                        <  4:102811/149‑1 (MQ=255)                                                                                                                      gAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGccctggccctgg                        <  1:84999/149‑1 (MQ=255)                                                                                                                          tCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCg                   <  3:165406/149‑1 (MQ=255)                                                                                                                          tCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCg                   <  2:201593/149‑1 (MQ=255)                                                                                                                            gACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCgg                  <  4:308914/148‑1 (MQ=255)                                                                                                                             aCAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa                <  4:103944/149‑1 (MQ=255)                                                                                                                              cAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg               <  7:268873/149‑1 (MQ=255)                                                                                                                              cAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg               <  4:252240/149‑1 (MQ=255)                                                                                                                              aagctgtCGACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa                <  5:218466/143‑1 (MQ=255)                                                                                                                                gCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGcc             <  6:171529/149‑1 (MQ=255)                                                                                                                                 cT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGcc             >  1:93377/1‑148 (MQ=255)                                                                                                                                        cTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAACGcc       >  2:250165/1‑148 (MQ=255)                                                                                                                                            gCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAACGCCGAACg  >  6:224363/1‑149 (MQ=255)                                                                                                                                            gCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAACGCCGAACg  >  4:178313/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                    CCTTGGTGCAGATCCTCGAAGAACTGGCATTC‑CCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGC‑GGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAGCCAACGCCGAACG  >  NC_002947/2258392‑2258677

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