Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,982,835	G→T	6.8%	G118G (GGG→GGT)	PP_3510 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,982,835	0	G	T	6.8%	259.8 / 6.4	83	G118G (GGG→GGT)	PP_3510	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/35);  new base T (2/4);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

GCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGT  >  NC_002947/3982689‑3982983                                                                                                                                                    |                                                                                                                                                      gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCa                                                                                                                                                     >  6:326457/1‑149 (MQ=255) gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCa                                                                                                                                                     <  1:238578/149‑1 (MQ=255)  cAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAt                                                                                                                                                    >  2:39299/1‑149 (MQ=255)   aaaGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATATAGCTGGACCATCCCGCGT‑GGCATGTACGCTAAGACAGCGAAGAGGCCTTTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTTCAGGGGCAtt                                                                                                                                                   >  1:104848/1‑149 (MQ=255)   aaaGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAACCGCCTGCAGGGGCAtt                                                                                                                                                   >  6:304624/1‑149 (MQ=255)      ggCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGAc                                                                                                                                                <  3:6293/149‑1 (MQ=255)           gTATGCTGGCACTGACCTGCAAGGCATTGGAACACGCCGTAATACAGCTGGACAAACCCGCAT‑GGCATGTACGCATTGACAGCGAAGAGGCCTGTGCACTGCTTCATGGCGTGGACTCGGGAAACCGCCTGCAGCGGCATTGACCCTGc                                                                                                                                           >  4:271871/1‑149 (MQ=255)              tGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCgg                                                                                                                                         <  5:326457/148‑1 (MQ=255)              tGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGt                                                                                                                                        <  6:284912/149‑1 (MQ=255)                cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                       <  3:181937/148‑5 (MQ=255)                cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                       <  4:197249/148‑5 (MQ=255)                cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                       <  4:190939/148‑5 (MQ=255)                cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                      <  7:334505/149‑1 (MQ=255)                cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                      <  3:135556/149‑1 (MQ=255)                cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                      >  2:110417/1‑149 (MQ=255)                 tGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACCGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGCAAACGCCTGCAGGGGCATTGACCCTGCGGgtt                                                                                                                                      >  7:269936/1‑148 (MQ=255)                   gCACTGACCTGCCAGTCATTGGACCAGGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGATGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTg                                                                                                                                     <  3:271871/147‑1 (MQ=255)                   gCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                    <  8:208785/148‑1 (MQ=255)                    cACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGAc                                                                                                                                   >  1:105642/1‑149 (MQ=255)                    cACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGAc                                                                                                                                   >  5:65594/1‑148 (MQ=255)                    cACTGACATGCCAGGCATTGGAGATCGAGGGAATTCATCTGGACCATCCCTCTT‑GGCATGTACGCTATGACAGCGTAGAGGCGTGTGCTCGGGTGCATGGCGGGGACTGGGGATAACGGCTGCAGGGGTATTGACCCTGCGGttttag                                                                                                                                   <  4:305761/148‑4 (MQ=255)                     aCTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGa                                                                                                                                 <  3:318687/149‑1 (MQ=255)                     aCTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGa                                                                                                                                 >  3:173617/1‑149 (MQ=255)                       tGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACg                                                                                                                                  <  6:278213/146‑1 (MQ=255)                       tGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACg                                                                                                                                  >  5:278213/1‑146 (MQ=255)                         aCCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATa                                                                                                                               <  1:39299/147‑1 (MQ=255)                           cTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGcag                                                                                                                                         <  6:102387/135‑3 (MQ=255)                           cTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGcag                                                                                                                                         >  5:102387/1‑133 (MQ=255)                            tGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATc                                                                                                                          >  1:249394/1‑149 (MQ=255)                                    aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAgggg                                                                                                                                                       <  4:384904/112‑1 (MQ=255)                                    aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAgggg                                                                                                                                                       >  3:384904/1‑112 (MQ=255)                                    aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTg                                                                                                                                            >  7:152970/1‑123 (MQ=255)                                    aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTg                                                                                                                                            <  8:152970/123‑1 (MQ=255)                                     ttGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCAc                                                                                                                >  4:83381/1‑149 (MQ=255)                                       ggACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCAcc                                                                                                               <  6:65594/148‑1 (MQ=255)                                                cGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCaaatcaaatc                                                                                                     >  4:297867/1‑149 (MQ=255)                                                cGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTAt                                                                                                           <  8:362657/143‑1 (MQ=255)                                                cGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTAt                                                                                                           >  7:362657/1‑143 (MQ=255)                                                cGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTc                                                                                                     <  2:105642/149‑1 (MQ=255)                                                  gAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGt                                                                                                   >  4:142000/1‑149 (MQ=255)                                                      aCAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACt                                                                                               >  2:394855/1‑149 (MQ=255)                                                       cAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTg                                                                                              >  6:35980/1‑149 (MQ=255)                                                       cAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTg                                                                                              >  6:41882/1‑149 (MQ=255)                                                       cAGATGGACCATACCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGTAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTg                                                                                              >  1:326411/1‑149 (MQ=255)                                                        aGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCa                                                                                                                 >  1:34199/1‑129 (MQ=255)                                                        aGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCa                                                                                                                 <  2:34199/129‑1 (MQ=255)                                                           tGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCg                                                                                          >  3:401749/1‑149 (MQ=255)                                                             gACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGt                                                                                         <  4:401749/148‑1 (MQ=255)                                                               ccATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                      <  7:179929/149‑1 (MQ=255)                                                                cATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGTCGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTccc                                                                                     >  7:253637/1‑149 (MQ=255)                                                                   cccGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGAt                                                                                  >  6:123668/1‑149 (MQ=255)                                                                      gCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCg                                                                               <  2:61189/149‑1 (MQ=255)                                                                               tGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGTGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                       <  5:41882/149‑1 (MQ=255)                                                                               tGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                       <  5:35980/149‑1 (MQ=255)                                                                                      tATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAtt                                                                  >  1:71577/1‑147 (MQ=255)                                                                                        ttaCAGCGAAGAGTCCTGTGCACGGGTGCATGGCGGGTACTTGGGAAAACGCCTGCAGGGGCATTGTCCCTGCGGTTTGACGATAGTATCGCATTCGCACCCTATCAATTCGTGACTGGGCGTTCCCGATCCGCCAAGGATGTATTAAg                                                               <  1:110417/147‑1 (MQ=255)                                                                                          aCAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGt                                                             <  2:71577/148‑1 (MQ=255)                                                                                                    ggCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgc                                                   <  1:394855/148‑1 (MQ=255)                                                                                                     gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcgc                                                 <  2:52387/149‑1 (MQ=255)                                                                                                     gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcg                                                  >  2:381811/1‑148 (MQ=255)                                                                                                        tgtgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCata                                              >  8:63754/1‑149 (MQ=255)                                                                                                            cACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                           <  1:381811/148‑1 (MQ=255)                                                                                                            cACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                           >  3:282356/1‑148 (MQ=255)                                                                                                            cACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAc                                          <  3:297867/149‑1 (MQ=255)                                                                                                               gggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcag                                        >  2:405525/1‑148 (MQ=255)                                                                                                                ggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagca                                      <  3:83381/149‑1 (MQ=255)                                                                                                                   gCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTc                                                                                                     <  3:272167/83‑1 (MQ=255)                                                                                                                   gCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTc                                                                                                     >  4:272167/1‑83 (MQ=255)                                                                                                                      tGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAgg                                                              <  2:277878/119‑1 (MQ=255)                                                                                                                      tGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAgg                                                              >  1:277878/1‑119 (MQ=255)                                                                                                                               cTGGGTAAAACTCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGAGCATATAACAGCAGCCCGAGCGATCggg                       <  2:326411/149‑1 (MQ=255)                                                                                                                                 ggggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGc                      >  7:335381/1‑148 (MQ=255)                                                                                                                                  gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                    >  3:298879/1‑149 (MQ=255)                                                                                                                                  gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                    >  1:137781/1‑149 (MQ=255)                                                                                                                                  gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGAGGCATGGTGCGCATATAACAGCAGCCCGAg                               >  2:34303/1‑138 (MQ=255)                                                                                                                                  gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGAGGCATGGTGCGCATATAACAGCAGCCCGAg                               <  1:34303/138‑1 (MQ=255)                                                                                                                                  gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTACCGATCAGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                    >  3:250950/1‑149 (MQ=255)                                                                                                                                    gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACtt                  >  7:164724/1‑149 (MQ=255)                                                                                                                                      aaaCGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTACCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCatat                                             >  5:318781/1‑120 (MQ=255)                                                                                                                                          gCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCATCCCGAGCGATCGGGCACTTGCAAGa            <  1:405525/149‑1 (MQ=255)                                                                                                                                          gCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGa            >  1:317847/1‑149 (MQ=255)                                                                                                                                          aaccgCAGGGTCAATGCCCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGa            >  3:113298/5‑149 (MQ=255)                                                                                                                                            cTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGAc           >  6:402723/1‑148 (MQ=255)                                                                                                                                                aGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa              >  4:210046/1‑141 (MQ=255)                                                                                                                                                aGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa              <  3:210046/141‑1 (MQ=255)                                                                                                                                                 ggggCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAg     >  7:106274/1‑149 (MQ=255)                                                                                                                                                  gggCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg    <  2:104848/149‑1 (MQ=255)                                                                                                                                                  gggCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg    <  2:116379/149‑1 (MQ=255)                                                                                                                                                    gCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt  <  5:402723/149‑1 (MQ=255)                                                                                                                                                    gCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt  <  3:142000/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGT  >  NC_002947/3982689‑3982983

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