Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,193,304	T→C	6.2%	N116D (AAC→GAC)	upxB ←	phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,193,304	0	T	C	6.2%	182.8 / 2.3	64	N116D (AAC→GAC)	upxB	phosphatase
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

TCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCACCGGTGG  >  NC_002947/1193162‑1193451                                                                                                                                                |                                                                                                                                                    tCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGtt                                                                                                                                                     >  5:348930/1‑143 (MQ=255) tCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGtt                                                                                                                                                     <  6:348930/143‑1 (MQ=255) tCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTgg                                                                                                                                               >  6:83215/1‑149 (MQ=255) tCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTgg                                                                                                                                               >  6:86021/1‑149 (MQ=255) tCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTgg                                                                                                                                               >  3:171204/1‑149 (MQ=255) tCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGGTCAGTgg                                                                                                                                               >  6:220665/1‑149 (MQ=255)   gAAATTGATCGCCAGCAGTCCTTCTTCCGATCTGCCGTTTATGGGGTAGAAGTGCATGCCGTCGTGGTTCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCa                                                                                                                                                   <  8:216291/143‑1 (MQ=255)   gAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGc                                                                                                                                              <  8:2518/148‑1 (MQ=255)   gAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGACCATGGAGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGttga                                                                                                                                                   >  7:216291/1‑141 (MQ=255)              ccAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGgc                                                                                                                                  <  5:220665/149‑1 (MQ=255)               cagcagGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGgcg                                                                                                                                 >  6:253320/1‑149 (MQ=255)                   agGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAg                                                                                                                             >  4:271508/1‑149 (MQ=255)                     gCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGc                                                                                                                                              <  1:137868/130‑1 (MQ=255)                             ccGAGGTGCCGTTGATGTGGAATAAGTGCATGGGGTGGTGGTTCATGCCCTTTGCGTTGGCCTGGTCGCTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCTGTGCGCCAGCACGTAg                                                                                                                      <  1:111459/146‑1 (MQ=255)                               gAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCt                                                                                                                 >  7:10979/1‑149 (MQ=255)                                aGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTg                                                                                                                <  8:108495/149‑1 (MQ=255)                                       gTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCcgg                                                                                                         >  1:148892/1‑149 (MQ=255)                                          gATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCcggcgg                                                                                                      >  6:100781/1‑149 (MQ=255)                                          gATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCcggcg                                                                                                       <  5:83215/148‑1 (MQ=255)                                          gATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCcggcg                                                                                                       <  5:253320/148‑1 (MQ=255)                                           aTGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCggcggc                                                                                                     >  3:170127/1‑149 (MQ=255)                                                  agaagTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACca                                                                                                <  3:327928/147‑1 (MQ=255)                                                  agaagTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACca                                                                                                >  4:327928/1‑147 (MQ=255)                                                  agaagTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCAGCGGCGACca                                                                                                <  8:22494/147‑1 (MQ=255)                                                  agaagTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCAGCGGCGACca                                                                                                >  7:22494/1‑147 (MQ=255)                                                   gaagTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACcacac                                                                                             >  4:173319/1‑149 (MQ=255)                                                     agggCATGCCGTCGTGGTGCATGCCCATGGCGGTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAg                                                                                           <  4:169867/146‑1 (MQ=255)                                                           tGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCg                                                                                     <  8:10979/149‑1 (MQ=255)                                                               gtcgtGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTgc                                                                                 <  5:86021/149‑1 (MQ=255)                                                                 cgtGGGGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTgcgc                                                                               <  5:96775/149‑1 (MQ=255)                                                                    ggtgCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGcgcgcg                                                                            >  7:52963/1‑147 (MQ=255)                                                                     gtgCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGGGc                                                                           >  3:131681/1‑149 (MQ=255)                                                                         aTGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTgg                                                                       <  3:271508/149‑1 (MQ=255)                                                                            cccATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGa                                                                     <  4:170127/148‑1 (MQ=255)                                                                             ccATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATg                                                                   >  5:132005/1‑149 (MQ=255)                                                                             ccATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATg                                                                   >  2:310303/1‑149 (MQ=255)                                                                                     ttGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGaa                                                           <  1:310303/149‑1 (MQ=255)                                                                                       ggCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGc                                                         >  2:34220/1‑149 (MQ=255)                                                                                         ccTGGTCGGTCGAGGTGTTGCTGCCGTTCGGACTTCCACGGGTTGCCGTTACTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGAAAGcc                                                        >  1:186818/1‑149 (MQ=255)                                                                                           tGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTCCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAg                                                     >  2:313758/1‑149 (MQ=255)                                                                                           tGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTCCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAg                                                     >  4:335880/1‑149 (MQ=255)                                                                                                  tcgAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGt                                                    <  3:232827/143‑1 (MQ=255)                                                                                                  tcgAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGt                                                    >  4:232827/1‑143 (MQ=255)                                                                                                  tcgAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGt                                              >  2:306825/1‑149 (MQ=255)                                                                                                    gAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGttt                                            <  7:325850/149‑1 (MQ=255)                                                                                                    gAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGttt                                            <  7:162790/149‑1 (MQ=255)                                                                                                    gAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGtt                                             <  6:132005/148‑1 (MQ=255)                                                                                                       gtgtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCgg                                         <  8:81015/149‑1 (MQ=255)                                                                                                       gtgtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCgg                                         >  7:153601/1‑149 (MQ=255)                                                                                                       gtgtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTCCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCgg                                         <  1:313758/149‑1 (MQ=255)                                                                                                       gtgtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTCCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCgg                                         <  3:335880/149‑1 (MQ=255)                                                                                                        tgtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGc                                        <  1:34220/149‑1 (MQ=255)                                                                                                         gtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTg                                                                        >  1:187510/1‑116 (MQ=255)                                                                                                         gtTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTg                                                                        <  2:187510/116‑1 (MQ=255)                                                                                                            gctgcCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGAc                                                                                                  >  5:315363/1‑87 (MQ=255)                                                                                                            gctgcCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGAc                                                                                                  <  6:315363/87‑1 (MQ=255)                                                                                                               gcCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGAc                                                                <  5:132253/118‑1 (MQ=255)                                                                                                               gcCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGAc                                                                >  6:132253/1‑118 (MQ=255)                                                                                                                ccG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGGGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCgg                                <  4:131681/149‑1 (MQ=255)                                                                                                                ccG‑TCGGACTTCCACGGGGTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTcc                                  <  3:243427/147‑1 (MQ=255)                                                                                                                    tCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCt                                                                                                                 >  5:326301/1‑65 (MQ=255)                                                                                                                    tCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCt                                                                                                                 <  6:326301/65‑1 (MQ=255)                                                                                                                              aCGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGc                   <  3:173319/149‑1 (MQ=255)                                                                                                                                      ccGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGt           <  5:295073/149‑1 (MQ=255)                                                                                                                                      ccGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGt           <  1:306825/149‑1 (MQ=255)                                                                                                                                        gTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCa         >  4:10687/1‑149 (MQ=255)                                                                                                                                        gTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCa         >  7:250045/1‑149 (MQ=255)                                                                                                                                               ttCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCACCGGTgg  <  8:153601/149‑1 (MQ=255)                                                                                                                                               ttCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCACCGGTg   <  5:100781/148‑1 (MQ=255)                                                                                                                                               ttCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGAAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCACCGGTg   <  2:186818/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    TCGAAATTGATCGCCAGCAGGCCGTCTTCCGAGCTGCCGTTGATGGGGAAGAAGTGCATGCCGTCGTGGTGCATGCCCATGGCGTTGGCCTGGTCGGTCGAGGTGTTGCTGCCG‑TCGGACTTCCACGGGTTGCCGTTGCTGTTCAGTGGCGTACCCCAGGGCGCCAGCACGTAGGCGCTGTAGCCGGCGGCGACCACACAGGCGTCGGTGCGCGAGCCTGGGATGGACTGGAAGCCCAGTGCCAGTTTCGGCACTTCCGGCTCGGTGACTGGCGGTTCGGTCACCGGTGG  >  NC_002947/1193162‑1193451

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑