Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	436,896	A→T	5.9%	L131Q (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,896	0	A	T	5.9%	248.8 / 4.5	85	L131Q (CTG→CAG)	PP_0358	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base T (2/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CCAGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGC  >  NC_002947/436750‑437038                                                                                                                                                   |                                                                                                                                                ccaGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCgccagc                                                                                                                                               <  6:177746/149‑1 (MQ=255) ccaGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCgccagc                                                                                                                                               <  1:293525/149‑1 (MQ=255)      ccGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTg                                                                                                                                          <  4:41865/149‑1 (MQ=255)       cGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACCCCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTg                                                                                                                                          >  6:211128/1‑148 (MQ=255)        gggCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCa                                                                                                                                        <  7:73034/149‑1 (MQ=255)            agcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAggcg                                                                                                                                    >  2:199786/1‑149 (MQ=255)              cgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggt                                                                                                                                  <  3:243798/149‑2 (MQ=255)               gcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggtg                                                                                                                                 <  8:279613/149‑3 (MQ=255)                 agcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTCCTGCGCCACACACACGGCAGCCATGCCCACCGCCAGCCATTGCAGccggcg                                                                                                                                 >  3:216044/1‑147 (MQ=255)                     accgcatccTGCTTACGGATAAGAACAAACAACCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCCGgcg                                                                                                                           >  8:191847/8‑149 (MQ=255)                     aCCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCGCCGCCTGCAATGGCTggcggtgggcat                                                                                                                            <  5:143462/148‑8 (MQ=255)                     aCCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCGCCGCCTGCAATGGCTggcggtgggcat                                                                                                                            >  6:143462/1‑141 (MQ=255)                      ccGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGgcgc                                                                                                                          <  6:162551/149‑1 (MQ=255)                       cGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCt                                                                                                                         <  5:211128/149‑1 (MQ=255)                         gcgcCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCa                                                                                                                       <  2:239918/149‑1 (MQ=255)                           gcCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCAc                                                                                                                      >  5:160745/1‑148 (MQ=255)                             cTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                   <  4:80607/149‑1 (MQ=255)                                cTTACGGATCAGCCCATAGAATCCCAAGCTCAACGCCAGCACCACTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                 <  3:191605/148‑1 (MQ=255)                                 ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAAt                                                                                                               >  5:208705/1‑149 (MQ=255)                                    cGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                            >  2:6328/1‑149 (MQ=255)                                     ggATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                            >  8:212596/1‑148 (MQ=255)                                          aGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                      <  6:160745/149‑1 (MQ=255)                                            cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                      >  1:195646/1‑147 (MQ=255)                                            cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                      >  1:44957/1‑147 (MQ=255)                                            cccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGcc                                                                                                      >  6:28049/1‑147 (MQ=255)                                             ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGTCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCcag                                                                                                   >  2:185965/1‑149 (MQ=255)                                             ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                   >  2:308483/1‑149 (MQ=255)                                             ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                   >  6:215423/1‑149 (MQ=255)                                               aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcagca                                                                                                 <  1:199786/149‑1 (MQ=255)                                               aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                   >  8:76/1‑147 (MQ=255)                                               aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                   >  4:58286/1‑147 (MQ=255)                                                tAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGACCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCCGCATGCCcagcag                                                                                                >  8:267066/1‑149 (MQ=255)                                                 aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCATCAAGAATCAGCATGCCcagcag                                                                                                >  6:215705/1‑148 (MQ=255)                                                 aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                   >  1:43668/1‑129 (MQ=255)                                                 aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                   <  2:43668/129‑1 (MQ=255)                                                   aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTa                                                                                              >  3:170196/1‑148 (MQ=255)                                                        gAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGAGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCac                                                                                                            >  6:300138/1‑128 (MQ=255)                                                             tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                   >  4:142620/1‑149 (MQ=255)                                                             tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                   >  7:65578/1‑149 (MQ=255)                                                               aaCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCgg                                                                                 <  1:6328/149‑1 (MQ=255)                                                                aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCggg                                                                                <  2:44957/149‑1 (MQ=255)                                                                     aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                           <  3:58286/149‑1 (MQ=255)                                                                     aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                           <  5:117085/149‑1 (MQ=255)                                                                           aGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATa                                                                     <  3:142620/149‑1 (MQ=255)                                                                            gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                   <  5:215423/149‑1 (MQ=255)                                                                            gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                   <  5:28049/149‑1 (MQ=255)                                                                            gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑g                                                                   <  6:166386/149‑1 (MQ=255)                                                                             tGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑Gc                                                                  >  8:215810/1‑149 (MQ=255)                                                                                 aCCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCGCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAg                                                              >  3:216056/1‑149 (MQ=255)                                                                                   ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGc                                                            >  7:256468/1‑149 (MQ=255)                                                                                   ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGc                                                            <  7:138556/149‑1 (MQ=255)                                                                                       ggCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGc                                                        <  8:65578/149‑1 (MQ=255)                                                                                        gCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                   >  4:95049/1‑122 (MQ=255)                                                                                        gCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                   <  3:95049/122‑1 (MQ=255)                                                                                            gCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGa                                                   <  7:215810/149‑1 (MQ=255)                                                                                             cTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGTCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAg                                                  <  8:48120/149‑1 (MQ=255)                                                                                             cTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAg                                                  <  8:256468/149‑1 (MQ=255)                                                                                                cccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCa                                                <  7:76/148‑1 (MQ=255)                                                                                                 ccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCCAGCAGTACATTGATCAGCGGGTTGa                                                                            >  8:12900/1‑120 (MQ=255)                                                                                                 ccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCCAGCAGTACATTGATCAGCGGGTTGa                                                                            <  7:12900/120‑1 (MQ=255)                                                                                                   aCTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACg                                            <  1:308483/149‑1 (MQ=255)                                                                                                    cTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGcc                                                       <  3:779/137‑1 (MQ=255)                                                                                                    cTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGcc                                                       >  4:779/1‑137 (MQ=255)                                                                                                    cTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCACAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACgg                                           <  6:208705/149‑1 (MQ=255)                                                                                                      tgggcagccTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGcc                                         <  1:17796/142‑1 (MQ=255)                                                                                                       gCCACACCTGCTGCGCCACTCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGcc                                         <  4:216056/148‑1 (MQ=255)                                                                                                       gCCACACCTGCTGCGCCACACCCACGGCAGGCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGcc                                         <  7:267066/148‑1 (MQ=255)                                                                                                       gCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCg                                        >  3:286022/1‑149 (MQ=255)                                                                                                        ccACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGTCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAACATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATATGCCCAGGCTGGCCTCGAGCATACGGCCg                                        <  1:185965/149‑1 (MQ=255)                                                                                                          acacCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCgttg                                     <  4:170196/149‑1 (MQ=255)                                                                                                             ccTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCATCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTc                                  <  5:215705/149‑1 (MQ=255)                                                                                                              ctgctgCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCa                                 >  5:75460/1‑149 (MQ=255)                                                                                                               tgctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                >  8:102976/1‑149 (MQ=255)                                                                                                               tgctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                >  1:41930/1‑149 (MQ=255)                                                                                                                    cgcCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATc                                                                                                              >  7:274052/1‑67 (MQ=255)                                                                                                                    cgcCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATc                                                                                                              <  8:274052/67‑1 (MQ=255)                                                                                                                      ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                          >  6:98289/1‑148 (MQ=255)                                                                                                                      ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                          >  1:284885/1‑148 (MQ=255)                                                                                                                       cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                        >  1:108410/1‑149 (MQ=255)                                                                                                                        acacCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATACCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCacac                       >  5:28854/1‑149 (MQ=255)                                                                                                                            gcACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCTTACGGCCGTTGTTCACCGACCACACGTAg                   <  4:216044/148‑1 (MQ=255)                                                                                                                                ggCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGt                 >  1:94379/1‑147 (MQ=255)                                                                                                                                    gCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGTCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCa           <  7:191847/149‑1 (MQ=255)                                                                                                                                       aTGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGtt        <  4:286022/149‑1 (MQ=255)                                                                                                                                           ccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCg    >  2:101190/1‑149 (MQ=255)                                                                                                                                           ccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACAGACCACACGTAGGTCAGCCAGTTGCCg    >  8:137618/1‑149 (MQ=255)                                                                                                                                             gccGCCTGCAATGGCTGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGc  >  8:3546/2‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                CCAGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATA‑GCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGC  >  NC_002947/436750‑437038

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