Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,921,592	A→C	10.3%	P667P (CCA→CCC)	PP_4330 →	toxin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,921,592	0	A	C	10.3%	180.1 / 12.1	68	P667P (CCA→CCC)	PP_4330	toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base C (3/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGA  >  NC_002947/4921446‑4921740                                                                                                                                                   |                                                                                                                                                     aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                                    >  5:192076/1‑149 (MQ=255) aaaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAgg                                                                                                                                                    >  1:49612/1‑149 (MQ=255)  aaGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAggg                                                                                                                                                   <  7:192960/149‑1 (MQ=255)     tCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTg                                                                                                                                                <  1:9142/149‑1 (MQ=255)       aTTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGc                                                                                                                                               <  6:23639/148‑1 (MQ=255)        ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                             <  5:85793/149‑1 (MQ=255)        ttAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAg                                                                                                                                             <  1:22321/149‑1 (MQ=255)           gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGAGAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTTCGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGt                                                                                                                                           <  4:3297/148‑1 (MQ=255)           gCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCACTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTa                                                                                                                                          <  4:98641/149‑1 (MQ=255)            cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaata                                                                                                                                          <  3:195491/148‑6 (MQ=255)            cGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCCGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCCCGGGCTGccaa                                                                                                                                            >  4:195491/1‑143 (MQ=255)               gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGGTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGc                                                                                                                                      >  4:43703/1‑149 (MQ=255)               gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCTCGGGCTGc                                                                                                                                               <  1:76740/140‑1 (MQ=255)               gCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCTGGGGCTCGGGCTGc                                                                                                                                               >  2:76740/1‑140 (MQ=255)                  aTTGCTGGCCAACCTTGCCGCCAACGCCGCGCGACTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTAGGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGc                                                                                                                                   <  2:102901/149‑1 (MQ=255)                   ttGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCt                                                                                                                                  <  1:65860/149‑1 (MQ=255)                     gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                                 <  3:167908/148‑1 (MQ=255)                     gCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCtt                                                                                                                                 >  2:153752/1‑148 (MQ=255)                              ccTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAg                                                                                                                       <  7:176905/149‑1 (MQ=255)                                    cgccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCACTTGAAGAGCAAGGcttct                                                                                                                 <  6:192076/149‑1 (MQ=255)                                      ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                                >  2:172387/1‑148 (MQ=255)                                      ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCttcttc                                                                                                               >  3:270315/1‑149 (MQ=255)                                      ccAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGCTATTGCCGCTTGAAGAGCAAGGcttctt                                                                                                                <  1:172387/148‑1 (MQ=255)                                         aggccccGCGCCTGCCCTACAGTGAACACCGCATGCGCCGGTTTCTGTCGAAGGCTGCGCCTACGGCTTACGACCAGACCCCAAGGCTATTGGCAGCCCGGGCCCCTGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCgg                                                                                                            <  6:184287/147‑1 (MQ=255)                                           gccgcGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAaga                                                                                                                            >  2:8728/1‑131 (MQ=255)                                           gccgcGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAaga                                                                                                                            <  1:8728/131‑1 (MQ=255)                                                gcgcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCt                                                                                                     <  7:32973/149‑1 (MQ=255)                                                  gcCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGc                                                                                                   <  5:277941/149‑1 (MQ=255)                                                   ccTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGcc                                                                                                  <  1:153752/149‑1 (MQ=255)                                                           aCAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                          <  4:270315/149‑1 (MQ=255)                                                             aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGTTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTa                                                                                        <  1:213735/149‑1 (MQ=255)                                                             aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                          <  5:86433/147‑1 (MQ=255)                                                             aGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGc                                                                                          >  6:86433/1‑147 (MQ=255)                                                               tGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTAt                                                                                       >  5:132525/1‑148 (MQ=255)                                                                   tACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGC                                                                                  >  3:91491‑M1/1‑149 (MQ=255)                                                                        cATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGCTCCGCTATTACGCGCAG                                                                              >  3:185551‑M1/1‑148 (MQ=255)                                                                        cATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAG                                                                              >  7:253410‑M1/1‑148 (MQ=255)                                                                           gcgcCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCGAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTAc                                                                           >  7:258576/1‑148 (MQ=255)                                                                                tGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTCGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTCAAGAGCAAGGCTTCTTCCGGCCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCa                                                                     >  5:30639/1‑149 (MQ=255)                                                                                           aaGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATgg                                                          >  3:76288/1‑149 (MQ=255)                                                                                               cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaatt                                                      >  8:87097‑M1/1‑125 (MQ=255)                                                                                               cTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAAc                                                      >  2:167260/1‑149 (MQ=255)                                                                                                ttCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTCCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACg                                                     >  2:141710/1‑149 (MQ=255)                                                                                                  tttttACCTCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCATGGTTGCAGGTATTGCCGTTTGAAGAGCAAGGCTTCTTCCGTGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGt                                                   <  7:111290/145‑1 (MQ=255)                                                                                                   gCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTg                                                  >  3:241799/1‑149 (MQ=255)                                                                                                      tACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCa                                               >  3:298453/1‑149 (MQ=255)                                                                                                       aCCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCATCGCTGTATGGCAACTGttata                                               >  8:111290/1‑144 (MQ=255)                                                                                                          gCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGc                                           <  1:167260/149‑1 (MQ=255)                                                                                                               cGAACAGTACGCAAGGCTATTGGGAGCCCGGTCGCCAGGGCTGCAGGTATTGGCGCTTTAAGAGCATGGCTTCTTGCGGGCCTTCTGCCAGAGCCGCTATTACGCGGAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATGGccacc                                       <  6:30639/148‑1 (MQ=255)                                                                                                                gAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCGATGccaccgt                                     <  2:49612/149‑3 (MQ=255)                                                                                                                    aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCG                                                                                 >  5:259281‑M1/1‑101 (MQ=255)                                                                                                                    aGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCG                                                                                 <  6:259281‑M1/101‑1 (MQ=255)                                                                                                                           aaGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTg                          >  4:210179/1‑149 (MQ=255)                                                                                                                             ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgc                        >  5:258070‑M1/1‑95 (MQ=255)                                                                                                                             ggCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGc                         >  7:11899/1‑148 (MQ=255)                                                                                                                               atacctgCAGCCCTGGGGCCCGGGCTGCCGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                           <  1:126748/138‑1 (MQ=255)                                                                                                                               atacctgCAGCCCTGGGGCCCGGGCTGCCGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCt                           >  2:126748/7‑144 (MQ=255)                                                                                                                                    ggCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAt                  <  4:241799/148‑1 (MQ=255)                                                                                                                                     gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt                >  3:248353/1‑149 (MQ=255)                                                                                                                                     gCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgt                >  8:11911/1‑149 (MQ=255)                                                                                                                                      cAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGAtgtg               <  4:248353/149‑1 (MQ=255)                                                                                                                                       aGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGc              >  4:91491‑M2/1‑148 (MQ=255)                                                                                                                                         cccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGctg            >  4:76288‑M2/1‑146 (MQ=255)                                                                                                                                          ccGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGctgc           <  8:78322‑M2/149‑5 (MQ=255)                                                                                                                                            gggCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgc         <  6:258070‑M1/149‑70 (MQ=255)                                                                                                                                            gggCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgc         <  7:87097‑M1/149‑70 (MQ=255)                                                                                                                                                cccAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTAcc     <  1:141710/149‑1 (MQ=255)                                                                                                                                                  cAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCg    >  1:295700/1‑148 (MQ=255)                                                                                                                                                   aGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGcacatcggtcgcagacgcgtaattggtggcgatgccaccgttgccatccagcgctgcctcgtactgcgcgtaatag  <  6:132525‑M1/149‑77 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AAAGTCATTAGCGAGCCATTGCTGGCCAACCTTGCCGCCAACGCCGCGCGCCTGAGCGACAGTGAATACAGCATGCGCCTGTTTCTGTCGAAGGCTTCGCCTACCGCTTACGAACAGAACGCAAGGCTATTGGCAGCCCGGGCCCCAGGGCTGCAGGTATTGCCGCTTGAAGAGCAAGGCTTCTTCCGGGCCTTCTGCCAGAGCCGCTATTACGCGCAGTACGAGGCAGCGCTGGATGGCAACGGTGGCATCGCCACCAATTACGCGTCTGCGACCGATGTGCTGCGCTACCGGA  >  NC_002947/4921446‑4921740

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