Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	436,905	A→T	6.4%	L128Q (CTG→CAG) ‡	PP_0358 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	436,905	0	A	T	6.4%	187.1 / 2.6	62	L128Q (CTG→CAG) ‡	PP_0358	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base T (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGCC  >  NC_002947/436758‑437048                                                                                                                                                     |                                                                                                                                                  ggCAGCGCAGCCACCGGCGCCTTCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTAGCTCCACACCCACGGCAGCCATGACCACCGCCAGCCATTGCa                                                                                                                                                   <  1:260360/149‑1 (MQ=255)     gcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAg‑gcgg                                                                                                                                             <  6:312368/149‑1 (MQ=255)     gcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAg‑gcgg                                                                                                                                             >  5:371727/1‑149 (MQ=255)       gcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTg‑gcggtg                                                                                                                                           <  5:332043/149‑3 (MQ=255)                   gcCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGTAG‑GCGGCGCAGGCGCTCAc                                                                                                                                >  6:59807/1‑148 (MQ=255)                        cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGaa                                                                                                                         >  5:337817/1‑149 (MQ=255)                        cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGa                                                                                                                          <  7:135227/148‑1 (MQ=255)                              gATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGAGCGGCGCAGGCGCTCACC‑AAGAATCAGc                                                                                                                    <  1:212746/149‑1 (MQ=255)                                    cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCca                                                                                                             <  6:40170/149‑1 (MQ=255)                                       aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCcagca                                                                                                          <  5:346650/149‑1 (MQ=255)                                       aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCcag                                                                                                            >  6:334661/1‑147 (MQ=255)                                           aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGTAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTa                                                                                                       >  3:118285/1‑148 (MQ=255)                                           aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTAc                                                                                                      <  3:7799/149‑1 (MQ=255)                                                    cTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGTAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAg                                                                                             <  4:118285/149‑1 (MQ=255)                                                      cAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCg                                                                                           >  5:104522/1‑149 (MQ=255)                                                             aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                    >  3:300850/1‑149 (MQ=255)                                                             aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                    <  3:310875/149‑1 (MQ=255)                                                               cACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGt                                                                                  <  4:300850/149‑1 (MQ=255)                                                                    gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAg                                                                             <  8:113047/149‑1 (MQ=255)                                                                     tGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGc                                                                            <  3:131989/149‑1 (MQ=255)                                                                         aCCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑TAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAg                                                                        >  1:305512/1‑149 (MQ=255)                                                                          cccACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                       >  4:98171/1‑149 (MQ=255)                                                                           ccACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGc                                                                      >  2:290277/1‑149 (MQ=255)                                                                           ccACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGAGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGc                                                                      >  5:152966/1‑149 (MQ=255)                                                                            cACGGCAGGTTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCt                                                                     >  2:23131/1‑149 (MQ=255)                                                                            cACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGACATGCCCACCGCCAGCCATTGCAG‑GCGCCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCACCGGGTTGATGTAATAGCCCAGGCt                                                                     >  7:344900/1‑149 (MQ=255)                                                                             aCGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGTAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTg                                                                    <  5:59807/149‑1 (MQ=255)                                                                             aCGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTg                                                                    <  1:290277/149‑1 (MQ=255)                                                                                    gCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCCGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTc                                                               >  3:200293/1‑147 (MQ=255)                                                                                    gCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGa                                                             <  1:169222/149‑1 (MQ=255)                                                                                     cTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGa                                                             <  6:337817/148‑1 (MQ=255)                                                                                         ccACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCTG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATa                                                        >  7:311776/1‑149 (MQ=255)                                                                                              tGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCTTGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGcc                                                   >  3:28463/1‑149 (MQ=255)                                                                                                  acacCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                               >  7:228786/1‑149 (MQ=255)                                                                                                  acacCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                               <  7:281366/149‑1 (MQ=255)                                                                                                   cacCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttgt                                              >  6:333543/1‑149 (MQ=255)                                                                                                     ccTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttgtt                                             >  5:201040/1‑148 (MQ=255)                                                                                                       tgct‑gcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                          >  2:80352/1‑149 (MQ=255)                                                                                                         ct‑gcGCCACACCCACTGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCTTGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAcc                                         <  6:371727/148‑1 (MQ=255)                                                                                                                 acacCCACGGCAGCCATGCCCACCGCCAGCCATTGCTG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCacac                                 <  8:311776/149‑1 (MQ=255)                                                                                                                 acacCCACGGCAGCCATGCCCACCGCCAGCCATTGCAT‑GCGGCGCAGGCGCTCAAC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCGAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCacac                                 <  3:6923/149‑1 (MQ=255)                                                                                                                 acacCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCacac                                 >  5:254715/1‑149 (MQ=255)                                                                                                                    cccACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACT‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACg                                >  6:175953/1‑147 (MQ=255)                                                                                                                     ccACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                             >  6:364080/1‑149 (MQ=255)                                                                                                                     ccACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                             <  6:270142/149‑1 (MQ=255)                                                                                                                     ccACGGCAGACATGCCCACCGCCAGCCATTGCAC‑GCGGCGCAGGCCCTCACC‑AACAATCAGCATGCCCACCAGAACATTCATCACCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATAAGGCCGTTGTTCACAGACCACACGAAg                             >  6:259200/1‑149 (MQ=255)                                                                                                                      cACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAgg                            >  8:54466/1‑149 (MQ=255)                                                                                                                       aCGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCAGACGGCCGTTGTTCACCGCCCACACGTAGGt                           >  3:301118/1‑149 (MQ=255)                                                                                                                            aGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACCGAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                       >  1:76910/1‑149 (MQ=255)                                                                                                                            aGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                       <  2:76910/148‑1 (MQ=255)                                                                                                                                aTGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCGCTCC‑TAGAAGCAGCATGCCCAGCAGTACAGTGATCGGCGGGGTGCTGTAATAGCCCAGGCTGGCCTCGCGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGtt                  >  3:212479/1‑149 (MQ=255)                                                                                                                                 cctgcgccGCCTGCAATGGCTG‑GCGGTGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTg                 >  2:316796/7‑149 (MQ=255)                                                                                                                                    ccACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCg              <  6:201040/149‑1 (MQ=255)                                                                                                                                    ccACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCg              >  6:367664/1‑149 (MQ=255)                                                                                                                                     cACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg             <  3:98171/149‑1 (MQ=255)                                                                                                                                      aCCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAcc                                         >  7:168626/1‑120 (MQ=255)                                                                                                                                      aCCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAcc                                         <  8:168626/120‑1 (MQ=255)                                                                                                                                        cGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGa          <  1:23131/149‑1 (MQ=255)                                                                                                                                           cagccaTTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt         >  2:17427/1‑147 (MQ=255)                                                                                                                                            agccaTTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAg      <  5:333543/149‑1 (MQ=255)                                                                                                                                              ccaTTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTg    <  1:80352/149‑1 (MQ=255)                                                                                                                                               caTTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGc   <  7:54466/149‑1 (MQ=255)                                                                                                                                                aTTCCAG‑GCGGCGCAGGCGCTCCCC‑AAGAATCAGCATGACCAGCAGTAAATTGATCAGCCGGATGATGTAATAGCCCAGGCTGGCCTCTCGCATACGGCCGTTCTTCCCCGACCACACCAAGCTCAGACACTTGCCGCCGAACAGTGcc  >  8:93935/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCT‑GCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAG‑GCGGCGCAGGCGCTCACC‑AAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGCC  >  NC_002947/436758‑437048

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