Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,597,389	A→G	6.0%	V807A (GTC→GCC)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,597,389	0	A	G	6.0%	233.7 / 3.5	83	V807A (GTC→GCC)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base G (3/2);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGC  >  NC_002947/4597241‑4597531                                                                                                                                                     |                                                                                                                                                cTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGa                                                                                                                                                 >  6:218992/1‑149 (MQ=255) cTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGa                                                                                                                                                 <  1:74829/149‑1 (MQ=255)  tCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGAc                                                                                                                                                >  2:63614/1‑149 (MQ=255)    aGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGc                                                                                                                                              <  4:152964/149‑1 (MQ=255)     ggCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGAATCGTGACGcc                                                                                                                                             <  6:60671/149‑1 (MQ=255)          gCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGa                                                                                                                                         <  1:250188/148‑1 (MQ=255)          gCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGa                                                                                                                                         <  6:233407/148‑1 (MQ=255)           cTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATc                                                                                                                                       <  3:242376/149‑1 (MQ=255)            tGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGa                                                                                                                                         >  1:1867/1‑146 (MQ=255)            tGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGa                                                                                                                                         <  2:1867/146‑1 (MQ=255)            tGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATc                                                                                                                                       >  5:277038/1‑148 (MQ=255)                gATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGt                                                                                                                                  >  5:184905/1‑149 (MQ=255)                gATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGt                                                                                                                                  >  1:406963/1‑149 (MQ=255)                gATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCATGGCGTCACGATccttg                                                                                                                                   >  2:344870/1‑145 (MQ=255)                gATCGAAACGGGAGAAGGGGCGTCATCCGCCACGCGAAAATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCGAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCGAAGTGTTGAAGGATCGTGAGGCCATGATCCAgtt                                                                                                                                  >  2:65165/1‑147 (MQ=255)                 aTCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGATCTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGt‑c                                                                                                                                >  6:106243/1‑149 (MQ=255)                  tCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGt‑ct                                                                                                                               >  3:321426/1‑149 (MQ=255)                  tCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGt‑ct                                                                                                                               >  2:240767/1‑149 (MQ=255)                  tCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCATGGCGTCACGATccttgaa                                                                                                                                 <  1:344870/148‑6 (MQ=255)                     aaaCGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCg                                                                                                                             >  8:330343/1‑148 (MQ=255)                     aaaCGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCg                                                                                                                             >  1:217895/1‑148 (MQ=255)                     aaaCGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCg                                                                                                                             >  1:101405/1‑148 (MQ=255)                      aaCGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGaa                                                                                                                           <  5:218992/149‑1 (MQ=255)                          ggCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATc                                                                                                                       <  4:93321/149‑1 (MQ=255)                           gCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATcc                                                                                                                      <  8:71591/149‑1 (MQ=255)                                ggggCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCAc                                                                                                                 <  6:184905/149‑1 (MQ=255)                                  ggCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACgg                                                                                                               <  2:217895/149‑1 (MQ=255)                                            gccacgccaACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGt                                                                                                      >  8:383128/1‑148 (MQ=255)                                             ccacgccaACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTc                                                                                                    <  4:123217/149‑1 (MQ=255)                                                     aCATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAggcgcg                                                                                            <  4:60803/149‑1 (MQ=255)                                                           aCCGGTGATAACAGCGATCTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGcc                                                                                       <  5:106243/148‑1 (MQ=255)                                                           aCCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGcc                                                                                       <  1:351712/148‑1 (MQ=255)                                                            ccGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGTGCCAg                                                                                     <  7:330343/149‑1 (MQ=255)                                                            ccGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAg                                                                                     >  7:104743/1‑149 (MQ=255)                                                            ccGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAg                                                                                     <  7:292951/149‑1 (MQ=255)                                                                  aTAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCg                                                                               <  2:60927/149‑1 (MQ=255)                                                                    aaCAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAggcgc                                                                                             <  4:57521/133‑1 (MQ=255)                                                                    aaCAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAggcgc                                                                                             >  3:57521/1‑133 (MQ=255)                                                                     aCAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGAt                                                                             <  8:104743/148‑1 (MQ=255)                                                                     aCAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAACGATCGTGACGCCATGATCCAGGT‑CTCGAATATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGGCGATg                                                                            <  1:65165/149‑1 (MQ=255)                                                                      cAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCa                                                                                                   >  6:351863/1‑125 (MQ=255)                                                                      cAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCa                                                                                                   <  5:351863/125‑1 (MQ=255)                                                                       aGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATg                                                                            >  1:52470/1‑147 (MQ=255)                                                                                ttgtCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGCCGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCt                                                                 >  8:303167/1‑149 (MQ=255)                                                                                   tCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTa                                                              >  3:338683/1‑149 (MQ=255)                                                                                        aaCCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGt                                                         <  2:406963/149‑1 (MQ=255)                                                                                           cGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTgc                                                      >  7:133329/1‑149 (MQ=255)                                                                                            gTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTgcg                                                     <  1:63614/149‑1 (MQ=255)                                                                                              ctcCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGGTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTgcg                                                     >  3:47309/1‑147 (MQ=255)                                                                                               tcCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAg                                                  >  6:403662/1‑149 (MQ=255)                                                                                               tcCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAg                                                  >  5:52243/1‑149 (MQ=255)                                                                                                ccAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGGCGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGt                                                                                                      <  1:124647/96‑1 (MQ=255)                                                                                                ccAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGGCGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGt                                                                                                      >  2:124647/1‑96 (MQ=255)                                                                                                ccAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAgg                                                 >  7:14224/1‑149 (MQ=255)                                                                                                ccAGCAATCGGTGCAAAGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑ATCGAAAATCCGCCACGGAGAATTCGTTCAAGGCGCGGCGCCAGCTGTCTATGTAGTCTGCGCTATAAAGGTTTCGAAgg                                                 >  8:306029/1‑149 (MQ=255)                                                                                                 cAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGgc                                                <  7:383128/149‑1 (MQ=255)                                                                                                   gCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGgcgc                                              >  4:233524/1‑149 (MQ=255)                                                                                                        cGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAg                                                  >  1:99394/1‑140 (MQ=255)                                                                                                        cGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAg                                                  <  2:99394/140‑1 (MQ=255)                                                                                                           tGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAg                                      >  2:103735/1‑149 (MQ=255)                                                                                                             cAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAgcg                                    >  8:293310/1‑149 (MQ=255)                                                                                                                  gcgcCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGcc                               <  6:277038/149‑1 (MQ=255)                                                                                                                  gcgcCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGcc                               <  8:133329/149‑1 (MQ=255)                                                                                                                      cgcTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTAAAGGTTGCGAAGGCGCTCGATCAGCGCAGCCTGGt                           >  4:337351/1‑149 (MQ=255)                                                                                                                       gcTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTc                          >  4:50987/1‑149 (MQ=255)                                                                                                                       gcTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTc                          <  2:290513/149‑1 (MQ=255)                                                                                                                           gCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCa                      >  1:16198/1‑149 (MQ=255)                                                                                                                              cGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCg                   <  7:293310/149‑1 (MQ=255)                                                                                                                               ggCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCg                   >  6:258144/1‑148 (MQ=255)                                                                                                                                gCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGa                 <  2:52470/149‑1 (MQ=255)                                                                                                                                    aCTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCTGa                 <  7:374737/145‑1 (MQ=255)                                                                                                                                    aCTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGa                 >  8:374737/1‑145 (MQ=255)                                                                                                                                     cTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAAt            <  3:257906/149‑1 (MQ=255)                                                                                                                                        ttCAAGGATCGTGACGCCATGATCCAGGTGCTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATc           >  1:30689/1‑148 (MQ=255)                                                                                                                                        ttCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATc           >  7:74134/1‑147 (MQ=255)                                                                                                                                         acctggATCATGGCGTCACGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCTGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAg        >  7:214542/5‑149 (MQ=255)                                                                                                                                           aaGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGtt      >  7:230393/1‑149 (MQ=255)                                                                                                                                           aaGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGtt      >  6:289315/1‑149 (MQ=255)                                                                                                                                           aaGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGtt      <  2:16198/149‑1 (MQ=255)                                                                                                                                           aaGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGtt      <  1:240767/149‑1 (MQ=255)                                                                                                                                            aGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGtt      >  8:185460/1‑148 (MQ=255)                                                                                                                                             ggATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGa    >  2:229272/1‑149 (MQ=255)                                                                                                                                             ggATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGa    <  5:289315/149‑1 (MQ=255)                                                                                                                                              gATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAg   >  5:165952/1‑149 (MQ=255)                                                                                                                                               aTCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGTAGAATTCGTTCATGGCGCGGCGCCAGCTGTCGATGTAGGCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGc  <  7:306029/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CTCAGGCTTGCTGTTGATCGAAACGGGCGAAGGGGCCTCATCCGCCACGCCAACATTAACCGGTGATAACAGCGAACTGTTGTCGCGAACCGTCTCCAGCAATCGGTGCAACGGCGCCGCTGGCCCGGCCAACTGTTCAAGGATCGTGACGCCATGATCCAGGT‑CTCGAAAATCCGCCACGGAGAATGCGTTCAAGGCGCGGCGCCAGCTGTCGATGTAGTCGGCGCTGTACAGGTTGCGCAGGCGCTCGTTCAGCGCAGCCTGGTCGGCATCGGAGTAATCCAGTTGAGC  >  NC_002947/4597241‑4597531

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