Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,982,838	T→A	7.2%	H119Q (CAT→CAA)	PP_3510 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,982,838	0	T	A	7.2%	243.7 / 7.4	83	H119Q (CAT→CAA)	PP_3510	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (2/4);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGT  >  NC_002947/3982690‑3982983                                                                                                                                                      |                                                                                                                                                   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                                                                                                                                       <  4:190939/148‑5 (MQ=255)               cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                       <  4:197249/148‑5 (MQ=255)               cTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                       <  3:181937/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:173617/1‑149 (MQ=255)                    aCTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGa                                                                                                                                 <  3:318687/149‑1 (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‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTg                                                                                                                                            <  8:152970/123‑1 (MQ=255)                                   aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTg                                                                                                                                            >  7:152970/1‑123 (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                                                                                                           >  7:362657/1‑143 (MQ=255)                                               cGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTAt                                                                                                           <  8:362657/143‑1 (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:41882/1‑149 (MQ=255)                                                      cAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTg                                                                                              >  6:35980/1‑149 (MQ=255)                                                      cAGATGGACCATACCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGTAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTg                                                                                              >  1:326411/1‑149 (MQ=255)                                                       aGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCa                                                                                                                 <  2:34199/129‑1 (MQ=255)                                                       aGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCa                                                                                                                 >  1:34199/1‑129 (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              <  3:210046/141‑1 (MQ=255)                                                                                                                                               aGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa              >  4:210046/1‑141 (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  <  3:142000/149‑1 (MQ=255)                                                                                                                                                   gCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt  <  5:402723/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTT‑GGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGC‑TTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGT  >  NC_002947/3982690‑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: ‑