Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	495,078	T→A	7.3%	Q319L (CAA→CTA) ‡	PP_0408 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	495,078	0	T	A	7.3%	240.1 / 5.6	82	Q319L (CAA→CTA) ‡	PP_0408	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/40);  new base A (3/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTGG  >  NC_002947/494933‑495225                                                                                                                                                  |                                                                                                                                                    gTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGct                                                                                                                                                  >  2:232203/1‑149 (MQ=255)    cGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCTCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCtt                                                                                                                                               >  5:270733/1‑149 (MQ=255)       gAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGt                                                                                                                                             >  4:34378/1‑148 (MQ=255)        aaGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAAGAGCAATTGTTc                                                                                                                                           >  5:40976/1‑149 (MQ=255)         aGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAAGAGCAATTGTTcc                                                                                                                                          <  6:40976/149‑1 (MQ=255)          gCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTccg                                                                                                                                         >  7:89357/1‑149 (MQ=255)            aaCCGGGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACATGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAAGAGCAATTGTTccgcc                                                                                                                                       <  8:174638/149‑1 (MQ=255)              ccGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCcgccgc                                                                                                                                     >  8:118257/1‑149 (MQ=255)              ccGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAAGAGCAATTGTTccgccgg                                                                                                                                     <  2:87716/149‑2 (MQ=255)                gcgcgTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCtt                                                                                                                                   <  5:116638/149‑1 (MQ=255)                 cgcgTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCttt                                                                                                                                  <  1:259022/149‑1 (MQ=255)                    ggaCCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCtttgtt                                                                                                                               <  6:88635/147‑1 (MQ=255)                     tACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCtttgttt                                                                                                                              <  3:34378/149‑1 (MQ=255)                                tGCACCAGGCATATGCCATTCTCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtggcatgg                                                                                                                    <  4:4473/148‑1 (MQ=255)                                 gCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtggcatggc                                                                                                                   >  2:101177/1‑148 (MQ=255)                                  cACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGcaa                                                                                                                 >  5:66961/1‑149 (MQ=255)                                   aCCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGcaat                                                                                                                >  5:296009/1‑149 (MQ=255)                                      aGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGcaatcaa                                                                                                             <  2:143926/149‑1 (MQ=255)                                         ccTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGa                                                                                                          <  7:118257/149‑1 (MQ=255)                                          cTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGa                                                                                                          >  3:66968/1‑148 (MQ=255)                                             tGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCtttgtttg                                                                                                                             >  7:352984/1‑126 (MQ=255)                                             tGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCtttgtttg                                                                                                                             <  8:352984/126‑1 (MQ=255)                                               ccATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAAtt                                                                                                     <  2:365119/148‑1 (MQ=255)                                                cATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGt                                                                                                   >  7:281635/1‑149 (MQ=255)                                                cATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGt                                                                                                   <  8:235389/149‑1 (MQ=255)                                                cATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTATGTTTGCTGGCAAGGCAATCAATGACAATTGt                                                                                                   >  7:331440/1‑149 (MQ=255)                                                 aTGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTc                                                                                                  <  3:294198/149‑1 (MQ=255)                                                 aTGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTc                                                                                                  <  4:66968/149‑1 (MQ=255)                                                    aTCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtgg                                                                                                                         >  8:128697/1‑123 (MQ=255)                                                    aTCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCtgg                                                                                                                         <  7:128697/123‑1 (MQ=255)                                                    aTCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCgg                                                                                                <  6:6261/148‑1 (MQ=255)                                                         cAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTgg                                                                                          <  6:66961/149‑1 (MQ=255)                                                          aaGCGGCGAGGCCTTTGTCGCCCCGCCACAGGGATCAAGCCTTACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTTCTGTAATGGCAATCAATGACTATTGTCTGTACTg                                                                                           <  1:282443/147‑1 (MQ=255)                                                            gcggcgcggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGTGCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTa                                                                                       >  1:305496/1‑149 (MQ=255)                                                               gcgcggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATcc                                                                                    <  7:333221/149‑1 (MQ=255)                                                               gcgcggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATcc                                                                                    <  2:226676/149‑1 (MQ=255)                                                                   ggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCt                                                                                <  1:232203/149‑1 (MQ=255)                                                                   ggcCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCt                                                                                >  3:372592/1‑149 (MQ=255)                                                                      cTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGcc                                                                             >  7:8846/1‑149 (MQ=255)                                                                      cTTTGACGCCCCGCCACAGGGAGCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGcc                                                                             <  7:327702/149‑1 (MQ=255)                                                                               cccGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTgg                                                                    <  6:340639/149‑1 (MQ=255)                                                                                     acaGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCATTGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCa                                                              >  8:253088/1‑149 (MQ=255)                                                                                     acaGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATc                                                                                     <  2:210708/126‑1 (MQ=255)                                                                                     acaGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATc                                                                                     <  2:217005/126‑1 (MQ=255)                                                                                     acaGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATc                                                                                     >  1:217005/1‑126 (MQ=255)                                                                                     acaGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATc                                                                                     >  1:210708/1‑126 (MQ=255)                                                                                           aTCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGtt                                                        >  5:336070/1‑149 (MQ=255)                                                                                              aaGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCt                                                     <  5:210630/149‑1 (MQ=255)                                                                                              aaGACTGACTAAGCCATGAGGGATCATCCACCGCCGCACCCCGCGGAACAATTGCTCGTCTGCCGCCGCTTTGTTTGCTGGCATGGCAATCAAGCACAATTGTACGTACTGGCTAACATTATGCCGATTGCTGGCCTGCAGCCGtagat                                                     >  4:23677/1‑145 (MQ=255)                                                                                               aGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTg                                                    >  6:316612/1‑149 (MQ=255)                                                                                                gCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTgg                                                   >  1:130135/1‑149 (MQ=255)                                                                                                gCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTgg                                                   >  2:151302/1‑149 (MQ=255)                                                                                                gCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTgg                                                   >  4:139752/1‑149 (MQ=255)                                                                                                   tGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGCAACAAATCCTCTTGTTCCTCCGCTTTGTTTGCTGGCATGGCAATCAATGACAAATGTCGGTACTGGATATCCTTCTGACGTTTGCTGGCCTTCAGCCGATGCTGGGc                                                 >  2:174706/1‑148 (MQ=255)                                                                                                     aCTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGc                                              <  2:229910/149‑1 (MQ=255)                                                                                                      cTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCg                                             >  4:32719/1‑149 (MQ=255)                                                                                                       tCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTTCTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGATGCTGGGCTGCGt                                            >  8:261984/1‑149 (MQ=255)                                                                                                             cTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGt                                      <  4:298865/149‑1 (MQ=255)                                                                                                                aGGGCTCATCCACCCCCGCACCCGGCGGCACAATTGATCTTGTTCAGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTATGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGGGGTCAg                                   >  4:97657/1‑149 (MQ=255)                                                                                                                    atCATCCACCCCCGCACCCGGCGGAACAATTGGTCTTGTTCCGCCGCGTTGTTTGCTGGCATGGCAATCATTGACAATTGTCGGTACTGGCTATCCTTCTTCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCt                               <  3:23677/148‑1 (MQ=255)                                                                                                                       aTCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGc                             <  8:281635/148‑1 (MQ=255)                                                                                                                           aCCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAgcgc                                <  2:53440/141‑1 (MQ=255)                                                                                                                           aCCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCtagatcgg                        >  1:53440/1‑142 (MQ=255)                                                                                                                             ccccGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGc                      >  8:359427/1‑149 (MQ=255)                                                                                                                              cccGCACCCGGGGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCg                     <  6:270733/149‑1 (MQ=255)                                                                                                                              cccGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCg                     <  7:359427/149‑1 (MQ=255)                                                                                                                               ccGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCg                     >  5:2704/1‑148 (MQ=255)                                                                                                                                  cACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtag                  >  5:137941/1‑148 (MQ=255)                                                                                                                                   aCCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtagta                <  2:130135/149‑1 (MQ=255)                                                                                                                                    cccGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGGCAGCGCTGCGGGTGGCGtagtag               <  6:178642/149‑1 (MQ=255)                                                                                                                                     ccGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCgtag                  >  3:382099/1‑145 (MQ=255)                                                                                                                                     ccGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATCACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGa              <  8:331440/149‑1 (MQ=255)                                                                                                                                       ggcggAACAAGAGCAATTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGata            >  7:179557/1‑149 (MQ=255)                                                                                                                                       ggcggAACAAGAGCAATTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGata            >  7:175568/1‑149 (MQ=255)                                                                                                                                        gcggAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGatat           <  5:316612/149‑1 (MQ=255)                                                                                                                                           gAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCg                                        <  4:378588/117‑1 (MQ=255)                                                                                                                                           gAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCg                                        >  3:378588/1‑117 (MQ=255)                                                                                                                                              cAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGATATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTTCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCg     >  2:63461/1‑149 (MQ=255)                                                                                                                                               aaTTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCTTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGt    <  1:101177/149‑1 (MQ=255)                                                                                                                                               aaTTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGt    <  6:336070/149‑1 (MQ=255)                                                                                                                                               aaTTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGt    <  1:151302/149‑1 (MQ=255)                                                                                                                                                aTTGCTGTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTg   <  3:139752/149‑1 (MQ=255)                                                                                                                                                 ttGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTgg  <  5:224799/149‑1 (MQ=255)                                                                                                                                                 ttGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTgg  <  1:63461/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTGCGTGAAGCAACCGCGCGTACCCGTGAGTTGCACCAGGCCTATGCCATGATCCGCAAGCGGCGCGGCCTTTGACGCCCCGCCACAGGGATCAAGCCTGACTCAGCCCTGAGGGCTCATCCACCCCCGCACCCGGCGGAACAATTGCTCTTGTTCCGCCGCTTTGTTTGCTGGCATGGCAATCAATGACAATTGTCGGTACTGGCTATCCTTCTGCCGTTTGCTGGCCTGCAGCCGTTGCTGGGCTGCGTTGCGGTCAGCGCTGCGGGTGGCGTAGTAGATATCGGCCGTGG  >  NC_002947/494933‑495225

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