Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,639,342	G→A	7.2%	A119T (GCC→ACC)	PP_3206 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,639,342	0	G	A	7.2%	200.0 / 4.3	69	A119T (GCC→ACC)	PP_3206	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base A (2/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

TGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCG  >  NC_002947/3639197‑3639490                                                                                                                                                  |                                                                                                                                                      tgtCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGc                                                                                                                                                      >  6:145367/1‑147 (MQ=255) tgtCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGc                                                                                                                                                      >  8:186092/1‑147 (MQ=255) tgtCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGCGTACACca                                                                                                                                                    <  8:117316/149‑1 (MQ=255)     gATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCcag‑cag                                                                                                                                               >  7:278999/1‑149 (MQ=255)          ggTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTcg                                                                                                                                          >  6:203650/1‑149 (MQ=255)           gTCTTCCATCTCGTCAGTGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTcg                                                                                                                                          <  5:97923/148‑1 (MQ=255)           gTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTcgc                                                                                                                                         >  8:332549/1‑149 (MQ=255)                ccATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTAACCAGGTAAACCTGCAAGCAAGCCTGGAGATGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCCgtg                                                                                                                                     >  5:219965/1‑148 (MQ=255)                   tctcGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTAcg                                                                                                                                 >  5:182186/1‑149 (MQ=255)                    ctcGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTAcgg                                                                                                                                >  3:165818/1‑149 (MQ=255)                     tcGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTAcggc                                                                                                                               >  8:322206/1‑149 (MQ=255)                       gTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCGTCTGTCCACAGTTGCTGGTGTACGCCAG‑CAGCGTCGCTGTGTAcggcgg                                                                                                                             <  8:354872/149‑1 (MQ=255)                       gTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAAGCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTAcggagg                                                                                                                             >  8:127009/1‑149 (MQ=255)                        tCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTGACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACggcggc                                                                                                                            <  7:75088/149‑1 (MQ=255)                          aGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCaa                                                                                                                          <  8:208829/149‑1 (MQ=255)                             gTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGtt                                                                                                                       <  5:145367/149‑1 (MQ=255)                                ccTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGcc                                                                                                                    >  5:36734/1‑149 (MQ=255)                                ccTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGc                                                                                                                     <  1:158150/148‑1 (MQ=255)                                 cTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCg                                                                                                                   >  8:16672/1‑149 (MQ=255)                                  tGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGa                                                                                                                  >  1:74090/1‑149 (MQ=255)                                   ggcggcGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGAc                                                                                                                 >  5:111163/1‑149 (MQ=255)                                       gcgcTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTTGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACgcgc                                                                                                             <  2:74090/149‑1 (MQ=255)                                          cTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATg                                                                                                          <  3:94013/149‑1 (MQ=255)                                           tGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGTTGTACGGCGGCAAGTTGCCGACGCGCATgg                                                                                                         >  3:302147/1‑149 (MQ=255)                                             ccGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATgg                                                                                                         <  4:165818/147‑1 (MQ=255)                                                 gaCGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAg                                                                                                   >  2:132460/1‑149 (MQ=255)                                                   cggattACGAACGGGGGTTCCATGTAAACCTGGAAGGAATGCTGGATCTGCTCGTACAACTTCGTAATCCCCTCTTTCCACCGGTGGTGGTGTTAGCCAGTCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAg                                                                                                   <  1:273360/143‑1 (MQ=255)                                                    gCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgc                                                                                                 >  6:321835/1‑148 (MQ=255)                                                           gAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCcagccagc                                                                                         >  8:109686/1‑149 (MQ=255)                                                             aCGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGt                                                                                       >  2:69936/1‑149 (MQ=255)                                                                          gTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTc                                                                          <  7:332549/149‑1 (MQ=255)                                                                          gTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTc                                                                          <  7:186092/149‑1 (MQ=255)                                                                             aaCCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTa                                                                       >  1:2688/1‑149 (MQ=255)                                                                              aCCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTAt                                                                      >  4:111734/1‑149 (MQ=255)                                                                               ccTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATg                                                                     >  4:320117/1‑149 (MQ=255)                                                                                cTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgc                                                                    >  8:112932/1‑149 (MQ=255)                                                                                      gcaagcCGTGAGCTGCTCGACCAACTGCGTAATCCCCTCTGGCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGca                                                              <  6:125299/149‑1 (MQ=255)                                                                                          gcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACaaa                                                          >  8:153145/1‑149 (MQ=255)                                                                                           ccTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAAc                                                         <  2:2688/149‑1 (MQ=255)                                                                                           ccTGGAGCTGCTCGACCAACTGCGTAATCACCTCTGTCCACCGGTGCTGGTCTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAACTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAAc                                                         >  5:62666/1‑149 (MQ=255)                                                                                               gAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAt                                                     <  4:333539/149‑1 (MQ=255)                                                                                                aGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtg                                                    <  2:250333/149‑1 (MQ=255)                                                                                                     ctCGACCAACTGCTTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGTCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTGCTATGCGGCGCACAAACGCAttgtg                                                 <  6:219965/147‑1 (MQ=255)                                                                                                     ctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGga                                               <  8:91695/149‑1 (MQ=255)                                                                                                            aaCTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcgcg                                        >  4:39822/1‑149 (MQ=255)                                                                                                               tGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTg                                     >  4:256949/1‑149 (MQ=255)                                                                                                                  gTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAg                                  <  7:16672/149‑1 (MQ=255)                                                                                                                    aaTCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAgg                                 <  6:36734/148‑1 (MQ=255)                                                                                                                        cccTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCtt                            >  8:234705/1‑149 (MQ=255)                                                                                                                         ccTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTg                           >  1:244730/1‑149 (MQ=255)                                                                                                                          ctctGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACACACGCATGGTGGAGATCGCGTTGCAGGATCTTGc                          >  2:238318/1‑149 (MQ=255)                                                                                                                          ctctGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGc                          >  1:160852/1‑149 (MQ=255)                                                                                                                             tGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCTTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCccg                       >  1:13686/1‑149 (MQ=255)                                                                                                                                 cACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCccgccg                    <  6:111163/148‑1 (MQ=255)                                                                                                                                 cACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGt                   >  2:30421/1‑149 (MQ=255)                                                                                                                                  gacgctGCTGGCGTACACCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATc                              <  1:349822/132‑1 (MQ=255)                                                                                                                                  gacgctGCTGGCGTACACCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATc                              >  2:349822/6‑137 (MQ=255)                                                                                                                                   ccGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTgg                 <  3:39822/149‑1 (MQ=255)                                                                                                                                     gctGCTGGCGTACACCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTg                  >  4:22716/3‑146 (MQ=255)                                                                                                                                     gctGCTGGCGTACACCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTg                  <  3:22716/144‑1 (MQ=255)                                                                                                                                      gTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGa              <  7:127009/149‑1 (MQ=255)                                                                                                                                       tGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtgg                                                   >  7:227861/1‑111 (MQ=255)                                                                                                                                       tGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtgg                                                   <  8:227861/111‑1 (MQ=255)                                                                                                                                       tGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAg             >  4:72335/1‑149 (MQ=255)                                                                                                                                        gCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAgg            <  7:153145/149‑1 (MQ=255)                                                                                                                                          tGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGtt          >  6:123943/1‑149 (MQ=255)                                                                                                                                             tgtACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGAc       <  7:109686/149‑1 (MQ=255)                                                                                                                                             tgtACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGAc       <  5:203650/149‑1 (MQ=255)                                                                                                                                               tACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACgg     <  3:256949/149‑1 (MQ=255)                                                                                                                                                  gCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCcg  >  6:333524/1‑149 (MQ=255)                                                                                                                                                  gCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCcg  >  5:284472/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      TGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAG‑CAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCG  >  NC_002947/3639197‑3639490

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