Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	470,955	A→T	5.7%	M1K (ATG→AAG) †	PP_0386 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	470,955	0	A	T	5.7%	215.0 / 2.8	70	M1K (ATG→AAG) †	PP_0386	sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base T (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

GAACGGTGGCCAAGAGGCGTCTACGCCAAGTTTCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAG  >  NC_002947/470810‑471102                                                                                                                                                    |                                                                                                                                                    gAACGGTGGCCAAGAGGCGTCTACGCCAAGTTTCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATgg                                                                                                                                                  >  5:282757/1‑149 (MQ=255)    cGGTGGCCAATTTGCGTTTTCGCCAATTTTTAGGTGGGGGTGGGCATCCATCCCTCCGCGTACCCCCCCCATGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGc                                                                                                                                                 <  1:279957/149‑1 (MQ=255)            aaGAGGCGTCTACGCCAAGTTTCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCACCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCAc                                                                                                                                       >  6:176030/1‑149 (MQ=255)                          ccAAGTTTCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGaa                                                                                                                         <  8:203131/149‑1 (MQ=255)                               tttCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACTTGGCAGCCTCGCACTCAGGTGGTCTGAAGTc                                                                                                                      >  7:25887/1‑147 (MQ=255)                               tttCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACTTGGCAGCCTCGCACTCAGGTGGTCTGAAGTc                                                                                                                      >  7:30996/1‑147 (MQ=255)                                  cAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAgt                                                                                                                  <  8:67161/148‑1 (MQ=255)                                  cAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAgt                                                                                                                  <  8:131977/148‑1 (MQ=255)                                    gCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCt                                                                                                                                            <  4:307810/120‑1 (MQ=255)                                    gCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCt                                                                                                                                            >  3:307810/1‑120 (MQ=255)                                     cTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAg                                                                                                              <  3:232854/149‑1 (MQ=255)                                      tGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGa                                                                                                             >  5:256335/1‑149 (MQ=255)                                          ggTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACggg                                                                                                         >  8:40236/1‑149 (MQ=255)                                           gTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACggg                                                                                                         >  1:92369/1‑148 (MQ=255)                                                 agccagcc‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGcc                                                                                                  >  8:37214/1‑149 (MQ=255)                                                     agcc‑‑agGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAg                                                                                              <  2:92369/149‑1 (MQ=255)                                                       cc‑‑agGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGgt                                                                                            <  3:212350/149‑1 (MQ=255)                                                           agGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGgtgt                                                                                          >  5:182632/1‑149 (MQ=255)                                                             gCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAgtgt                                                                                                                >  7:40598/1‑125 (MQ=255)                                                             gCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAgtgt                                                                                                                <  8:40598/125‑1 (MQ=255)                                                                      tcATCTGTCAGGGTCAGGGCGCCGGCCTTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAgg                                                                               >  8:208575/1‑149 (MQ=255)                                                                      tcATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCAtt                                                                                                     >  5:116916/1‑127 (MQ=255)                                                                      tcATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCAtt                                                                                                     <  6:116916/127‑1 (MQ=255)                                                                        aTCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCg                                                                             <  6:282757/149‑1 (MQ=255)                                                                             tcagggtcagggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACTTGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGc                                                                          <  8:30996/147‑1 (MQ=255)                                                                             tcagggtcagggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACTTGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGc                                                                          <  8:25887/147‑1 (MQ=255)                                                                               agggtcagggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTa                                                                      <  6:182632/149‑1 (MQ=255)                                                                                gggtcagggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTAt                                                                     <  5:309655/149‑1 (MQ=255)                                                                                gggtcagggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTAt                                                                     <  1:235564/149‑1 (MQ=255)                                                                                       ggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACgg                                                              >  8:108473/1‑149 (MQ=255)                                                                                       ggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACgg                                                              <  6:256335/149‑1 (MQ=255)                                                                                       ggCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACgg                                                              >  4:165699/1‑149 (MQ=255)                                                                                          gccggccgTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGc                                                            >  4:57896/1‑148 (MQ=255)                                                                                                cgTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTa                                                       >  5:87316/1‑147 (MQ=255)                                                                                                    gtgtTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGaa                                                  >  2:327401/1‑148 (MQ=255)                                                                                                         ccACAAAATCAGCATCACCAACAGAATGGCCGGCTTCAAGGACAAGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGACTGTAGACGGGCAATGCCGAAGGAGAGCAACGGAAGGCGGGCATTATCTCACGGGCCACTAAAGAAAAAccc                                            >  8:232328/1‑149 (MQ=255)                                                                                                           acaAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCg                                          <  7:101704/149‑1 (MQ=255)                                                                                                           acaAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCg                                          >  5:249169/1‑149 (MQ=255)                                                                                                            caAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCgg                                         >  3:43050/1‑149 (MQ=255)                                                                                                             aaaaTCAGGATAAGCAACAGAATGGCCGCCTTCAAGGACATGGCAGCCACGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCgg                                         >  5:6414/1‑148 (MQ=255)                                                                                                                    ggATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGcc                                 >  5:85336/1‑149 (MQ=255)                                                                                                                     gATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTACCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCg                                >  1:330096/1‑149 (MQ=255)                                                                                                                      aTCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgg                               >  6:73667/1‑149 (MQ=255)                                                                                                                       tCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              <  7:37214/149‑1 (MQ=255)                                                                                                                       tCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              <  3:165699/149‑1 (MQ=255)                                                                                                                       tCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              <  5:176030/149‑1 (MQ=255)                                                                                                                       tCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCggg                              <  6:85336/149‑1 (MQ=255)                                                                                                                           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACgg                                                              >  4:254604/1‑113 (MQ=255)                                                                                                                           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACgg                                                              <  3:254604/113‑1 (MQ=255)                                                                                                                           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              >  1:54769/1‑145 (MQ=255)                                                                                                                           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              <  2:54769/145‑1 (MQ=255)                                                                                                                           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGttt                          >  8:269079/1‑149 (MQ=255)                                                                                                                           cAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCGGGGttt                          <  1:327401/149‑1 (MQ=255)                                                                                                                            aaCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                         <  7:208575/149‑1 (MQ=255)                                                                                                                            aaCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                         <  7:269079/149‑1 (MQ=255)                                                                                                                             aCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTg                        >  5:312076/1‑149 (MQ=255)                                                                                                                              cAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCCAGGCCGGGGtttttc                       <  3:57896/149‑3 (MQ=255)                                                                                                                               aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtc                      <  1:153371/149‑1 (MQ=255)                                                                                                                               aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtc                      <  4:43050/149‑1 (MQ=255)                                                                                                                                gAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtca                     <  7:232328/149‑1 (MQ=255)                                                                                                                                           ttCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACtcgtc            >  2:323396/1‑147 (MQ=255)                                                                                                                                             cAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGag         >  1:340615/1‑148 (MQ=255)                                                                                                                                              aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAgg                                                                                             <  6:347469/63‑1 (MQ=255)                                                                                                                                              aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              <  3:283280/126‑1 (MQ=255)                                                                                                                                              aaGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCggg                              >  4:283280/1‑126 (MQ=255)                                                                                                                                               aGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCTCTCGTCGaggaa      <  6:312076/149‑1 (MQ=255)                                                                                                                                                  aCATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                         <  4:290375/127‑1 (MQ=255)                                                                                                                                                  aCATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtttt                         >  3:290375/1‑127 (MQ=255)                                                                                                                                                  aCATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaag     >  6:123820/1‑147 (MQ=255)                                                                                                                                                   cATGGCAGCCTCGCAGTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGGCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAg  <  6:87316/149‑1 (MQ=255)                                                                                                                                                    aTGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaa      <  7:263930/144‑1 (MQ=255)                                                                                                                                                    aTGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaa      >  8:263930/1‑144 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GAACGGTGGCCAAGAGGCGTCTACGCCAAGTTTCAGCTGCGGGTGGGCAGCCAGCC‑‑AGGCGTTCTCCTCATCTGTCAGGGTCAGGGCGCCGGCCGTTGTGTTCCACAAAATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAG  >  NC_002947/470810‑471102

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