Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	464,320	G→C	7.8%	H259Q (CAC→CAG) ‡	davB ←	lysine 2‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	464,320	0	G	C	7.8%	182.5 / 4.7	62	H259Q (CAC→CAG) ‡	davB	lysine 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGGCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGCG  >  NC_002947/464172‑464465                                                                                                                                                        |                                                                                                                                                  tggcggcgcaccgcgtaccggggtcaagcgcattgcccgcgccgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAgg                                                                                                                                                   >  7:86570‑M2/47‑149 (MQ=255)   GCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTg                                                                                                                                                 <  2:269758‑M1/149‑1 (MQ=255)     GGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCGGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACcagggt                                                                                                                                                 <  3:280962‑M1/147‑3 (MQ=255)      GCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCtg                                                                                                                                              >  7:268710‑M1/1‑149 (MQ=255)         CGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCtggtg                                                                                                                                           >  4:230090‑M1/1‑149 (MQ=255)          gggCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCtggtgg                                                                                                                                          >  7:160120/1‑149 (MQ=255)             cAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcg                                                                                                                                       <  2:295496/149‑1 (MQ=255)                  cgcTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTGCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCa                                                                                                                                  <  5:231539/149‑1 (MQ=255)                        aCCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCTCTTGGt                                                                                                                            >  7:298716/1‑149 (MQ=255)                         ccgcGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtc                                                                                                                           <  5:212838/146‑1 (MQ=255)                          cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGtcgtcg                                                                                                                                    <  8:187125/139‑1 (MQ=255)                          cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGtcgtcg                                                                                                                                    >  7:187125/1‑139 (MQ=255)                          cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtca                                                                                                                          >  5:118436/1‑149 (MQ=255)                           ccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtcat                                                                                                                         <  3:52871/149‑1 (MQ=255)                               tACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATcacc                                                                                                                     <  6:241007/149‑1 (MQ=255)                                 ttgcccgcgccgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATcaccac                                                                                                                   <  8:86570‑M2/135‑1 (MQ=255)                                  gcgGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcg                                                                                                                  >  4:119187/1‑149 (MQ=255)                                        cgccgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGaa                                                                                                             <  3:82582‑M2/141‑1 (MQ=255)                                               atggcgCGTGCTCAGGGGTGGGCCCTCTGGCCAAATCACGGATCGGTTCGGTACGTCGGGCCATATTCCTTGTT‑GGACCAAAGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTc                                                                                                       <  1:100135/144‑1 (MQ=255)                                                    gCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                                <  6:118436/149‑1 (MQ=255)                                                      gTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                                >  6:238487/1‑147 (MQ=255)                                                         cTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCa                                                                                                     >  4:207322/1‑139 (MQ=255)                                                         cTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCa                                                                                                     <  3:207322/139‑1 (MQ=255)                                                           cAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTgg                                                                                         >  3:22759/1‑149 (MQ=255)                                                            aGGCTTGTGCCCGCTGGCCAATGCACGCAGCATTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCTACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGCGTTgg                                                                                         <  2:233197/148‑1 (MQ=255)                                                            aGGCTTGTGCACTGTGTCCAATGAACTCATCTTTCCGGTAAGTCGCGTCAGATGCCTTGTG‑GCAC‑‑TTGTTGCAGGCTCCCGAGGACCTGTTGCTGGTGGTCTTCGCAGTTGGTCTTTACCACTTGTAAGATTTCCAGCATCTAGTTgg                                                                                         <  8:160120/148‑1 (MQ=255)                                                                   tggccTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCGTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTcc                                                                                 <  4:280962/146‑1 (MQ=255)                                                                     cccTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGa                                                                               <  1:269758/149‑1 (MQ=255)                                                                          tGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTc                                                                            <  4:17058/147‑1 (MQ=255)                                                                             ccAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAc                                                                       <  5:235386/149‑1 (MQ=255)                                                                                     cGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTg                                                               >  5:260768/1‑149 (MQ=255)                                                                                      gCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGTGCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTg                                                               >  1:48880/1‑149 (MQ=255)                                                                                      gCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAccgcc                                                                   >  8:129999/1‑144 (MQ=255)                                                                                      gCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAccgcc                                                                   <  7:129999/144‑1 (MQ=255)                                                                                      gCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGc                                                              >  4:112858/1‑149 (MQ=255)                                                                                                ggTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGAc                                                    >  4:9645/1‑149 (MQ=255)                                                                                                 gTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGt                                                                             <  1:305700/123‑1 (MQ=255)                                                                                                 gTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGt                                                                             >  2:305700/1‑123 (MQ=255)                                                                                                         cgcCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                           <  3:112858/149‑1 (MQ=255)                                                                                                         cgcCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                           <  7:45357/149‑1 (MQ=255)                                                                                                                tGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTcg                                    <  3:119187/149‑1 (MQ=255)                                                                                                                   cTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTGGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCgg                                 >  1:225980/1‑149 (MQ=255)                                                                                                                          gCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGcc                                                             <  6:195942/115‑1 (MQ=255)                                                                                                                          gCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGcc                                                             >  5:195942/1‑115 (MQ=255)                                                                                                                          gCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTg                            <  6:260768/148‑1 (MQ=255)                                                                                                                            aC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGa                                                                               <  3:35946/95‑1 (MQ=255)                                                                                                                            aC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGa                                                                               >  4:35946/1‑95 (MQ=255)                                                                                                                            aC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTgg                           >  3:94217/1‑147 (MQ=255)                                                                                                                             c‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAg                        >  4:125691/1‑149 (MQ=255)                                                                                                                                   ttCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCt                      >  2:173467/1‑147 (MQ=255)                                                                                                                                   ttCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGAAGTTGGTCATCACCACGCGGAAGATTTCCATCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCa                    <  8:298716/149‑1 (MQ=255)                                                                                                                                     ccACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGc                  >  2:32470/1‑149 (MQ=255)                                                                                                                                      cACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGaagctaagct                 >  7:109879/1‑149 (MQ=255)                                                                                                                                      cACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCt                 >  2:14835/1‑149 (MQ=255)                                                                                                                                      cACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCt                 <  2:48880/149‑1 (MQ=255)                                                                                                                                      cACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCt                 >  3:135699/1‑149 (MQ=255)                                                                                                                                         gCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg              <  1:173467/149‑1 (MQ=255)                                                                                                                                         gCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg              >  8:216038/1‑149 (MQ=255)                                                                                                                                           cgacGACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg              <  5:180145/144‑1 (MQ=255)                                                                                                                                           cgacGACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg              >  6:180145/4‑147 (MQ=255)                                                                                                                                             ccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCaa           >  8:81842/1‑148 (MQ=255)                                                                                                                                               gACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGt                                                                <  6:203049/93‑1 (MQ=255)                                                                                                                                               gACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGt                                                                >  5:203049/1‑93 (MQ=255)                                                                                                                                                     caGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgcgc   >  7:75487/1‑148 (MQ=255)                                                                                                                                                        gTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGcgcg  <  1:280995/146‑1 (MQ=255)                                                                                                                                                        gTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGcgcg  <  1:285486/146‑1 (MQ=255)                                                                                                                                                        gTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGcgcg  >  2:280995/1‑146 (MQ=255)                                                                                                                                                        gTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGcgcg  >  2:285486/1‑146 (MQ=255)                                                                                                                                                        |                                                                                                                                                  CGGCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTG‑GCAC‑‑TTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGCG  >  NC_002947/464172‑464465

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