Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,388,834	G→A	8.5%	T636I (ACC→ATC) ‡	PP_3865 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,388,834	0	G	A	8.5%	197.0 / 6.4	70	T636I (ACC→ATC) ‡	PP_3865	tail protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (4/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

TTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTG  >  NC_002947/4388688‑4388979                                                                                                                                                   |                                                                                                                                                   ttCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGt                                                                                                                                                   <  6:161388/148‑1 (MQ=255)    gTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACg                                                                                                                                              >  1:47914/1‑149 (MQ=255)    gTGCAGCTCGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGACTACCATGGTTCGCCCCCCGACATTCGGCGTCATGGAC‑ACg                                                                                                                                              >  8:224950/1‑149 (MQ=255)         gCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGACCATC‑ACGCCgaa                                                                                                                                         <  6:375171/149‑3 (MQ=255)              gTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGAc                                                                                                                                    <  2:182298/149‑1 (MQ=255)                 cgcCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTg                                                                                                                                  <  5:12265/148‑1 (MQ=255)                  gccggccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGACTCGCGTCGCAAATCCCGCCATTCACTCATAGGTGCCTCCCATCGTTCTCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTgct                                                                                                                                >  8:10877/1‑147 (MQ=255)                   ccggccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGtt                                                                                                                               <  1:271088/149‑1 (MQ=255)                       ccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCa                                                                                                                           <  6:111888/149‑1 (MQ=255)                       ccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCa                                                                                                                           <  6:162479/149‑1 (MQ=255)                        cgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAc                                                                                                                          >  7:125434/1‑149 (MQ=255)                         gACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAcc                                                                                                                         >  8:267125/1‑149 (MQ=255)                           ccGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGa                                                                                                                       >  4:199250/1‑149 (MQ=255)                            cGGCACGCTGTCGCTGTCAACCCAGAACGGGACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGAc                                                                                                                      >  1:35441/1‑149 (MQ=255)                            cGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGAc                                                                                                                      >  2:229868/1‑149 (MQ=255)                            cGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGAc                                                                                                                      <  3:294879/149‑1 (MQ=255)                                aCGCTGTCGCTGTCAACCCAGAACGGCACGGCGCGAAACGACGCCTGGCGTCGCAAATCACGCCATTCACTGATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑TCGCGGGGCTGACTGGTTTTCACCGAGTCGa                                                                                                                  <  7:224950/149‑1 (MQ=255)                                   ctgtcgctgtcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGacac                                                                                                               <  7:281896/149‑1 (MQ=255)                                      tcgctgtcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACg                                                                                                              >  1:150436/1‑147 (MQ=255)                                        gctgtcAACCCAGAACGGAACTCCGCGAAACGACGCCTCGCGTGGCAAATCACGCCATTCACTCTTAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGcag                                                                                                           <  7:10877/148‑1 (MQ=255)                                           gtcAACCAAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGAGCCTCCCATCGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGCCTGGTTTTCCACGACTCGACACGcaggc                                                                                                         >  7:143066/1‑147 (MQ=255)                                            tcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGcaggccag                                                                                                      <  2:19277/149‑1 (MQ=255)                                            tcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGcaggcca                                                                                                       <  8:156536/148‑1 (MQ=255)                                               aCCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAc                                                                                                                          <  6:36614/126‑1 (MQ=255)                                               aCCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAc                                                                                                                          >  5:36614/1‑126 (MQ=255)                                                 ccAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGGCTGGTTTTCGCCGACTCGACCCGCAGGCCCggcggc                                                                                                 >  8:276063/1‑149 (MQ=255)                                                  cAGAGCGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGACCATC‑ACGcc                                                                                                                                            <  3:358686/105‑1 (MQ=255)                                                  cAGAGCGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGACCATC‑ACGcc                                                                                                                                            >  4:358686/1‑105 (MQ=255)                                                   aGAACGGCCCGCCGCCAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGa                                                                                               <  2:324036/149‑1 (MQ=255)                                                       cGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTa                                                                                           >  3:18799/1‑149 (MQ=255)                                                          cACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATc                                                                                         >  8:190140/1‑148 (MQ=255)                                                             gccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAg                                                                                     <  2:47914/149‑1 (MQ=255)                                                             gccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAg                                                                                     <  4:15585/149‑1 (MQ=255)                                                              ccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGCCATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAgg                                                                                    >  5:375065/1‑149 (MQ=255)                                                                    aaCGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCAccc                                                                              >  4:332386/1‑149 (MQ=255)                                                                        acgCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCATCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCAc                                                                          <  1:289685/149‑1 (MQ=255)                                                                          gCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCt                                                                        >  2:386683/1‑149 (MQ=255)                                                                           ccTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTc                                                                       <  6:19602/149‑1 (MQ=255)                                                                                    gCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAg                                                              <  2:150436/149‑1 (MQ=255)                                                                                         tCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgct                                                         >  4:204264/1‑149 (MQ=255)                                                                                           aCGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgctgc                                                       <  7:190140/149‑1 (MQ=255)                                                                                             gCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCTTCATTCGGCGTGATGGTC‑AGGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAATGCCAGCCGCACCTCCCCTTCCATCTGCTGCtt                                                     <  2:14330/149‑1 (MQ=255)                                                                                             gCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTCGACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgctgct                                                      >  1:305142/1‑149 (MQ=255)                                                                                             gCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgctgcat                                                     >  4:254382/1‑147 (MQ=255)                                                                                               cATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGc                                                   >  5:21241/1‑149 (MQ=255)                                                                                                aTTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGcc                                                  >  4:367685/1‑149 (MQ=255)                                                                                                aTTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGc                                                   >  7:287988/1‑148 (MQ=255)                                                                                                 ttCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGcca                                                 <  3:204264/149‑1 (MQ=255)                                                                                                    aCTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGccacc                                               >  6:325617/1‑148 (MQ=255)                                                                                                        aTAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCCCCCACCt                                          >  4:96950/1‑149 (MQ=255)                                                                                                        aTAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCt                                          <  8:287988/149‑1 (MQ=255)                                                                                                        aTAGGTGCCTCCCATGGTTAGACGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTTCTTGccacccac                                            <  6:159498/147‑1 (MQ=255)                                                                                                           ggTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTgc                                        <  4:18799/148‑1 (MQ=255)                                                                                                               ccTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCCCCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGa                                   >  2:89819/1‑149 (MQ=255)                                                                                                                   ccATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGTCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGAtt                                 <  7:267125/147‑1 (MQ=255)                                                                                                                   ccATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGAtt                                 <  3:332386/147‑1 (MQ=255)                                                                                                                   ccATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGa                               >  8:95579/1‑149 (MQ=255)                                                                                                                   ccATGGTTCGCCGCCC‑ACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGAc                              >  5:176904/1‑149 (MQ=255)                                                                                                                    cATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGAc                              >  4:49666/1‑149 (MQ=255)                                                                                                                        gTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAg                          <  1:229868/149‑1 (MQ=255)                                                                                                                         ttCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAg                          >  8:137667/1‑148 (MQ=255)                                                                                                                           cgccgcCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTa                        >  6:250311/1‑148 (MQ=255)                                                                                                                            gccgccCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTaaa                      <  5:325617/149‑1 (MQ=255)                                                                                                                                     agcccGGCGTGACCATC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGc             >  2:131461/5‑149 (MQ=255)                                                                                                                                      gcccGGCGTGACCATC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCt            >  5:78849/4‑149 (MQ=255)                                                                                                                                      gcccGGCGTGACCATC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCt            >  3:160214/4‑149 (MQ=255)                                                                                                                                      aTTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCt            <  3:301552/149‑1 (MQ=255)                                                                                                                                         cGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCgg         >  6:142497/1‑149 (MQ=255)                                                                                                                                         cGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCg          <  8:125434/148‑1 (MQ=255)                                                                                                                                                aTGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTg  >  6:28400/1‑149 (MQ=255)                                                                                                                                                   gTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTaaa                      >  2:213810/1‑126 (MQ=255)                                                                                                                                                   gTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTaaa                      <  1:213810/126‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTG  >  NC_002947/4388688‑4388979

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