Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	464,320	G→C	6.9%	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	6.9%	275.9 / 8.1	87	H259Q (CAC→CAG) ‡	davB	lysine 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (38/43);  new base C (2/4);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GGCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGCGAGG  >  NC_002947/464173‑464468                                                                                                                                                    |                                                                                                                                                     ggcggcGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGt                                                                                                                                                     >  3:280064/1‑149 (MQ=255)    ggcgCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCt                                                                                                                                                  <  8:145587/149‑1 (MQ=255)     gcgcGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtg                                                                                                                                                 <  6:301744/149‑1 (MQ=255)       gcgGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtggt                                                                                                                                               >  3:42416/1‑149 (MQ=255)        cgGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtggtg                                                                                                                                              <  6:277866/149‑1 (MQ=255)           gCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtc                                                                                                                                           <  1:204412/149‑1 (MQ=255)             aaTGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcgt                                                                                                                                         >  4:196376/1‑149 (MQ=255)             aaTGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcgt                                                                                                                                         <  2:223382/149‑1 (MQ=255)             aaTGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                         <  4:183888/149‑2 (MQ=255)             aaTGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                         <  1:301516/149‑2 (MQ=255)              aTGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcggg                                                                                                                                        <  5:302738/149‑3 (MQ=255)               tGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGtcgtcg                                                                                                                                       >  6:27100/1‑149 (MQ=255)                  gcTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAg                                                                                                                                    >  3:289341/1‑149 (MQ=255)                   cTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGt                                                                                                                                   <  1:169422/149‑1 (MQ=255)                      gACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGtt                                                                                                                                  <  7:76991/147‑1 (MQ=255)                      gACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGtt                                                                                                                                  <  4:280064/147‑1 (MQ=255)                      gACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTg                                                                                                                                 <  2:131867/148‑1 (MQ=255)                       aCCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCTCGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTgg                                                                                                                                <  4:42416/148‑1 (MQ=255)                       aCCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGt                                                                                                                               <  7:294921/149‑1 (MQ=255)                       aCCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTAGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTgg                                                                                                                                <  1:242481/148‑1 (MQ=255)                        ccccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCTTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtc                                                                                                                              <  4:204551/149‑1 (MQ=255)                        ccccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACACAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGt                                                                                                                               >  6:322201/1‑148 (MQ=255)                         cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtca                                                                                                                             >  5:246248/1‑149 (MQ=255)                           cGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtcat                                                                                                                            >  6:170164/1‑148 (MQ=255)                             gTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATcac                                                                                                                         <  5:170164/149‑1 (MQ=255)                                 gcgGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcg                                                                                                                     >  5:99951/1‑149 (MQ=255)                                 gcgGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcg                                                                                                                     <  2:133170/149‑1 (MQ=255)                                   ggTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcgcg                                                                                                                   >  8:311309/1‑149 (MQ=255)                                       cgccgcCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGaa                                                                                                                <  6:194693/148‑1 (MQ=255)                                       cgccgcCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGaa                                                                                                                <  4:289341/148‑1 (MQ=255)                                          cgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGAtt                                                                                                            >  4:81027/1‑149 (MQ=255)                                            ccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTc                                                                                                          <  3:107208/149‑1 (MQ=255)                                            ccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTc                                                                                                          <  4:182957/149‑1 (MQ=255)                                            ccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTCGTGTTGTACGACGCAGTTGTAAATCACCACCCGGAACATTTc                                                                                                          >  2:62779/1‑149 (MQ=255)                                              aTGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCa                                                                                                        <  1:177839/149‑1 (MQ=255)                                                  aGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGCAAGATTTcatgcaa                                                                                                    >  4:77678/1‑143 (MQ=255)                                                   gCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATc                                                                                                   <  2:113120/149‑1 (MQ=255)                                                     gTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGa                                                                                                 <  2:23861/149‑1 (MQ=255)                                                              cTGGTGCCCTCTGGCCAATGCACGCAGCTTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCAAGCGCGCTACCTCCAGGTGCTGGTGGTCGTCGCATTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGaa                                                                                        <  6:199369/149‑1 (MQ=255)                                                               tGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAg                                                                                       >  4:94491/1‑149 (MQ=255)                                                               tGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCAAGCATCGAGTTGGGGAAg                                                                                       >  7:74005/1‑149 (MQ=255)                                                                ggTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACAACGCGGAAGATTTCCAGCATCGAGTTGGGGaatt                                                                                      >  7:279835/1‑147 (MQ=255)                                                                    cccTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGa                                                                                  <  3:81027/149‑1 (MQ=255)                                                                    cccTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGa                                                                                  <  6:246248/149‑1 (MQ=255)                                                                     actctGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAg                                                                                 >  8:195128/2‑149 (MQ=255)                                                                         ttggcAATGCAAGCAGCGTTCCGGTAAGTGGCGCCATTTTCCTTGTGGGACTTTTTCCAAGGCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTcc                                                                              <  1:62779/144‑1 (MQ=255)                                                                           gCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGCAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCaa                                                                           >  3:40651/1‑149 (MQ=255)                                                                              aaTGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAccg                                                                        >  6:128009/1‑149 (MQ=255)                                                                                   aCGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCACGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGt                                                                   >  4:75532/1‑149 (MQ=255)                                                                                     gCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGATGAATTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAAACGCCGGTGc                                                                 >  8:13552/1‑149 (MQ=255)                                                                                          gTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGcc                                                                >  1:340314/1‑145 (MQ=255)                                                                                          gTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGcc                                                                <  2:340314/145‑1 (MQ=255)                                                                                               ggTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGAc                                                       <  7:311309/149‑1 (MQ=255)                                                                                               ggTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGAc                                                       <  3:75532/149‑1 (MQ=255)                                                                                                    gTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGc                                                  >  5:153201/1‑149 (MQ=255)                                                                                                        cgcCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                              >  1:164910/1‑149 (MQ=255)                                                                                                             gATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                              >  6:189846/1‑144 (MQ=255)                                                                                                             gATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                              <  5:189846/144‑1 (MQ=255)                                                                                                                 ccTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTcgcg                                     >  1:249444/1‑149 (MQ=255)                                                                                                                      tgGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCt                                <  2:249444/149‑1 (MQ=255)                                                                                                                      tgGCACTTGTTCCACGCCCCCGACAACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCt                                <  8:74005/149‑1 (MQ=255)                                                                                                                          aCTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGaa                            <  6:153201/149‑1 (MQ=255)                                                                                                                          aCTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGaa                            <  5:128009/149‑1 (MQ=255)                                                                                                                          aCTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGaa                            <  3:196376/149‑1 (MQ=255)                                                                                                                           cTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAg                           <  6:99951/149‑1 (MQ=255)                                                                                                                           cTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAg                           >  2:231351/1‑149 (MQ=255)                                                                                                                                 ccACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGc                     >  6:116895/1‑149 (MQ=255)                                                                                                                                  cACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCATCATCGAGTTGGCGAAGTCCGAGTCCCAACCGCCGGTGCAGAAACCGCCCTGGCCGAACACGTCGAGGTGCTTGAAGCTCAGCt                    >  4:221272/1‑149 (MQ=255)                                                                                                                                   ccGCCCCCGACCACCAAGTGCTGGTGGTCGTCGCAGTTGGTCATCACCCCGCGGAAGATTTCCAGAATCGAGGTGGGGCAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGCCCGTGCCGACCCCTGAGCGGTGCTGCAATCTCAGCtt                   >  6:267031/2‑149 (MQ=255)                                                                                                                                     gCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg                 <  2:164910/149‑1 (MQ=255)                                                                                                                                     gCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg                 <  4:196920/149‑1 (MQ=255)                                                                                                                                     gCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg                 <  5:27100/149‑1 (MQ=255)                                                                                                                                       ccccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTcc                                                                                    <  2:333814/80‑1 (MQ=255)                                                                                                                                       ccccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTcc                                                                                    >  1:333814/1‑80 (MQ=255)                                                                                                                                       ccccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGc                     >  6:105990/1‑143 (MQ=255)                                                                                                                                       ccccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGc                     <  5:105990/143‑1 (MQ=255)                                                                                                                                       ccccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGc                >  7:204993/1‑148 (MQ=255)                                                                                                                                       ccccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCAGCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGAAAGCTCAGCTTGGc                >  5:293986/1‑148 (MQ=255)                                                                                                                                         ccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGaa                                                                                                                >  7:122878/1‑50 (MQ=255)                                                                                                                                         ccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGaa                                                                                                                <  8:122878/50‑1 (MQ=255)                                                                                                                                         ccGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTAATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGAGCCCAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCaa              >  7:160198/1‑148 (MQ=255)                                                                                                                                         acGACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCaa              >  7:286485/2‑148 (MQ=255)                                                                                                                                           gACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCaa              >  3:231054/1‑146 (MQ=255)                                                                                                                                              caccaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAaga          <  1:73012/147‑1 (MQ=255)                                                                                                                                              caccaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAaga          >  2:73012/1‑147 (MQ=255)                                                                                                                                               accaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTgg                                                                                            <  3:267986/64‑1 (MQ=255)                                                                                                                                               accaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTgg                                                                                            >  4:267986/1‑64 (MQ=255)                                                                                                                                                    gTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGCGAgg  >  8:111917/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGCGAGG  >  NC_002947/464173‑464468

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