Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,905,625	A→T	6.7%	E627V (GAG→GTG)	nirB →	nitrite reductase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,905,625	0	A	T	6.7%	176.3 / 3.2	70	E627V (GAG→GTG)	nirB	nitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (2/3);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

GAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCC  >  NC_002947/1905483‑1905771                                                                                                                                               |                                                                                                                                                    gAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCgg                                                                                                                                               >  1:222137/1‑149 (MQ=255) gAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCgg                                                                                                                                               >  5:364622/1‑149 (MQ=255)  aaaTAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑t                                                                                                                                             <  5:23319/149‑1 (MQ=255)    aTAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑ttt                                                                                                                                           <  8:331691/149‑1 (MQ=255)        cccccGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAgg                                                                                                                                                   <  2:9762/138‑1 (MQ=255)        cccccGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAgg                                                                                                                                                   >  1:9762/1‑138 (MQ=255)               cAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGcc                                                                                                                                <  7:210599/149‑1 (MQ=255)               acaaatgatcgtcatcgGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGAGCGTGGCCGG‑TTTCGAAACCGGcc                                                                                                                                <  4:138872/142‑1 (MQ=255)                aaGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGcct                                                                                                                               <  1:60323/149‑2 (MQ=255)                aaGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCa                                                                                                                               >  4:251574/1‑149 (MQ=255)                      atcgtcatcgGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTa                                                                                                                         >  6:188039/1‑149 (MQ=255)                            atcgGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGc                                                                                                                     >  1:105457/1‑147 (MQ=255)                             tcgGCCAGGTGACCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCTTACGGCAAg                                                                                                                  >  1:230715/1‑149 (MQ=255)                               ggCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGt                                                                                                                 <  5:14860/148‑1 (MQ=255)                               ggCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGt                                                                                                                 <  3:302678/148‑1 (MQ=255)                                 ccAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGa                                                                                                              >  3:119315/1‑149 (MQ=255)                                  cAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGAc                                                                                                             <  2:222137/149‑1 (MQ=255)                                  cAGGTGGCCAAGAAATACCACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGAc                                                                                                             <  1:12136/149‑1 (MQ=255)                                       ggacAAGAAATCCCACCTGTACCCCATTCGCCCCGGTGGCCAGGGCATCGCCTTGTGCGCCGCCCAGTGCCACGCGCGGCCGCTTATCGGGGGCGAATTGAGCGTGGCCGG‑TTTCTAAACCGGCCACGCCTACGGCAAGTCGACCCGCa                                                                                                        <  4:173931/146‑1 (MQ=255)                                          cAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGAccc                                                                                                           >  8:206149/1‑143 (MQ=255)                                          cAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGAccc                                                                                                           <  7:206149/143‑1 (MQ=255)                                          cAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACgg                                                                                                     >  5:28329/1‑149 (MQ=255)                                              aaaTACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACgactt                                                                                                  >  1:272622/1‑144 (MQ=255)                                                    gACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTcgtg                                                                                           >  4:296443/1‑149 (MQ=255)                                                    gACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTcgtg                                                                                           <  6:364622/149‑1 (MQ=255)                                                             acCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAACCCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGgcag                                                                                  <  7:24005/149‑1 (MQ=255)                                                             acCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGAccc                                                                                                           <  4:260389/124‑1 (MQ=255)                                                             acCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGAccc                                                                                                           >  3:260389/1‑124 (MQ=255)                                                             acCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGgcag                                                                                  >  8:24005/1‑149 (MQ=255)                                                             acCAAGATCACCGGTGGCCAGCGCATCGACATGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGCGGCCGG‑TTTCGAAACCGGCCCCGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGgcag                                                                                  >  2:286086/1‑149 (MQ=255)                                                               cAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGgcagca                                                                                >  2:2567/1‑149 (MQ=255)                                                               cAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGACCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGgcagca                                                                                >  2:210729/1‑149 (MQ=255)                                                                  gATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCt                                                                             <  5:188039/149‑1 (MQ=255)                                                                   aTCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTg                                                                            <  2:105457/149‑1 (MQ=255)                                                                      aCCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTGCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTg                                                                         <  3:203904/149‑1 (MQ=255)                                                                      aCCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGACCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTg                                                                         <  1:210729/149‑1 (MQ=255)                                                                        cGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTgg                                                                           >  7:330360/1‑145 (MQ=255)                                                                        cGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTgg                                                                           <  8:330360/145‑1 (MQ=255)                                                                        cGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTg                                                                            <  6:24699/144‑1 (MQ=255)                                                                        cGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTg                                                                            >  5:24699/1‑144 (MQ=255)                                                                            ggCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTa                                                                   <  1:2567/149‑1 (MQ=255)                                                                               cAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACgg                                                                >  5:359590/1‑149 (MQ=255)                                                                                aGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGc                                                               <  1:232555/149‑1 (MQ=255)                                                                                aGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGc                                                               >  3:63154/1‑149 (MQ=255)                                                                                         aCCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGa                                                       >  8:268592/1‑148 (MQ=255)                                                                                          ccTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGTGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGc                                                                        >  7:49807/1‑130 (MQ=255)                                                                                          ccTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACa                                                     <  3:296443/149‑1 (MQ=255)                                                                                           cTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGTGGCCGG‑TTTCGAAACCGGcc                                                                                                                                <  5:18262/73‑1 (MQ=255)                                                                                           cTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGTGGCCGG‑TTTCGAAACCGGcc                                                                                                                                >  6:18262/1‑73 (MQ=255)                                                                                              ttCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAgtg                                                  <  4:119315/148‑1 (MQ=255)                                                                                                       cAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGc                                        >  7:100082/1‑149 (MQ=255)                                                                                                            gCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCCGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTgc                                   >  2:134719/1‑149 (MQ=255)                                                                                                            gCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACg                                                                 <  2:241611/119‑1 (MQ=255)                                                                                                            gCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACg                                                                 >  1:241611/1‑119 (MQ=255)                                                                                                             cACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTgcg                                  >  2:239522/1‑149 (MQ=255)                                                                                                             cACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTgcg                                  >  6:257943/1‑149 (MQ=255)                                                                                                                   cTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGTTTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTg                              >  1:239823/1‑148 (MQ=255)                                                                                                                    tGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTgg                             >  6:143346/1‑147 (MQ=255)                                                                                                                            aTCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCt                    <  1:239522/148‑1 (MQ=255)                                                                                                                            aTCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCt                    <  5:257943/148‑1 (MQ=255)                                                                                                                            aTCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTa                   >  4:328183/1‑149 (MQ=255)                                                                                                                                ggggCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTc                                                >  3:159484/1‑116 (MQ=255)                                                                                                                                ggggCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCAACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGCCGTGCAGGACAGTGACGCCATGCCACAGCGcc                                >  4:285199/1‑132 (MQ=255)                                                                                                                                 gggCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTTCGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAgg              >  1:292395/1‑149 (MQ=255)                                                                                                                                 gggCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTAc                  <  7:268592/145‑1 (MQ=255)                                                                                                                                   gCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGc            >  8:228600/1‑149 (MQ=255)                                                                                                                                      aaCTGATGGAGGCCGT‑TTTCGAAACGGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTg         <  6:28329/149‑1 (MQ=255)                                                                                                                                         tGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTgcg       >  2:107821/1‑148 (MQ=255)                                                                                                                                             cGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGcc  >  8:207256/1‑149 (MQ=255)                                                                                                                                              gAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGcc  <  2:292395/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    GAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGG‑TTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCC  >  NC_002947/1905483‑1905771

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