Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,903	C→A	9.8%	intergenic (‑20/+21)	PP_1879 ← / ← PP_1880	hypothetical protein/outer membrane autotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,104,903	0	C	A	9.8%	182.5 / 9.4	61	intergenic (‑20/+21)	PP_1879/PP_1880	hypothetical protein/outer membrane autotransporter
Reads supporting (aligned to +/- strand):  ref base C (33/22);  new base A (3/3);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCACCGA  >  NC_002947/2104760‑2105048                                                                                                                                                |                                                                                                                                                   cATCGTCCACCACGCTTGGTAGGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCAGCGTCCATCTCCATGTGTCACGCAGGAATCCAGACGGGAAGgcgccgtgg                                                                                                                                               >  6:107948/1‑144 (MQ=255) cATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTga                                                                                                                                               >  8:76133/1‑148 (MQ=255)        caccacGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCCAATGCGCCTTccc                                                                                                                                        <  2:172387/149‑1 (MQ=255)           cacGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCg                                                                                                                                     >  2:178422/1‑149 (MQ=255)            acGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCCAATGCGCCTTCCCttct                                                                                                                                    <  4:215659/149‑4 (MQ=255)                   gTACTGTTGGTGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTGCTGCGGTTAACGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGTGAAGGCGCATTGGCCTTCCCTCGTTTCttt                                                                                                                              <  6:36097/148‑1 (MQ=255)                    tACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTAc                                                                                                                            >  1:336428/1‑149 (MQ=255)                     aCGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTAca                                                                                                                           >  8:142070/1‑148 (MQ=255)                       ggTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAg                                                                                                                         <  1:178422/149‑1 (MQ=255)                         tttggggtACATGAGAGGCCGAGTCCAGTTCAGGGTTTTGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGtt                                                                                                                       <  4:277116/146‑1 (MQ=255)                         ttGGGGGGACATGTGAGGCCGAGTCCAGTTCAGGGGTTGGAGGTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGtt                                                                                                                       <  4:135041/149‑1 (MQ=255)                                    agagagGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGgcgc                                                                                                            >  8:247986/1‑149 (MQ=255)                                           ccGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGc                                                                                                     >  3:177119/1‑148 (MQ=255)                                             gAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCCAATGCGCCTTccc                                                                                                                                        >  1:241087/1‑112 (MQ=255)                                             gAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCCAATGCGCCTTccc                                                                                                                                        <  2:241087/112‑1 (MQ=255)                                                tCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCAATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGc                                                                                               >  8:36333/1‑149 (MQ=255)                                                  cAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGccc                                                                                             >  6:91310/1‑149 (MQ=255)                                                   aGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCCCGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCGGCGTTTCTTTACCACTTGTAACGGGCGCCGA‑CGCTCACGccc                                                                                             >  7:240612/1‑148 (MQ=255)                                                        aGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATgg                                                                                        <  4:22772/148‑1 (MQ=255)                                                        aGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATgg                                                                                        <  7:142070/148‑1 (MQ=255)                                                        aGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATgg                                                                                        <  7:217709/148‑1 (MQ=255)                                                        aGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGc                                                                                       >  2:81527/1‑149 (MQ=255)                                                        aGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGc                                                                                       >  6:221165/1‑149 (MQ=255)                                                        aGGCGTTCGAGATCTGCGGTTACCGGAGCGCTTGTTCTGTGACATCCTCAATCACCATGTGTCACGCATGAATACAGAAGGGAAGGCGCAATCGCCATGCCTCGTTTCTATAACAGAAGTAACGGGCGCCGA‑CGCTCACGCCCCATgca                                                                                       >  4:204305/1‑147 (MQ=255)                                                            gTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGtt                                                                                   >  5:323583/1‑149 (MQ=255)                                                                gAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCaa                                                                                <  1:81527/148‑1 (MQ=255)                                                                 aGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGc                                                                              <  7:36333/149‑1 (MQ=255)                                                                      tGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCAc                                                                                                 <  3:213960/125‑1 (MQ=255)                                                                      tGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCAc                                                                                                 >  4:213960/1‑125 (MQ=255)                                                                              ccGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATgg                                                                  <  5:221165/148‑1 (MQ=255)                                                                               cGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCAt                                                                                          <  8:208838/123‑1 (MQ=255)                                                                               cGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCAt                                                                                          >  7:208838/1‑123 (MQ=255)                                                                               cGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCt                                                                >  6:80734/1‑149 (MQ=255)                                                                                 caGCGCTAGTTCAGTGTCATCGTCCAACTCCATGTGTCACGCATCAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGCGCGCCGA‑CGCACACGCCCCATGGCTGTACAAGCTTGCTGCCATGGCTg                                                               >  4:299281/2‑148 (MQ=255)                                                                                               tgtCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCg                                                >  5:322587/1‑149 (MQ=255)                                                                                                gtCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGACCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGt                                               <  7:76133/149‑1 (MQ=255)                                                                                                   aTCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTggg                                            >  2:272071/1‑149 (MQ=255)                                                                                                   aTCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTgg                                             <  6:323583/148‑1 (MQ=255)                                                                                                      gTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCt                                         <  7:247986/149‑1 (MQ=255)                                                                                                      gTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCt                                         <  5:80734/149‑1 (MQ=255)                                                                                                       tCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTg                                        >  2:334080/1‑149 (MQ=255)                                                                                                       tCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTg                                        >  5:161377/1‑149 (MQ=255)                                                                                                       tCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTg                                        >  2:32284/1‑149 (MQ=255)                                                                                                           tctcCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCAt                                    >  4:184211/1‑149 (MQ=255)                                                                                                               cATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCAttttt                                >  1:298662/1‑149 (MQ=255)                                                                                                                aTGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTc                               <  6:322587/149‑1 (MQ=255)                                                                                                                aTGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTc                               <  1:272071/149‑1 (MQ=255)                                                                                                                aTGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTc                               <  6:161377/149‑1 (MQ=255)                                                                                                                   tgtCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGcc                            >  8:9224/1‑149 (MQ=255)                                                                                                                    gtCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGccc                           >  2:289345/1‑149 (MQ=255)                                                                                                                     tCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCCAAGTCGGCGTGGGCCTGCCATTTTTc                               >  1:153826/1‑144 (MQ=255)                                                                                                                     tCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCCAAGTCGGCGTGGGCCTGCCATTTTTc                               <  2:153826/144‑1 (MQ=255)                                                                                                                         gCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGc                      >  7:90227/1‑149 (MQ=255)                                                                                                                                    gaagggaaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCAccc           >  7:292989/1‑149 (MQ=255)                                                                                                                                      agggaaggCCAATGCGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCt         >  5:65265/1‑149 (MQ=255)                                                                                                                                        ggaaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATc                                                          <  6:298114/98‑1 (MQ=255)                                                                                                                                        ggaaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATc                                                          >  5:298114/1‑98 (MQ=255)                                                                                                                                        ggaaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTg        >  8:269545/1‑148 (MQ=255)                                                                                                                                        ggaaggCCAATGCGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGACCGATCACCCCCCATGGCTGTTAAAGCTGCATCCCATGGCTGAAATCGAATTCGGCTTGGGCCTGCAATTTTACGCCCCAGGCTACCCCCAcccatc        >  1:262626/1‑146 (MQ=255)                                                                                                                                          aaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTgg                                             >  2:9208/1‑109 (MQ=255)                                                                                                                                          aaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTgg                                             <  1:9208/109‑1 (MQ=255)                                                                                                                                          aaggCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCAc     >  8:252286/1‑149 (MQ=255)                                                                                                                                             gcgcATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGc                                       <  4:255606/112‑1 (MQ=255)                                                                                                                                             gcgcATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGc                                       >  3:255606/1‑112 (MQ=255)                                                                                                                                             gcgcATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCACCGa  >  4:329456/1‑149 (MQ=255)                                                                                                                                               gcATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCa                                     >  6:84164/1‑112 (MQ=255)                                                                                                                                               gcATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCa                                     <  5:84164/112‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATGCAGAAGGGAAGGCGCATTGGCCTTCCCTCGTTTCTTTACCAGTTGTAACGGGCGCCGA‑CGCTCACGCCCCATGGCTGTTCAAGCTTGCTGCCATGGCTGTAATCGAAGTCGGCGTGGGCCTGCCATTTTTCGCCCCAGGCGACCGCCACCCCTGCACCGA  >  NC_002947/2104760‑2105048

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