Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,666,641	T→A	6.9%	L835Q (CTG→CAG)	acnA‑II →	aconitate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,666,641	0	T	A	6.9%	197.2 / 3.2	67	L835Q (CTG→CAG)	acnA‑II	aconitate hydratase
Reads supporting (aligned to +/- strand):  ref base T (25/37);  new base A (2/3);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TTCAAGCCGGGCACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAA  >  NC_002947/2666493‑2666787                                                                                                                                                     |                                                                                                                                                   ttCAAGCCGGGCACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCt                                                                                                                                                    >  7:308075/1‑149 (MQ=255)   cAAGCCGGGCACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTg                                                                                                                                                   <  4:59614/148‑1 (MQ=255)   cAAGCCGGGCACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTg                                                                                                                                                   <  2:236443/148‑1 (MQ=255)         gggCACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCAGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGTCCAGACGGCAGGTCAccg                                                                                                                                            <  5:342806/149‑1 (MQ=255)            cACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGTCCAGACGGCAGGTCAccg                                                                                                                                            <  5:97586/146‑1 (MQ=255)            cACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGTCCAGACGGCAGGTCAccg                                                                                                                                            >  6:97586/1‑146 (MQ=255)             accgaccgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGa                                                                                                                                        <  6:117636/149‑1 (MQ=255)              ccgaccgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAg                                                                                                                                       <  2:325181/149‑1 (MQ=255)               cgaccgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGTCCAGACGGCAGGTCAccggtactt                                                                                                                                      <  7:135382/149‑7 (MQ=255)                gaccgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCCCGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAgg                                                                                                                                      >  8:347523/1‑148 (MQ=255)                gaccgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAgg                                                                                                                                      >  2:398790/1‑148 (MQ=255)                 accgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    <  3:350095/149‑1 (MQ=255)                 accgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    <  1:244198/149‑1 (MQ=255)                 accgCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    <  7:68338/149‑1 (MQ=255)                    gCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCTGTCTGGACACCGCCGAGGAAgtg                                                                                                                                 >  1:28376/1‑149 (MQ=255)                       aGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTc                                                                                                                               >  8:227917/1‑148 (MQ=255)                         aCCCTGGGGCTGGATGGCAGCGAAACCTATGACGTGCTGTGTGAGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCGAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGAt                                                                                                                            <  7:273856/149‑1 (MQ=255)                                 ccTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    <  6:97806/133‑1 (MQ=255)                                 ccTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    >  5:97806/1‑133 (MQ=255)                                 ccTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    <  6:103499/133‑1 (MQ=255)                                 ccTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGaa                                                                                                                                    >  5:103499/1‑133 (MQ=255)                                 ccTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAg                                                                                                                    <  2:251277/149‑1 (MQ=255)                                      aTGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAgcgggc                                                                                                               <  2:79285/149‑1 (MQ=255)                                           aGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGt                                                                                                          >  2:232765/1‑149 (MQ=255)                                           aGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGt                                                                                                          >  4:205308/1‑149 (MQ=255)                                            gCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGtg                                                                                                         >  1:356144/1‑149 (MQ=255)                                                            gCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGcc                                                                                          >  3:792/1‑148 (MQ=255)                                                              tGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGGGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAg                                                                                       >  4:182318/1‑149 (MQ=255)                                                                ggTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGa                                                                                     <  4:792/149‑1 (MQ=255)                                                                  tGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACt                                                                                   <  8:308075/149‑1 (MQ=255)                                                                    cgcgTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCTGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTc                                                                                 <  2:28376/149‑1 (MQ=255)                                                                       gaaCGCCACGGGCGACCCTCACCCTGCTGGTTACCCGTGCCAATGGCGAGTGCCTCGAAGTACAGCACACCTGCCGTCAGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACATCCTc                                                                              >  8:48103/3‑149 (MQ=255)                                                                           gCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAg                                                                          >  7:308425/1‑149 (MQ=255)                                                                           gCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAg                                                                          >  4:73041/1‑149 (MQ=255)                                                                            ccACGGGCGACCCTGCCCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGc                                                                         <  6:13524/149‑1 (MQ=255)                                                                             cACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCttt                                                                                             >  5:270695/1‑128 (MQ=255)                                                                             cACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCttt                                                                                             <  6:270695/128‑1 (MQ=255)                                                                               cGGTCGACACTGACCCTGGTGCTTACCCGTGCCAAGGGCGAGTGCCTGGAAGTACCGGTGACATGCCGGCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCTGGCGGGGTGCTGCAGCGCTTTTCCCAGGACTTCCTCGAAGCGAc                                                                      <  4:300351/149‑1 (MQ=255)                                                                                  gCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGc                                                                   >  1:354926/1‑149 (MQ=255)                                                                                        ctgaccctgGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAg                                                              >  8:185681/1‑148 (MQ=255)                                                                                         tgagcgtggtggTTACCCGTGCCAATGGCGAGTGCCTGGAAGTAGGGGTGACCTGGGGTCTGGACTGCGGGGTGGTAGTGTGGTTGTACGATGCGGGCGGGGTGCTGCTGCGCTTTGCCGAGGTCTTGCTCGAAGGGACCGCCTGAGCa                                                            <  7:48103/143‑1 (MQ=255)                                                                                               tggtggTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGaaag                                                      >  1:144673/1‑147 (MQ=255)                                                                                                 gtggtTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAAt                                                       <  2:144673/146‑1 (MQ=255)                                                                                                   ggtTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCgggg                                                  <  2:354926/149‑1 (MQ=255)                                                                                                    gtTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCgggg                                                  <  8:308425/148‑1 (MQ=255)                                                                                                     ttACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGc                                                <  7:227917/149‑1 (MQ=255)                                                                                                            tGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAg                                         <  1:398790/149‑1 (MQ=255)                                                                                                              ccAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCg                                       <  3:73041/149‑1 (MQ=255)                                                                                                               cAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCAt                                                           >  6:18635/1‑128 (MQ=255)                                                                                                               cAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCAt                                                           <  5:18635/128‑1 (MQ=255)                                                                                                                 aTGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGTAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAg                                         >  8:181574/1‑144 (MQ=255)                                                                                                                 aTGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGTAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAg                                         <  7:181574/144‑1 (MQ=255)                                                                                                                    gCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGTGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCa                                 >  4:98232/1‑149 (MQ=255)                                                                                                                      gAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAgg                               <  3:205308/149‑1 (MQ=255)                                                                                                                      gAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAgg                               <  2:356144/149‑1 (MQ=255)                                                                                                                         tGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGACGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAATGCGGCCCCAGGGTa                            <  7:347523/149‑1 (MQ=255)                                                                                                                            cTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTAttt                         >  6:347265/1‑149 (MQ=255)                                                                                                                               gAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCttt                                                                                             <  5:414360/78‑1 (MQ=255)                                                                                                                               gAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCttt                                                                                             >  6:414360/1‑78 (MQ=255)                                                                                                                               gAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCAc                      <  1:232765/149‑1 (MQ=255)                                                                                                                                 aGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCg                    <  1:311909/149‑1 (MQ=255)                                                                                                                                      cGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACTGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCACGGTATTTCACCGGACTa               >  1:80428/1‑149 (MQ=255)                                                                                                                                      cGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAATCGGCCCCAAGGTATTTAACCGGACTa               >  1:579/1‑149 (MQ=255)                                                                                                                                             cTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAg                                <  1:323406/125‑1 (MQ=255)                                                                                                                                             cTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGAt        <  3:182318/149‑1 (MQ=255)                                                                                                                                              tGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCCCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATg       >  7:333676/1‑149 (MQ=255)                                                                                                                                                  gTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCTTTGCGGccccc                                 <  2:178292/119‑2 (MQ=255)                                                                                                                                                   tCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCaa  <  5:347265/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TTCAAGCCGGGCACCGACCGCAAGACCCTGGGCCTGGATGGCAGCGAAACCTATGACGTGCTGGGTGCGCGTACGCCACGGGCGACCCTGACCCTGGTGGTTACCCGTGCCAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAA  >  NC_002947/2666493‑2666787

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