Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,000,817	A→G	6.0%	Q451R (CAG→CGG)	PP_5256 →	cyclic nucleotide‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,000,817	0	A	G	6.0%	196.4 / 2.5	66	Q451R (CAG→CGG)	PP_5256	cyclic nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base G (2/2);  total (32/34)
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.22e-01

ACGAAGAAGATTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAAC  >  NC_002947/6000676‑6000965                                                                                                                                              |                                                                                                                                                      acgaAGAAGATTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCgg                                                                                                                                                <  3:305415/149‑1 (MQ=255) acgaAGAAGATTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCgg                                                                                                                                                <  2:103077/149‑1 (MQ=255)      gaagaTTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTc                                                                                                                                           >  7:54289/1‑149 (MQ=255)      gaagaTTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTc                                                                                                                                           >  6:216800/1‑149 (MQ=255)      gaagaTTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCAGGGAGTc                                                                                                                                           >  7:255674/1‑149 (MQ=255)         gaTTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAGCAGGTGAAGAGCTGTTTGCAGCAGCGTTGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCAGTTCAGGCGGCAGAGCCGTGAGTCgtt                                                                                                                                        >  7:59793/1‑149 (MQ=255)                tGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGttgttgttg                                                                                                                                 >  2:297404/1‑149 (MQ=255)                 ggCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGttgttgttg                                                                                                                                 >  4:27158/1‑148 (MQ=255)                     gAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGaa                                                                                                                            <  6:138076/149‑1 (MQ=255)                        ttCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGcc                                                                                                                         >  6:294048/1‑149 (MQ=255)                            gCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCt                                                                                                                     >  7:3586/1‑149 (MQ=255)                                    aCCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTtcaatc                                                                                                               >  5:6210/1‑147 (MQ=255)                                     ccAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTTCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCATAAGCCGGCTtcaatcaa                                                                                                             <  5:24695/148‑1 (MQ=255)                                        gTTGGGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCgtt                                                                                                                                        >  3:187192/1‑118 (MQ=255)                                        gTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCgtt                                                                                                                                        <  4:187192/118‑1 (MQ=255)                                             gTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTgg                                                                                                    >  4:137249/1‑149 (MQ=255)                                              tGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGc                                                                                                   >  6:62481/1‑149 (MQ=255)                                                cTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCg                                                                                                                        >  4:18105/1‑126 (MQ=255)                                                cTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCg                                                                                                                        <  3:18105/126‑1 (MQ=255)                                                   tACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTcc                                                                                               <  5:62481/148‑1 (MQ=255)                                                   tACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCt                                                                                              >  8:1006/1‑149 (MQ=255)                                                   tACCGGATCGACAAGGAACAGGTGAAGAGCTGTTGGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGCGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCAGGCTTCAATCAAGAAGGGTGGCTTCCt                                                                                              >  8:298084/1‑149 (MQ=255)                                                    aCCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCt                                                                                              >  6:211928/1‑148 (MQ=255)                                                     ccGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGa                                                                                            <  2:212922/149‑1 (MQ=255)                                                     ccGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGa                                                                                            <  8:3586/149‑1 (MQ=255)                                                       ggATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGc                                                                                                                          >  6:96342/1‑117 (MQ=255)                                                       ggATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGc                                                                                                                          <  5:96342/117‑1 (MQ=255)                                                          tCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCGGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttg                                                                                        <  3:120198/148‑1 (MQ=255)                                                          tCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCGGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttg                                                                                        <  3:115531/148‑1 (MQ=255)                                                          tCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCGGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttg                                                                                        >  4:115531/1‑148 (MQ=255)                                                          tCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCGGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttg                                                                                        >  4:120198/1‑148 (MQ=255)                                                            aaCAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttggtt                                                                                     <  1:4656/148‑1 (MQ=255)                                                             aCAAGGATAAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTTAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAg                                                                                                                              <  7:177436/107‑1 (MQ=255)                                                             aCAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttggttg                                                                                    >  5:181776/1‑149 (MQ=255)                                                             aCAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGAGAAGTGCAGACGGCGTTTAGCAAATTGCAGAGGTTCAGGCGGCAGAGCCGGGAGTAGTTGTTGTTGCAg                                                                                                                              >  8:177436/1‑107 (MQ=255)                                                              cAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAgttggttg                                                                                    >  6:307614/1‑148 (MQ=255)                                                               aaGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGca                                                                                  >  4:287001/1‑149 (MQ=255)                                                               aaGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGca                                                                                  >  2:274939/1‑149 (MQ=255)                                                                  gAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCgttg                                                                                                                                       <  2:218323/93‑1 (MQ=255)                                                                      aGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGa                                                                           <  8:54289/149‑1 (MQ=255)                                                                                  gTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGa                                                                                            <  4:179886/120‑1 (MQ=255)                                                                                  gTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGc                                                               >  1:168550/1‑149 (MQ=255)                                                                                   tttGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGc                                                               <  2:168550/148‑1 (MQ=255)                                                                                            gcGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑gg                                                    <  3:137249/149‑1 (MQ=255)                                                                                                   gAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAAcc                                             >  1:54741/1‑149 (MQ=255)                                                                                                       tGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGa                                                                           <  2:285693/116‑1 (MQ=255)                                                                                                       tGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGa                                                                           >  1:285693/1‑116 (MQ=255)                                                                                                         cAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTc                                        >  8:110150/1‑148 (MQ=255)                                                                                                           gACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTCGTTGCACAAGTCAGCCTGTGCCGGCCTCTTCGC‑GGCTAAACACGCTCATa                                     >  6:173118/1‑149 (MQ=255)                                                                                                              ggCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACaa                                  <  3:287001/149‑1 (MQ=255)                                                                                                                    gAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCACGAAGGCTGGCTTCCTGAGTTGGGTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAagtagt                             <  7:83543/148‑1 (MQ=255)                                                                                                                     aGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTg                            <  1:297404/148‑1 (MQ=255)                                                                                                                      gCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCa                          >  2:116612/1‑149 (MQ=255)                                                                                                                       cAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCATAGCGGGGTTAAACCCCCCCTAAAagtatttac                          >  3:164307/1‑144 (MQ=255)                                                                                                                          aTTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGc                      <  7:110150/149‑1 (MQ=255)                                                                                                                           ttGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGcc                     >  5:241936/1‑149 (MQ=255)                                                                                                                             gCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCt                    <  6:6210/148‑1 (MQ=255)                                                                                                                                gCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTTCAGAATCAGGCTTCAATCAAGAAGGTTGGCTTCATTATTTGTTTGCACAAGTGAGCCTGTGCCTGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTgaga                <  5:173118/149‑1 (MQ=38)                                                                                                                                gCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTAGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTgaga                <  8:255674/149‑1 (MQ=255)                                                                                                                                 cGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGTAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGagag               <  3:27158/149‑1 (MQ=255)                                                                                                                                      ccGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCt           <  2:54741/148‑1 (MQ=255)                                                                                                                                        ggcggcAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCtgtg        <  6:241936/149‑1 (MQ=255)                                                                                                                                         gcggcAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTgtgg       <  1:130790/149‑1 (MQ=255)                                                                                                                                             cAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGaa   <  5:211928/149‑1 (MQ=255)                                                                                                                                             cAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGaa   <  1:116612/149‑1 (MQ=255)                                                                                                                                             cAGAGCCGGGAGTCGTTATTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGaa   <  5:216800/149‑1 (MQ=255)                                                                                                                                              agagCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAAc  <  7:1006/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                      ACGAAGAAGATTCGCTGGCCGAATTCCGCACGCTGACCAGTTGCGTGCTGTACCGGATCGACAAGGAACAGGTGAAGAGCTGTTTGCAGCAGCGTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAAC  >  NC_002947/6000676‑6000965

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