Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,379,889	G→T	6.3%	L213M (CTG→ATG)	PP_3856 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,379,889	0	G	T	6.3%	192.6 / 3.5	64	L213M (CTG→ATG)	PP_3856	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base T (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

CAACCGGCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCG  >  NC_002947/4379742‑4380033                                                                                                                                                     |                                                                                                                                                 cAACCGGCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAg                                                                                                                                                  <  7:305683/148‑1 (MQ=255)  aaCCTGCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGc                                                                                                                                                <  7:146428/149‑1 (MQ=255)  aaCCGGCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTCTCAGACAGCGCGCTGTTGAACCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGAAGTGAC‑AAGGc                                                                                                                                                >  7:75797/1‑149 (MQ=255)       gCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCAcc                                                                                                                                            >  8:170882/1‑148 (MQ=255)       gCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCAca                                                                                                                                            >  5:152986/1‑147 (MQ=255)              tCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAAc                                                                                                                                      >  8:101665/1‑147 (MQ=255)                ttGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAgg                                                                                                                                  >  3:119400/1‑149 (MQ=255)                  gTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGcc                                                                                                                                >  8:49088/1‑149 (MQ=255)                   tttAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCt                                                                                                                               >  2:125577/1‑149 (MQ=255)                     tagtagAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTgg                                                                                                                             >  5:329291/1‑149 (MQ=255)                      agtagAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTgg                                                                                                                             <  4:72943/148‑1 (MQ=255)                          gagagTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CATGCCACCCATAACCAGGCCTGGGGc                                                                                                                          >  8:195832/1‑147 (MQ=255)                          gagagTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CATGCCACCCATAACCAGGCCTGGGGc                                                                                                                          >  8:22472/1‑147 (MQ=255)                                ttttGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAg                                                                                                                       >  8:268319/1‑144 (MQ=255)                                ttttGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAg                                                                                                                       <  7:268319/144‑1 (MQ=255)                                     aCAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTc                                                                                                             >  6:77255/1‑149 (MQ=255)                                      cAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCt                                                                                                            >  6:222289/1‑149 (MQ=255)                                       aaGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAgg                                                                                                                                  <  5:172289/126‑1 (MQ=255)                                       aaGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAgg                                                                                                                                  >  6:172289/1‑126 (MQ=255)                                       aaGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTaga                                                                                                                >  3:60188/1‑144 (MQ=255)                                       aaGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTaga                                                                                                                <  4:60188/144‑1 (MQ=255)                                       aaGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTAGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTgggg                                                                                                                           >  5:203034/1‑133 (MQ=255)                                       aaGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTAGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTgggg                                                                                                                           <  6:203034/133‑1 (MQ=255)                                            gCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAAc                                                                                                       <  2:71235/148‑1 (MQ=255)                                                  tACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCGAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGcatcat                                                                                                >  1:257150/1‑149 (MQ=255)                                                    ccAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTg                                                                                              >  3:180888/1‑149 (MQ=255)                                                    ccAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTg                                                                                              <  6:152986/149‑1 (MQ=255)                                                      aaGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTgg                                                                                             >  6:53850/1‑148 (MQ=255)                                                                     tCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAAAGCATCATTGGTCAGCTTGGCTtgctg                                                                             >  1:239870/1‑149 (MQ=255)                                                                       aGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTtgctgc                                                                            <  5:53850/148‑1 (MQ=255)                                                                       aGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCa                                                                           >  1:2199/1‑149 (MQ=255)                                                                              gcgcTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCAtt                                                                    >  7:135560/1‑149 (MQ=255)                                                                                gcTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTa                                                                  >  3:274473/1‑149 (MQ=255)                                                                                  tGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACg                                                                >  7:152796/1‑149 (MQ=255)                                                                                  tGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTAAACGCATCATTGGTCAGCTTGGCTTGCTGGACGCCATTTACg                                                                >  4:236801/1‑149 (MQ=255)                                                                                      gAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGATGCACGCCATTTACGTa                                                              >  7:228024/1‑147 (MQ=255)                                                                                         gCCCATCGGCGAGCGCATCAAGCAACGCGGCGAAGTAGGCGTTGCAGTGGAAGTGAC‑CACGCCACCCATAACAAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGACAGCTTGGCCTGCTCCACGCCATTTACGTACAGCg                                                         >  7:196094/1‑149 (MQ=255)                                                                                            cATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGACGCAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGt                                                        >  1:85107/1‑148 (MQ=255)                                                                                            cATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCa                                                      >  8:81027/1‑149 (MQ=255)                                                                                             aTCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCAt                                                     >  4:256654/1‑149 (MQ=255)                                                                                                   gAGCGCATCAAGCATGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCATCGCATCTTTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTacc                                               <  8:75797/149‑1 (MQ=255)                                                                                                   gAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTacc                                               <  8:141998/149‑1 (MQ=255)                                                                                                         aTCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACcgc                                         <  5:222289/149‑1 (MQ=255)                                                                                                              gcaagGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTCGGGCGAGGCGTAGCGTCTTCATCGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGCTACCGACCGCGCGGt                                    >  1:258735/1‑149 (MQ=255)                                                                                                                aagGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCg                                  >  2:27884/1‑149 (MQ=255)                                                                                                                 agGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTa                                                                                                                  >  8:263846/1‑68 (MQ=255)                                                                                                                 agGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTa                                                                                                                  <  7:263846/68‑1 (MQ=255)                                                                                                                 agGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCgg                                 <  5:77255/149‑2 (MQ=255)                                                                                                                      gcgAAGTAGGCGTTGCAGTGGTAGTGAC‑CATGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTg                            <  7:195832/149‑1 (MQ=255)                                                                                                                      gcgAAGTAGGCGTTGCAGTGGTAGTGAC‑CATGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTg                            <  7:22472/149‑1 (MQ=255)                                                                                                                       cgAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGa                           <  1:309066/149‑1 (MQ=255)                                                                                                                       cgAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGa                           <  7:170882/149‑1 (MQ=255)                                                                                                                        gAAGTAGGGGTTGGAGTGGTAGTGAC‑CTGGGCAGCCATAAGGAGTCTTGGGGAGAGGGTTAGAGTCTTCAACGCATCATTGGTCTGGTTGGCGTGCTGCTCGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAg                          <  8:196094/149‑1 (MQ=255)                                                                                                                         aaGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAgg                         <  4:274473/149‑1 (MQ=255)                                                                                                                          aGTAGGCGGTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCGGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCCACCGCGCGGTCGTATTTGAGGa                        <  7:49088/149‑1 (MQ=255)                                                                                                                             aGGCGTTGCAGTGGTAGTGAC‑CAGGCCCCCCATAACCAGGACTGGTGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATc                      >  8:162966/1‑148 (MQ=255)                                                                                                                                cGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAgc                  >  6:48184/1‑149 (MQ=255)                                                                                                                                 gTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAgcg                 >  4:69300/1‑149 (MQ=255)                                                                                                                                    gCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCa               <  2:257150/148‑1 (MQ=255)                                                                                                                                    gCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCa               <  1:125577/148‑1 (MQ=255)                                                                                                                                            aGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGATGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCAtgcg       <  4:119400/148‑1 (MQ=255)                                                                                                                                              tGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCAtgcggtg    <  2:2199/149‑1 (MQ=255)                                                                                                                                               gAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTTCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCAtgcggtg    <  7:101665/148‑1 (MQ=255)                                                                                                                                                aC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTTGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCAtgcggtgc   <  7:180872/148‑1 (MQ=255)                                                                                                                                                aC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCAtgcggtgcg  <  6:329291/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CAACCGGCTGACATCTTGTTTAGTAGAGAGTTTTTGACAAGTAGCGGTCTACCAAGGCCAGTTCGCTGTCAGACAGCGCGCTGTTGAAGCCCATCGGCGAGCGCATCAAGCAAGGCGGCGAAGTAGGCGTTGCAGTGGTAGTGAC‑CAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCG  >  NC_002947/4379742‑4380033

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