Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,949	A→G	5.0%	Q60Q (CAA→CAG)	PP_3574 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,054,949	0	A	G	5.0%	296.1 / 4.0	99	Q60Q (CAA→CAG)	PP_3574	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (48/46);  new base G (3/2);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ACCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGCC  >  NC_002947/4054803‑4055096                                                                                                                                                   |                                                                                                                                                      aCCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACt                                                                                                                                                     <  8:59703/149‑1 (MQ=255)   ccTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGc                                                                                                                                                   >  8:315183/1‑149 (MQ=255)   ccTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGc                                                                                                                                                   <  3:72776/149‑1 (MQ=255)     tGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCtt                                                                                                                                                 >  5:164811/1‑149 (MQ=255)      ggACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGCCCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCtt                                                                                                                                                 >  4:260823/1‑148 (MQ=255)       gACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCtt                                                                                                                                                 <  7:258224/147‑1 (MQ=255)       gACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTg‑‑a                                                                                                                                             <  7:155230/149‑1 (MQ=255)          aaCGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAGTTGCTTg‑‑acgt                                                                                                                                          <  5:294976/149‑4 (MQ=255)          aaCGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTg‑‑aag                                                                                                                                           >  1:399279/1‑148 (MQ=255)            cGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑aagaag                                                                                                                                        >  1:144576/1‑149 (MQ=255)             gACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGc                                                                                                                                       <  3:424465/149‑1 (MQ=255)             gACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGc                                                                                                                                       <  4:202416/149‑1 (MQ=255)               cTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCg                                                                                                                                     >  2:187756/1‑149 (MQ=255)                  gtgtCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCt                                                                                                                                  >  2:116199/1‑149 (MQ=255)                   tgtCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTc                                                                                                                                 >  2:359568/1‑149 (MQ=255)                     tCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCgg                                                                                                                               <  3:203539/149‑1 (MQ=255)                      cAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGa                                                                                                                              >  1:154923/1‑149 (MQ=255)                       aGGCCGTGCGGGCGGGCAATACCCTCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGAc                                                                                                                             <  8:216987/149‑1 (MQ=255)                       aGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGAc                                                                                                                             >  5:195545/1‑149 (MQ=255)                            gTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTc                                                                                                                        <  3:402206/149‑1 (MQ=255)                             tGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCg                                                                                                                       <  6:195545/149‑1 (MQ=255)                              gcgggcggGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGcc                                                                                                                                      <  4:369026/133‑1 (MQ=255)                              gcgggcggGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGcc                                                                                                                                      >  3:369026/1‑133 (MQ=255)                               cgggcgggCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGca                                                                                                                     >  3:271551/1‑149 (MQ=255)                                gggcgggcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAATCCGTCTCGGACCTGTCGca                                                                                                                     >  4:292861/1‑148 (MQ=255)                                gggcgggcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGca                                                                                                                     >  6:45693/1‑148 (MQ=255)                                    gggcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCg                                                                                                                >  1:417692/1‑149 (MQ=255)                                      gcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGt                                                                                                               >  3:210921/1‑148 (MQ=255)                                      gcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTc                                                                                                              <  2:7030/149‑1 (MQ=255)                                      gcAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTc                                                                                                              <  1:187756/149‑1 (MQ=255)                                       cAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTc                                                                                                              >  5:169557/1‑148 (MQ=255)                                       cAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTc                                                                                                              >  3:331148/1‑148 (MQ=255)                                        aaTACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCaa                                                                                                            <  2:399279/149‑1 (MQ=255)                                             cATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAAAAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGacca                                                                                                       >  2:366805/1‑149 (MQ=255)                                                  aCGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccaccacc                                                                                                  >  7:258932/1‑149 (MQ=255)                                                   cGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccaccac                                                                                                   <  6:164811/147‑1 (MQ=255)                                                     tGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccaccac                                                                                                   <  1:385976/145‑1 (MQ=255)                                                     tGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGaccaccac                                                                                                   >  2:385976/1‑145 (MQ=255)                                                     tGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTAc                                                                                               <  1:296068/149‑1 (MQ=255)                                                     tGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTAc                                                                                               <  3:260823/149‑1 (MQ=255)                                                      gcgcGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGAACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTAc                                                                                               >  7:287252/1‑148 (MQ=255)                                                            ccAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCAcc                                                                                         >  7:57786/1‑148 (MQ=255)                                                                 tCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCcg                                                                                   <  8:57786/149‑1 (MQ=255)                                                                      aCCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCtt                                                                              >  6:392365/1‑149 (MQ=255)                                                                       ccGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCtt                                                                              <  5:392365/148‑1 (MQ=255)                                                                        cGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTcc                                                                            >  6:363465/1‑149 (MQ=255)                                                                          aCTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGCCTCGCACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCCCTTCcgc                                                                          >  8:177943/1‑149 (MQ=255)                                                                            ttCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCcgc                                                                          >  7:185833/1‑147 (MQ=255)                                                                            ttCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCcgc                                                                          >  6:173933/1‑147 (MQ=255)                                                                              cGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCcgc                                                                                  <  1:237297/137‑1 (MQ=255)                                                                              cGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCcgc                                                                                  >  2:237297/1‑137 (MQ=255)                                                                              cGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGc                                                                      >  6:199760/1‑149 (MQ=255)                                                                                aaGGCAAGCTGGTGGGCCTGTGTGACCCGCGCGCACAAACCGAACAGGCGATGAATAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACACCACCGACCCGCGCTTCCGCGAGCCg                                                                    <  3:292861/149‑1 (MQ=255)                                                                                aaGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCg                                                                    <  6:169557/149‑1 (MQ=255)                                                                                  ggCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGgt                                                                  <  5:199760/149‑1 (MQ=255)                                                                                         tggtggGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAag                                                           >  7:99311/1‑149 (MQ=255)                                                                                         tggtggGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAag                                                           >  5:159710/1‑149 (MQ=255)                                                                                                ccTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTgg                                                    <  7:301668/149‑1 (MQ=255)                                                                                                 cTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTggg                                                   <  2:144576/149‑1 (MQ=255)                                                                                                 cTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAACAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTggg                                                   >  1:387628/1‑149 (MQ=255)                                                                                                    ggTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGc                                                  <  8:41837/147‑1 (MQ=255)                                                                                                     gttccccTCGCGAACAAACAGAACAGGTGATGAAGAACGTCAAGCAACT‑CTTGCTAAGAAGCGGTCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCaa                                                <  1:359568/145‑1 (MQ=255)                                                                                                            gcgcgcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTc                                         >  5:26568/1‑148 (MQ=255)                                                                                                              gcgcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAg                                      <  4:210921/149‑1 (MQ=255)                                                                                                               cgcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAgg                                     <  6:159710/149‑1 (MQ=255)                                                                                                                gcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAggg                                    <  1:116199/149‑1 (MQ=255)                                                                                                                gcACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAggg                                    >  8:164751/1‑149 (MQ=255)                                                                                                                   caAACCGAACAGGCTATGAAGAACTTCAAGCTACTGCTTG‑‑AAGAAGGCGGCTGGGTCCTGTCGCACATCGTCTAGACCACCACGTTCATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCgt                                 <  7:177943/149‑1 (MQ=255)                                                                                                                    aaaCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTc                                                                             <  8:39239/104‑1 (MQ=255)                                                                                                                    aaaCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTc                                                                             >  7:39239/1‑104 (MQ=255)                                                                                                                    aaaCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCgtg                                >  4:233636/1‑149 (MQ=255)                                                                                                                    aaaCCGAACAGGCGATGAAGAACGTCAAGAAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGAACATCGTCAAGACAACCACATACATAACCGACCCTCGCTTCCGAGAGCCGGTGTACAAGGAGGTGGGCAAGTGGATCAAGGACgtg                                >  4:239123/1‑149 (MQ=255)                                                                                                                      aCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGgt                                                                  <  3:264539/113‑1 (MQ=255)                                                                                                                      aCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGgt                                                                  >  4:264539/1‑113 (MQ=255)                                                                                                                      aCCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGtt                              <  5:363465/149‑1 (MQ=255)                                                                                                                        cGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTcc                            >  2:140357/1‑149 (MQ=255)                                                                                                                         gAACAGGCGAGTAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTccc                           <  8:307966/149‑1 (MQ=255)                                                                                                                              ggCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGAttt                      <  4:331148/149‑1 (MQ=255)                                                                                                                               gCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTc                     >  5:92239/1‑149 (MQ=255)                                                                                                                               gCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTc                     <  2:154923/149‑1 (MQ=255)                                                                                                                               gCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGACGGTGTACAAGGAGGTGGGAAAGTGGCTCAAGGGCGTGTTCCCGATTTc                     <  1:366805/149‑1 (MQ=255)                                                                                                                                   tGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCAcc                 >  3:166895/1‑149 (MQ=255)                                                                                                                                     aagaaCGTCAAGCAACTTCTTG‑‑AAGAAGACGGATCTGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCTCGAGCCGGTGTACAATGATGTGGGCAAGTTGCTCAAGGGCGTGTTCCCGATTTCCACCg                <  3:239123/148‑1 (MQ=255)                                                                                                                                     aagaaCGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCg                <  3:233636/148‑1 (MQ=255)                                                                                                                                       gaaCGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGt             >  3:145715/1‑149 (MQ=255)                                                                                                                                          cttCAAGCAGTTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtgg          >  4:244273/3‑149 (MQ=255)                                                                                                                                          cGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtgg          <  8:99311/149‑1 (MQ=255)                                                                                                                                          cGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtgg          <  7:164751/149‑1 (MQ=255)                                                                                                                                           ttCAAGCAGTTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggt         >  6:303995/2‑149 (MQ=255)                                                                                                                                           ttCAAGCAGTTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggt         >  2:257854/2‑149 (MQ=255)                                                                                                                                             cAAGCAGTTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtg        <  5:303995/148‑1 (MQ=255)                                                                                                                                             cAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtgg       >  3:116012/1‑149 (MQ=255)                                                                                                                                              aaGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGt                                              <  7:284349/109‑1 (MQ=255)                                                                                                                                              aaGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtggt      >  6:28707/1‑149 (MQ=255)                                                                                                                                              aaGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTAGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGATTCCGCGATCCGGTGGACAAGGACGTGGGCAAGt                                              >  8:284349/1‑109 (MQ=255)                                                                                                                                              aaGCAACTGCGTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACCTCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtgg       <  4:271551/148‑1 (MQ=255)                                                                                                                                               aGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtggtg     >  7:402925/1‑149 (MQ=255)                                                                                                                                               aGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTtggtggtg     <  5:173933/149‑1 (MQ=255)                                                                                                                                                gCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAg                                      >  2:317767/1‑115 (MQ=255)                                                                                                                                                  aaCTGCTTG‑‑AAGAAGCCGGCTCGGACCTGGCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGcc  <  2:369359/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ACCCTGGACAACGACTTGTGTCAGGCCGTGCGGGCGGGCAATACCATCTACGTGCGCGGCCAGGTCGGTACCGACTTCGAAGGCAAGCTGGTGGGCCTGGGTGACCCGCGCGCACAAACCGAACAGGCGATGAAGAACGTCAAGCAACTGCTTG‑‑AAGAAGCCGGCTCGGACCTGTCGCACATCGTCAAGACCACCACCTACATCACCGACCCGCGCTTCCGCGAGCCGGTGTACAAGGAGGTGGGCAAGTGGCTCAAGGGCGTGTTCCCGATTTCCACCGGGTTGGTGGTGGCC  >  NC_002947/4054803‑4055096

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