Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,148,793	G→T	8.7%	T151K (ACG→AAG)	PP_1007 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,148,793	0	G	T	8.7%	229.1 / 9.1	77	T151K (ACG→AAG)	PP_1007	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base T (2/5);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

CCTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCA  >  NC_002947/1148647‑1148940                                                                                                                                                    |                                                                                                                                                    ccTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGt                                                                                                                                                    >  3:97465/1‑148 (MQ=255) ccTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGt                                                                                                                                                    >  5:326512/1‑148 (MQ=255) ccTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGt                                                                                                                                                    >  6:348715/1‑148 (MQ=255)     ggggACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAggt                                                                                                                                               <  5:44720/149‑1 (MQ=255)      gggACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAggtg                                                                                                                                              <  7:77218/149‑1 (MQ=255)        gACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCTTACggtggt                                                                                                                                            <  7:54415/149‑1 (MQ=255)         aCCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTc                                                                                                                                           <  3:70865/149‑1 (MQ=255)         aCCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTc                                                                                                                                           <  5:294001/149‑1 (MQ=255)          ccTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCTTACGGTGGTCg                                                                                                                                          <  4:398329/149‑1 (MQ=255)          ccTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCg                                                                                                                                          <  6:326512/149‑1 (MQ=255)           cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCgg                                                                                                                                         >  2:236137/1‑149 (MQ=255)              gAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCt                                                                                                                                      <  1:365378/149‑1 (MQ=255)                  tAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAg                                                                                                                                  >  5:317693/1‑149 (MQ=255)                   aaGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTAG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAgg                                                                                                                                 >  1:380841/1‑149 (MQ=255)                    aGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                <  1:236137/149‑1 (MQ=255)                     gCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                <  1:72489/149‑1 (MQ=255)                     gCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgcg                                                                                                                               <  7:174101/149‑1 (MQ=255)                         cGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCTTACGGTGGTCGGc                                                                                                                                        <  6:193999/136‑1 (MQ=255)                         cGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCTTACGGTGGTCGGc                                                                                                                                        >  5:193999/1‑136 (MQ=255)                           cccTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTACAGGCGCATggg                                                                                                                         >  3:24520/1‑149 (MQ=255)                            ccTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATggg                                                                                                                         >  1:18717/1‑148 (MQ=255)                              tGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      >  8:169453/1‑149 (MQ=255)                               gTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGc                                                                                                                                    <  1:189924/134‑1 (MQ=255)                                tAGCCGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAAGGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      <  8:261467/147‑1 (MQ=255)                                      gggggggACGTGGGGGCGGTCGTGGGGTTGGCTGGCTAATGCCTAATTGGTTTTGAGAAGCATTTTTGCGCGGTTTTCCAGTTGCTGGCGATAG‑CGCTGGCCGACCACGGTATGGTGGTGGGCCTGCAGGCGCATGGGCAGGTTACCGa                                                                                                              <  2:380841/147‑1 (MQ=255)                                         ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACgg                                                                                                           >  2:205884/1‑149 (MQ=255)                                          gcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGt                                                                                                          <  4:97465/149‑1 (MQ=255)                                             aCCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGaa                                                                                                       >  7:393268/1‑149 (MQ=255)                                                  gcggcCGTCCTGGCGTTCGGTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGcc                                                                                                  <  6:317693/149‑1 (MQ=255)                                                  gcggcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAg                                                                                                                                                 >  8:266919/1‑102 (MQ=255)                                                  gcggcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAg                                                                                                                                                 <  7:266919/102‑1 (MQ=255)                                                  gcggcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGcc                                                                                                  <  3:270066/149‑1 (MQ=255)                                                   cggccggcCTGGCGTTCGCTGGCGAATGCCGAATGGGGGTTGAGCAGCACTTTGGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGccc                                                                                                 <  6:337358/149‑1 (MQ=255)                                                    ggcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGccc                                                                                                 >  8:304362/1‑148 (MQ=255)                                                      ccGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGAc                                                                                              >  7:334679/1‑149 (MQ=255)                                                      ccGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGAc                                                                                              <  7:169453/149‑1 (MQ=255)                                                              gCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcagcag                                                                                        >  1:229610/1‑147 (MQ=255)                                                                gTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCTTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                <  2:233349/105‑1 (MQ=255)                                                                gTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                >  1:233349/1‑105 (MQ=255)                                                                gTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCagcagcagc                                                                                    <  7:304362/149‑1 (MQ=255)                                                                     cTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGc                                                                               >  8:40326/1‑149 (MQ=255)                                                                      tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCATCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGc                                                                               >  3:121943/1‑148 (MQ=255)                                                                        gCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCATCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGc                                                                               <  4:121943/146‑1 (MQ=255)                                                                         ctaaTGCCTATTTGGTGTTGAGCAGCAGTTTTGCGCCGTCTTCCTTTTGCCGGCGCTGG‑CGCTCGCCGACCACCTTATGGTGGTGGGCCTACAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTg                                                                           <  4:24520/147‑1 (MQ=255)                                                                         cGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTg                                                                           <  8:334679/149‑1 (MQ=255)                                                                          gAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      >  2:47067/1‑105 (MQ=255)                                                                          gAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      <  1:47067/105‑1 (MQ=255)                                                                               ccGATTTGGTGTGGCCCACCACGTTTGCGCCGTCTGCCAGTCCCAGGCGCTGG‑CGCGCGCCGACCACCGTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAgg                                                                     <  7:139132/149‑1 (MQ=255)                                                                               ccGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAg                                                                      <  3:168586/148‑1 (MQ=255)                                                                                cGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                    <  7:40326/149‑1 (MQ=255)                                                                                cGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                    >  1:207726/1‑149 (MQ=255)                                                                                cGAATTGGTGTTGAGCAGCACTTTTGCGCCGGCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                    <  4:95096/149‑1 (MQ=255)                                                                                      ggTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTc                                                               <  1:205884/148‑1 (MQ=255)                                                                                       gtgtTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTc                                                               >  4:203421/1‑147 (MQ=255)                                                                                                agcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTgggtg                                                      >  1:351185/1‑143 (MQ=255)                                                                                                agcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTgcg                                                    <  2:18717/149‑1 (MQ=255)                                                                                                   aCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAg                                                 <  8:393268/149‑1 (MQ=255)                                                                                                    cTTTTGCGCCGTGTTCCAGTTGCCGGGGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAgg                                                <  2:207726/149‑1 (MQ=255)                                                                                                     ttttGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgc                                               >  3:203483/1‑149 (MQ=255)                                                                                                       ttGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgcgc                                             >  2:224082/1‑149 (MQ=255)                                                                                                          cgcCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCcg                                          >  8:183668/1‑149 (MQ=255)                                                                                                                  tCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGa                                  >  8:141827/1‑149 (MQ=255)                                                                                                                    cAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGAc                                 <  7:141827/148‑1 (MQ=255)                                                                                                                      gTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                    <  7:115911/111‑1 (MQ=255)                                                                                                                      gTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                    >  8:115911/1‑111 (MQ=255)                                                                                                                         gCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGTGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGt                           >  2:50932/1‑149 (MQ=255)                                                                                                                            ggcgctgg‑cgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTAACGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgc                        >  8:331906/1‑149 (MQ=255)                                                                                                                             gcgctgg‑cgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCAGGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgct                       <  4:96181/149‑1 (MQ=255)                                                                                                                                   g‑cgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTgg                  <  5:348715/148‑1 (MQ=255)                                                                                                                                           ccGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgc                                               <  1:65578/112‑1 (MQ=255)                                                                                                                                           ccGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgc                                               >  2:65578/1‑112 (MQ=255)                                                                                                                                              accaccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCg                                                                              >  6:308285/1‑78 (MQ=255)                                                                                                                                              accaccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCg                                                                              <  5:308285/78‑1 (MQ=255)                                                                                                                                                 accGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGc    >  6:133748/1‑149 (MQ=255)                                                                                                                                                  ccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGGGGCCAGGGGTTTCCAGTGACTGCGCCGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGcc   >  6:321336/1‑149 (MQ=255)                                                                                                                                                   cGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGGCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCa  >  1:128687/1‑149 (MQ=255)                                                                                                                                                   cGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCATGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCa  >  4:318811/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CCTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGG‑CGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCA  >  NC_002947/1148647‑1148940

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