Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	784,240	G→A	8.0%	G562G (GGC→GGT)	PP_0672 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	784,240	0	G	A	8.0%	176.4 / 5.1	62	G562G (GGC→GGT)	PP_0672	sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/3);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCTCTGTACCGAGGGGCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTG  >  NC_002947/784093‑784386                                                                                                                                                    |                                                                                                                                                    gCTCTGTACCGAGGGGCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGc                                                                                                                                                    <  8:339592/149‑1 (MQ=255)       tACCGAGGGGCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCacac                                                                                                                                              >  3:106255/1‑149 (MQ=255)          cGAGGGGCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCga                                                                                                                                           >  5:369225/1‑149 (MQ=255)             ggggCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGAGCTGAACGGTCACATCCTTGCCCACGCCgagga                                                                                                                                        >  8:151000/1‑149 (MQ=255)                gCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATgc                                                                                                                                     <  5:277250/149‑1 (MQ=255)                gCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATgc                                                                                                                                     <  6:369225/149‑1 (MQ=255)                  tGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATgcgc                                                                                                                                   <  4:291325/149‑1 (MQ=255)                      aTAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCACTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCa                                                                                                                               >  5:358344/1‑149 (MQ=255)                         gTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑a                                                                                                                           <  8:25593/149‑1 (MQ=255)                            gagcttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCAcc                                                                                                                             <  8:291335/145‑1 (MQ=255)                            gagcttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCAcc                                                                                                                             >  7:291335/1‑145 (MQ=255)                            gagcttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTACCCACACCGATGATGCGCGTCACC‑AGGc                                                                                                                        >  6:181474/1‑149 (MQ=255)                               cttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGccc                                                                                                                                                  >  1:2470/1‑121 (MQ=255)                               cttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGccc                                                                                                                                                  <  2:2470/121‑1 (MQ=255)                               cttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTACCCACACCGATGATGCGCGTCACC‑AGGCCtt                                                                                                                     <  5:181474/149‑1 (MQ=255)                               ctagagctGCTTGACGGTGGAAAAGAACACGTCGCTCATGCTTTCGGCGAGCAAGCGGTAGCGTTGCTCGCTGTCACGTACCGACTGGCAGGCCTCGATCTCAACGGTCACAACCTTGCCCACACCGATGAAGCGCGTCACC‑AGGccgt                                                                                                                     >  4:268452/4‑147 (MQ=255)                                ttgagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTc                                                                                                                    >  2:355470/1‑149 (MQ=255)                                  gagctGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCg                                                                                                                   <  6:126314/148‑1 (MQ=255)                                        gTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTg                                                                                                             <  7:7955/148‑1 (MQ=255)                                           gTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCGTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCgg                                                                                                          <  1:254189/148‑1 (MQ=255)                                           gTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCgg                                                                                                          <  4:113272/148‑1 (MQ=255)                                              ggTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCa                                                                                                      <  7:105185/149‑1 (MQ=255)                                                    aaaGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCt                                                                                                >  5:219710/1‑149 (MQ=255)                                                     aaGATCACGGCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTg                                                                                               <  7:136710/149‑1 (MQ=255)                                                        aTCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTc                                                                                            >  7:155460/1‑149 (MQ=255)                                                         tCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTcg                                                                                           <  8:73437/149‑1 (MQ=255)                                                          cACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTcgc                                                                                          >  3:288965/1‑149 (MQ=255)                                                           aCGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTcgcg                                                                                         >  7:122817/1‑149 (MQ=255)                                                           aCGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTcgc                                                                                          <  8:155460/148‑1 (MQ=255)                                                                gCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATg                                                                                     <  6:223376/148‑1 (MQ=255)                                                                    aTGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTACCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTa                                                                                <  5:9572/149‑1 (MQ=255)                                                                      gCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCg                                                                                                                   >  4:59666/1‑112 (MQ=255)                                                                      gCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCg                                                                                                                   <  3:59666/112‑1 (MQ=255)                                                                            ggCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACCAAGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAgcgc                                                                         <  1:271094/149‑1 (MQ=255)                                                                              cGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACCCCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAg                                                                      >  6:297291/1‑149 (MQ=255)                                                                                  cATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTc                                                                  >  3:233670/1‑149 (MQ=255)                                                                                         tAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTcgcg                                                           >  7:158523/1‑149 (MQ=255)                                                                                                  tCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCgg                                                    >  1:25061/1‑147 (MQ=255)                                                                                                    gCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTACCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAgc                                                >  8:86466/1‑149 (MQ=255)                                                                                                     cTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAgcg                                               <  1:355470/149‑1 (MQ=255)                                                                                                        tCACGTAGCGACTGGCTGGCCTCGAGCTGAAGGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGTGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAg                                            <  3:268452/149‑1 (MQ=255)                                                                                                          aCGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCt                                          >  6:127055/1‑149 (MQ=255)                                                                                                                cGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTagc                                                                              >  4:135932/1‑107 (MQ=255)                                                                                                                cGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTagc                                                                              <  3:135932/107‑1 (MQ=255)                                                                                                                      gctggcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCg                              >  6:123854/1‑149 (MQ=255)                                                                                                                       ctggcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCgg                             >  3:13098/1‑149 (MQ=255)                                                                                                                          gcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTg                          <  6:358344/149‑1 (MQ=255)                                                                                                                          gcCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTg                          <  4:106255/149‑1 (MQ=255)                                                                                                                          gcCTCGATCTGAACGGTCACATCCTTACCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGtt                           <  7:86466/148‑1 (MQ=255)                                                                                                                             tCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTc                                                                  >  1:2308/1‑106 (MQ=255)                                                                                                                             tCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTc                                                                  <  2:2308/106‑1 (MQ=255)                                                                                                                               gATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTc                       <  3:169148/147‑1 (MQ=255)                                                                                                                                 tCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAgcgc                                              <  4:108128/122‑1 (MQ=255)                                                                                                                                 tCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAgcgc                                              >  3:108128/1‑122 (MQ=255)                                                                                                                                    gAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAg                >  4:141751/1‑149 (MQ=255)                                                                                                                                     aaCGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAg                <  4:288965/148‑1 (MQ=255)                                                                                                                                      aCGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGtt              <  3:141751/149‑1 (MQ=255)                                                                                                                                       cGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCGGCTGGCAGTGCAGCGGTTGGTCATGACAGTTg             >  6:198508/1‑149 (MQ=255)                                                                                                                                          tCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTcg          <  5:123854/149‑1 (MQ=255)                                                                                                                                            acaTCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGAGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTcgcg        >  5:370649/1‑149 (MQ=255)                                                                                                                                                 cTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGtt   <  8:158523/149‑1 (MQ=255)                                                                                                                                                  ttGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTg  >  2:257559/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCTCTGTACCGAGGGGCTGACATAGTTGAGCTTGAGCTGGTTGTCGGTGGAAAAGATCACGTCGCTGATGCTTTCGGCGAGCATGCGGTAGCGTTGCTCGCTGTCACGTAGCGACTGGCTGGCCTCGATCTGAACGGTCACATCCTTGCCCACACCGATGATGCGCGTCACC‑AGGCCTTCGGCGTTGCGGGTCAGCACCTGCTCGCGGATGTCGTAGCAGCGCCAGCTTCCGTCGCGGTGGCGGAAGCGCAGCTGGCAGTGCAGCGGTTGGTCATGACAGTTGTCGCGCTGTTG  >  NC_002947/784093‑784386

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