Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,402,311	A→T	5.3%	L19Q (CTG→CAG)	PP_2103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,402,311	0	A	T	5.3%	231.8 / 2.6	75	L19Q (CTG→CAG)	PP_2103	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base T (2/2);  total (36/39)
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.13e-01

GCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGG  >  NC_002947/2402164‑2402450                                                                                                                                                      |                                                                                                                                            gCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGcagcag                                                                                                                                            >  3:27119/1‑149 (MQ=255)    ccGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAgcgc                                                                                                                                         <  3:68032/149‑1 (MQ=255)     cGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGcc                                                                                                                                        <  7:31567/149‑1 (MQ=255)       aTGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                      <  4:138140/149‑1 (MQ=255)           cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                   >  7:140991/1‑148 (MQ=255)             cacaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCat                                                                                                                                <  8:289316/149‑2 (MQ=255)             cacaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAg                                                                                                                                <  4:27119/149‑1 (MQ=255)              acaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatg                                                                                                                               <  2:342914/149‑3 (MQ=255)              acaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGt                                                                                                                               >  3:57338/1‑149 (MQ=255)                 gTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcg                                                                                                                            >  7:256325/1‑149 (MQ=255)                    gcgcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                         >  4:48312/1‑149 (MQ=255)                    gcgcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgc                                                                                                                         <  8:71313/149‑1 (MQ=255)                       cATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCggg                                                                                                                      <  8:205621/149‑1 (MQ=255)                            ttGATACACAGGCCTTCAAGTAACGAGCTCCCAATACCGCGCACGCGGCGCAGCGCTTCACGCAACTCACCGGCCTCATTACGGATTTCATGGCGGAGCAT‑TTCAATTCTC‑AGGCGCAGCAGCGCCTTGCGCAGAACGCTCGGCTTAcg                                                                                                                 >  4:272539/1‑149 (MQ=255)                                 aCCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAcgc                                                                                                                <  5:264741/145‑1 (MQ=255)                                 aCCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAcgc                                                                                                                >  6:264741/1‑145 (MQ=255)                                  cccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGtt                                                                                                           >  1:361841/1‑149 (MQ=255)                                   ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAg                                                                                                                                <  8:221046/127‑1 (MQ=255)                                   ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAg                                                                                                                                >  7:221046/1‑127 (MQ=255)                                    cAGGCCTTCATGTAGCGATCTGCCCATACCGCGCACGCGGCGCATCGGGTCAACCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATCTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGtt                                                                                                                    <  1:276081/139‑1 (MQ=255)                                    cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGtt                                                                                                                    >  2:276081/1‑138 (MQ=255)                                     aGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                        <  8:256325/149‑1 (MQ=255)                                       gCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCa                                                                                                      <  6:343167/149‑1 (MQ=255)                                          ttCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGtt                                                                                                   >  2:171693/1‑149 (MQ=255)                                             aTGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTg                                                                                                >  4:256122/1‑149 (MQ=255)                                              tGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTgg                                                                                               <  4:143163/149‑1 (MQ=255)                                                  gCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAgg                                                                                           >  2:129393/1‑149 (MQ=255)                                                   cGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGc                                                                                          <  2:361841/149‑1 (MQ=255)                                                     aGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGGTTGGCAGTGCTGGCAGGcg                                                                                         >  8:124832/1‑147 (MQ=255)                                                      gCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCa                                                                                       >  3:166755/1‑149 (MQ=255)                                                         gCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATgg                                                                                    >  4:217680/1‑149 (MQ=255)                                                          cccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATg                                                                                     <  4:57338/147‑1 (MQ=255)                                                           ccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGc                                                                                   >  3:203439/1‑148 (MQ=255)                                                                ccGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAgg                                                                             <  1:129393/149‑1 (MQ=255)                                                                ccGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAg                                                                              <  4:380276/148‑1 (MQ=255)                                                                 cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACCGgc                                                                            >  8:8761/1‑149 (MQ=255)                                                                 cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgc                                                                            >  5:146009/1‑149 (MQ=255)                                                                   tgcACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgcgc                                                                          <  6:255406/148‑1 (MQ=255)                                                                     cACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTc                                                                        <  6:146009/149‑1 (MQ=255)                                                                      aCGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTcc                                                                       >  2:13735/1‑149 (MQ=255)                                                                       cgcgGCGCAGCGGTTCAAGCAACTGAACGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTccc                                                                      >  7:138249/1‑149 (MQ=255)                                                                          ggcgCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGt                                                                   >  8:348003/1‑149 (MQ=255)                                                                          ggcgCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGt                                                                   <  3:217680/149‑1 (MQ=255)                                                                           gcgcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGt                                                                   <  4:192459/148‑1 (MQ=255)                                                                           gcgcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGt                                                                   >  3:192459/1‑148 (MQ=255)                                                                              cAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCtt                                                               <  1:28592/149‑1 (MQ=255)                                                                               aGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGAGCAT‑TTCCATCCGC‑AGGCGCAGCAGCCCCTTCCGCAGTTCGCGCGGGTTACGCGTCTTAGGCAGTTCTGGCACGCTCATGGCAACAGACGCTCCCTGTCCTTg                                                              >  7:154934/1‑149 (MQ=255)                                                                                gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTgg                                                             <  4:41940/149‑1 (MQ=255)                                                                                 cGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTgg                                                             >  5:108948/1‑148 (MQ=255)                                                                                 cGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                            <  4:203439/149‑1 (MQ=255)                                                                                         gCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgt                                                    >  8:339282/1‑149 (MQ=255)                                                                                          cAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgtc                                                   <  8:138249/149‑1 (MQ=255)                                                                                                aCCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGGGGGTTACGCGTGTTTGGCAGTTCTTGCAGGCTCATGGCAACAGGCGCTCCCTGTCCtt                                                               <  3:38569/131‑1 (MQ=255)                                                                                                aCCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAgg                                               <  4:166755/147‑1 (MQ=255)                                                                                                aCCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGAGC‑AGGCGCATCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAgg                                               <  1:65528/147‑1 (MQ=255)                                                                                                aCCGGACTCATGACGGATTTCCTGGCGGAGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGGGTTTGGGAGTTCTGGCAGGATCATGGCAACATGCGCTCCCTGTACtt                                                               >  4:38569/1‑131 (MQ=255)                                                                                                 ccGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCGATGCGC‑AGGCGCAGCAGGGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCATGTGCGTCGTGGGTg                                            <  8:140991/149‑1 (MQ=255)                                                                                                 ccGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAgg                                               <  3:56479/146‑1 (MQ=255)                                                                                                 ccGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAgg                                               >  4:56479/1‑146 (MQ=255)                                                                                                     aCTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgcgc                                        <  7:348003/149‑1 (MQ=255)                                                                                                         aTGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                            >  4:130755/1‑125 (MQ=255)                                                                                                         aTGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGc                                                            <  3:130755/125‑1 (MQ=255)                                                                                                          tGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAg                                   >  6:45716/1‑149 (MQ=255)                                                                                                            aCGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAggg                                 <  3:48312/149‑1 (MQ=255)                                                                                                                atgtCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                             <  6:84489/146‑1 (MQ=255)                                                                                                                        gCGGGGCAT‑TTCCATGCTC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtcg                       <  7:363671/147‑1 (MQ=255)                                                                                                                                t‑ttCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccg             >  1:362546/1‑149 (MQ=255)                                                                                                                                       tGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCt       >  5:191845/1‑149 (MQ=255)                                                                                                                                         cgc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt      >  4:33354/1‑148 (MQ=255)                                                                                                                                          gc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt      >  1:113695/1‑147 (MQ=255)                                                                                                                                          gc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt      >  7:168023/1‑147 (MQ=255)                                                                                                                                          gc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt      >  5:237262/1‑147 (MQ=255)                                                                                                                                           c‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg   <  2:215360/149‑1 (MQ=255)                                                                                                                                           c‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg   <  1:171693/149‑1 (MQ=255)                                                                                                                                            aaggcgctgCTGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTGGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg   >  2:311148/2‑149 (MQ=255)                                                                                                                                             aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAgg  <  7:124832/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                            GCACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCAT‑TTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGG  >  NC_002947/2402164‑2402450

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