Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,970,525	A→G	6.0%	W257R (TGG→CGG)	PP_5236 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,970,525	0	A	G	6.0%	244.4 / 4.5	84	W257R (TGG→CGG)	PP_5236	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/34);  new base G (3/2);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

TCGGCGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGG  >  NC_002947/5970377‑5970669                                                                                                                                                      |                                                                                                                                                 tCGGCGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCa                                                                                                                                                  <  7:995/149‑1 (MQ=255)     cGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCgg                                                                                                                                              <  3:306935/149‑1 (MQ=255)     cGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCgg                                                                                                                                              >  5:131335/1‑149 (MQ=255)            cTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGgcagca                                                                                                                                       <  3:113841/149‑1 (MQ=255)            cTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGgcag                                                                                                                                         >  8:32130/1‑147 (MQ=255)            cTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAACAAGGCTGCCGCCCCAGCGGGgcag                                                                                                                                         >  1:126086/1‑147 (MQ=255)               cACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCTGCCCCGCTGGGGcggcagcc                                                                                                                                    <  5:31310/149‑3 (MQ=255)                 ccTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTac                                                                                                                                  <  5:336287/149‑1 (MQ=255)                  cTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTaca                                                                                                                                 >  3:298753/1‑149 (MQ=255)                  cTTGTCAGCCTGTGGCAGCTTCAGCCGGCCATCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTaca                                                                                                                                 <  8:225138/149‑1 (MQ=255)                   ttGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGcagcag                                                                                                                                      <  8:193217/143‑1 (MQ=255)                   ttGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGcagcag                                                                                                                                      >  7:193217/1‑143 (MQ=255)                      tCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGTGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAg                                                                                                                              >  1:285569/1‑148 (MQ=255)                      tCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCTGCCCCGCTGGGGcggcagcc                                                                                                                                    <  5:8863/142‑3 (MQ=255)                      tCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCTGCCCCGCTGGGGcggcagcc                                                                                                                                    >  6:8863/1‑140 (MQ=255)                                gCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGGTTGCGCACACCGGGCCTTTGACAAACACCTCGAAGGCAACAAGCCTGCCGCCCCCCCGGGGCAGAAGGTACACAGCAAGGGCCTg                                                                                                                    >  1:110586/1‑149 (MQ=255)                                gCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTg                                                                                                                    >  2:185809/1‑148 (MQ=255)                                gCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTg                                                                                                                    >  2:273285/1‑148 (MQ=255)                                  aGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGAt                                                                                                                 <  4:86655/149‑1 (MQ=255)                                      tCAGCCGTCCTTCGACATCAATCCATCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCTAACGCATCAAGTCTGCCGCCCCAGCTGGGCAGCAGGTACACAGCAAGGGc                                                                                                                       <  2:1264/139‑1 (MQ=255)                                      tCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCACCAGCGGGGCAGCAGGTACACAGCAAGGGc                                                                                                                       >  1:1264/1‑139 (MQ=255)                                            ggCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAag                                                                                                       >  8:97681/1‑149 (MQ=255)                                                ttCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGcagcag                                                                                                                                      >  8:211708/1‑114 (MQ=255)                                                ttCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGcagcag                                                                                                                                      <  7:211708/114‑1 (MQ=255)                                                   gACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCg                                                                                                 <  1:273285/148‑1 (MQ=255)                                                      aGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggc                                                                                             >  1:223218/1‑149 (MQ=255)                                                      aGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGAGCGCTTCGAAGCAGCCGgggc                                                                                             >  8:312954/1‑149 (MQ=255)                                                        ccATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcgg                                                                                           >  4:297188/1‑149 (MQ=255)                                                              agcgcagcgcATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGAt                                                                                      >  7:251708/1‑148 (MQ=255)                                                               gcgcagcgcaTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCa                                                                                    <  3:234374/149‑1 (MQ=255)                                                               gcgcagcgcaTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCa                                                                                    <  8:251708/149‑1 (MQ=255)                                                                cgcagcgcaTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAg                                                                                   <  6:131335/149‑1 (MQ=255)                                                                 acagcgcaTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAgc                                                                                  >  7:122825/2‑149 (MQ=255)                                                                     cgcaTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCAt                                                                              <  3:297188/149‑1 (MQ=255)                                                                        aTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGgc                                                                           >  4:97939/1‑149 (MQ=255)                                                                        aTTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGgc                                                                           >  7:90742/1‑149 (MQ=255)                                                                             gCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCTCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGc                                                                      <  7:242906/149‑1 (MQ=255)                                                                             gCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGc                                                                      <  4:305740/149‑1 (MQ=255)                                                                             gCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGc                                                                      <  7:97681/149‑1 (MQ=255)                                                                             gCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGc                                                                      <  4:287220/149‑1 (MQ=255)                                                                                   gCTCGTCCCA‑GGTTGTGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTaa                                                                <  2:285569/149‑1 (MQ=255)                                                                                       gTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGc                                                            >  7:86102/1‑149 (MQ=255)                                                                                         cccA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCG‑CCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCtt                                                                                                            <  2:21790/98‑1 (MQ=255)                                                                                         cccA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCG‑CCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCtt                                                                                                            >  1:21790/1‑98 (MQ=255)                                                                                          ccA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAg                                                         >  1:305118/1‑149 (MQ=255)                                                                                          ccA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAg                                                         >  5:266802/1‑149 (MQ=255)                                                                                               gTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGaa                                                                                                        >  4:106035/1‑98 (MQ=255)                                                                                               gTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGaa                                                                                                        <  3:106035/98‑1 (MQ=255)                                                                                                    gcACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGt                                                >  4:100786/1‑149 (MQ=255)                                                                                                        aGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTc                                            <  7:32130/149‑1 (MQ=255)                                                                                                            ggCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATgg                                        <  2:223218/149‑1 (MQ=255)                                                                                                                       aaCAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCtt                             >  4:19350/1‑149 (MQ=255)                                                                                                                        aCAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGc                                                            >  1:220550/1‑117 (MQ=255)                                                                                                                        aCAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGc                                                            <  2:220550/117‑1 (MQ=255)                                                                                                                        aCAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTg                            >  6:202896/1‑149 (MQ=255)                                                                                                                         cAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGc                           <  5:202896/149‑1 (MQ=255)                                                                                                                            cTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGcc                          >  7:319316/1‑147 (MQ=255)                                                                                                                            cTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCa                        >  7:319321/1‑149 (MQ=255)                                                                                                                            cTCCAACCCATCAAGGCTGCCGCCCCAGCGCGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTACGAAGCAGCCGGCGCGGGCGATCAGCCCATCGCCCACCCAGTAAAAGCGAGCGTACTCGTTTTCATGGTCGACCCGCATGccc                         >  6:120557/1‑148 (MQ=255)                                                                                                                             tCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGcc                          >  7:250913/1‑146 (MQ=255)                                                                                                                             tCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGcc                          >  3:211499/1‑146 (MQ=255)                                                                                                                             tCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGcc                          >  1:216457/1‑146 (MQ=255)                                                                                                                               gAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGa                      <  2:305118/148‑1 (MQ=255)                                                                                                                               gAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGa                      <  8:122825/148‑1 (MQ=255)                                                                                                                                   gCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAgcg                 <  3:97939/149‑1 (MQ=255)                                                                                                                                   gCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTACGCAGCAGCCGGGTCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGTGTCCTCGGTTTCATCGTCGACCCGCTTGCCCAGATAgcg                 >  2:66067/1‑149 (MQ=255)                                                                                                                                    cATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAgcgc                <  2:110586/149‑1 (MQ=255)                                                                                                                                     aTCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATa                    >  6:48623/1‑144 (MQ=255)                                                                                                                                     aTCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATa                    <  5:48623/144‑1 (MQ=255)                                                                                                                                     aTCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGcgcg               >  4:283636/1‑149 (MQ=255)                                                                                                                                     aTCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGcgcg               >  2:256755/1‑149 (MQ=255)                                                                                                                                        aaGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCgg              >  8:147331/1‑147 (MQ=255)                                                                                                                                        aaGGCTGCCGCCCCAGAGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCAGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCgg              >  2:206082/1‑147 (MQ=255)                                                                                                                                         aGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCAt                                                                              >  2:1621/1‑82 (MQ=255)                                                                                                                                          ggCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGc          <  2:126086/149‑1 (MQ=255)                                                                                                                                            ctgctgCCCCGCTGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGa        >  1:134038/1‑149 (MQ=255)                                                                                                                                            ctgctgCCCCGCTGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGa        >  1:344762/1‑149 (MQ=255)                                                                                                                                              gccgccCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACt      >  3:249818/1‑149 (MQ=255)                                                                                                                                               cccccccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATGGCTTCGAAGCAGCCGGGGCGGGCGATCTGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtg     <  5:120557/146‑1 (MQ=255)                                                                                                                                                 gccCCAGCGGGGCAGCAGGGACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtgt    <  1:317092/148‑1 (MQ=255)                                                                                                                                                  ccccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGGCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTgg  <  8:90742/149‑1 (MQ=255)                                                                                                                                                  ccccAGCGGGGCAGCAGGGACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtgt    >  2:317092/1‑147 (MQ=255)                                                                                                                                                    ccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCtt                             <  1:257625/120‑1 (MQ=255)                                                                                                                                                    ccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCtt                             >  2:257625/1‑120 (MQ=255)                                                                                                                                                    ccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTgg  >  5:156090/1‑147 (MQ=255)                                                                                                                                                      |                                                                                                                                                 TCGGCGACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCA‑GGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGG  >  NC_002947/5970377‑5970669

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