Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,489,941	T→C	5.7%	Q498R (CAG→CGG)	clpV ←	protein ClpV1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,489,941	0	T	C	5.7%	206.7 / 2.2	70	Q498R (CAG→CGG)	clpV	protein ClpV1
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAG  >  NC_002947/3489794‑3490087                                                                                                                                                    |                                                                                                                                                   ccgcTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTgg                                                                                                                                                  <  2:177083/150‑1 (MQ=255) ccgcTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTgg                                                                                                                                                  <  2:423936/150‑1 (MQ=255)  cgcTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGACTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTggg                                                                                                                                                 <  2:2163424/150‑1 (MQ=255)      gTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTgg                                                                                                                                                  <  2:2005160/145‑1 (MQ=255)      gTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTgg                                                                                                                                                  >  1:2005160/1‑145 (MQ=255)       tCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACCGGGCATTg                                                                                                                                            >  1:1341068/1‑150 (MQ=255)       tCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACCGGGCATTg                                                                                                                                            >  1:1341073/1‑150 (MQ=255)        cGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGa                                                                                                                                           <  1:975845/150‑1 (MQ=255)         gACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAg                                                                                                                                          <  1:679804/150‑1 (MQ=255)         gACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGATCCAGGTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAg                                                                                                                                          >  1:555036/1‑150 (MQ=255)           ccAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgct                                                                                                                                        <  2:771055/150‑1 (MQ=255)            cAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctg                                                                                                                                       >  1:2032164/1‑150 (MQ=255)             aGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgc                                                                                                                                         >  2:1145203/1‑147 (MQ=255)             aGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgc                                                                                                                                         <  1:1145203/147‑1 (MQ=255)              gcgcGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctgct                                                                                                                                     <  1:1800351/150‑1 (MQ=255)              gcgcGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCACCGTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAgctgct                                                                                                                                     <  1:624020/150‑1 (MQ=255)                gcgTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGt                                                                                                                                   <  1:857910/150‑1 (MQ=255)                  gTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTc                                                                                                                                 <  1:80927/150‑1 (MQ=255)                     gACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCctg                                                                                                                              <  1:195672/150‑1 (MQ=255)                        cAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTAAGCTGCTGTTCctgctg                                                                                                                           <  2:2342404/150‑1 (MQ=255)                               tGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCg                                                                                                                    >  1:1923444/1‑150 (MQ=255)                                 gTCGCCTTGCACGCGCGCCAGCTCGGCAGGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGc                                                                                                                  <  2:2032164/150‑1 (MQ=255)                                     ccTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGc                                                                                                              <  2:2432793/150‑1 (MQ=255)                                     ccTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGc                                                                                                              >  2:430729/1‑150 (MQ=255)                                       ttGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAg                                                                                                            <  2:1740107/150‑1 (MQ=255)                                           ccgcgcgCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCa                                                                                                        <  2:393033/149‑1 (MQ=255)                                           aCGCGCGCCAGCGCGGCCCGCAGCGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCa                                                                                                        <  2:179791/150‑1 (MQ=255)                                             gcgcgcCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCCGCTCGGCCTGCATGTCAGc                                                                                                      <  2:1474309/150‑1 (MQ=255)                                                 gcCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGc                                                                                                  >  1:29030/1‑150 (MQ=255)                                                   cAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAg                                                                                                <  1:1746482/150‑1 (MQ=255)                                                        cGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGAtt                                                                                           >  1:834379/1‑150 (MQ=255)                                                            aCGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAg                                                                                       >  1:2138452/1‑150 (MQ=255)                                                             ctcAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGATTCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAgg                                                                                      <  1:500265/148‑1 (MQ=255)                                                              gCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTACCGGCTCCTTCTCGGTGTGGAGGGCGGCGTGCAGCGGTTACAggg                                                                                     >  2:2230088/1‑149 (MQ=255)                                                                aGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTTTTCCCGCTGCTGCTCGGCCTGCAGGGCAGCCTCCAGCGTTTCCAggcaa                                                                                   >  2:500265/1‑148 (MQ=255)                                                                           ctgctgGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACCGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTc                                                                        <  2:1341068/150‑1 (MQ=255)                                                                           ctgctgGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACCGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTc                                                                        <  2:1341073/150‑1 (MQ=255)                                                                               tgGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGc                                                                    <  2:940968/150‑1 (MQ=255)                                                                                  gCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGAAGACACCGATTCCAGGCGGCGGGCATGGTCGTGcccc                                                                 >  2:189799/1‑150 (MQ=255)                                                                                      cGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGc                                                             <  2:2320425/150‑1 (MQ=255)                                                                                      cGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGc                                                             <  1:1886746/150‑1 (MQ=255)                                                                                       ggCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGcc                                                            <  2:773985/150‑1 (MQ=255)                                                                                        gCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCg                                                           <  2:653060/150‑1 (MQ=255)                                                                                          aTCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCATGCGGCGGGCATGGTCGTGCCCCTTGCCGGc                                                         >  2:1308515/1‑150 (MQ=255)                                                                                          aTCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGc                                                         >  2:899362/1‑150 (MQ=255)                                                                                             gTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTc                                                      >  1:1712035/1‑150 (MQ=255)                                                                                              ttGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCg                                                     >  1:825734/1‑150 (MQ=255)                                                                                               tGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGt                                                    <  2:795304/150‑1 (MQ=255)                                                                                                   ggcgTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGcc                                                >  2:2024040/1‑150 (MQ=255)                                                                                                   ggcgTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGcc                                                >  2:2238791/1‑150 (MQ=255)                                                                                                     cgTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGa                                              <  2:449173/150‑1 (MQ=255)                                                                                                          aGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACg                                         >  1:2430624/1‑150 (MQ=255)                                                                                                           gcgcCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACgg                                        >  1:2246281/1‑149 (MQ=255)                                                                                                                  gAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTc                                 >  2:56358/1‑150 (MQ=255)                                                                                                                     cTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGc                              <  2:1923444/150‑1 (MQ=255)                                                                                                                      tGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAg                                                                                       >  1:1334369/1‑92 (MQ=255)                                                                                                                      tGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAg                                                                                       <  2:1334369/92‑1 (MQ=255)                                                                                                                          cGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCa                                                                                                 <  1:2110019/78‑1 (MQ=255)                                                                                                                          cGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCa                                                                                                 >  2:2110019/1‑78 (MQ=255)                                                                                                                           gACNAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGAAACCATCGCCTAGCACGTCGATTTCCGCCTGCAg                        >  2:1496175/1‑150 (MQ=255)                                                                                                                                gTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGt                   >  1:948073/1‑150 (MQ=255)                                                                                                                                gTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGt                   >  2:109747/1‑150 (MQ=255)                                                                                                                                 ttCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGc                                                             <  2:312536/107‑1 (MQ=255)                                                                                                                                 ttCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGc                                                             >  1:312536/1‑107 (MQ=255)                                                                                                                                   ccAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCgc                >  2:409856/1‑149 (MQ=255)                                                                                                                                   ccAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGa                <  2:2266588/150‑1 (MQ=255)                                                                                                                                   ccAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGa                >  2:1287084/1‑150 (MQ=255)                                                                                                                                      gCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCt             >  2:1615183/1‑150 (MQ=255)                                                                                                                                               gcCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGCTCTGCCggcggc    >  2:714108/1‑150 (MQ=255)                                                                                                                                                 cACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAg  >  1:2383436/1‑150 (MQ=255)                                                                                                                                                    tGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGc                                                         >  1:1686816/1‑92 (MQ=255)                                                                                                                                                    tGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGc                                                         <  2:1686816/92‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCGCTGTCGACCAGCGCGTGGACCAGTGGCTGGTCGCCTTGCACGCGCGCCAGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCGGGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAG  >  NC_002947/3489794‑3490087

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