Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,388,832	T→C	8.4%	I637V (ATC→GTC)	PP_3865 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,388,832	0	T	C	8.4%	196.4 / 7.6	71	I637V (ATC→GTC)	PP_3865	tail protein
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (4/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

GGTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTG  >  NC_002947/4388684‑4388979                                                                                                                                                     |                                                                                                                                                     ggTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGAt                                                                                                                                                      >  8:281896/1‑149 (MQ=255)  gTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATg                                                                                                                                                     <  5:53153/149‑1 (MQ=255)     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:162479/149‑1 (MQ=255)                           ccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCa                                                                                                                           <  6:111888/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                                                                                     <  4:15585/149‑1 (MQ=255)                                                                 gccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAg                                                                                     <  2:47914/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)                                                                                                                                                     |                                                                                                                                                     GGTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTG  >  NC_002947/4388684‑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: ‑