Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	762,983	T→A	7.2%	I10N (ATT→AAT)	PP_0655 →	fimbrial‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	762,983	0	T	A	7.2%	181.1 / 2.7	64	I10N (ATT→AAT)	PP_0655	fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base A (3/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AGGCGATGTTCGACCACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGTCAATCAAAG  >  NC_002947/762837‑763130                                                                                                                                                   |                                                                                                                                                    aGGCGATGTTCGACCACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTg                                                                                                                                                   <  7:295251/149‑1 (MQ=255) aGGCGATGTTCGACCACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTg                                                                                                                                                   <  2:246358/149‑1 (MQ=255)              ccACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTg                                                                                                                                                   >  8:295299/1‑136 (MQ=255)              ccACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTg                                                                                                                                                   <  7:295299/136‑1 (MQ=255)              ccACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGaccca                                                                                                                                      >  2:337231/1‑149 (MQ=255)                aCTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGacccacc                                                                                                                                    >  3:144088/1‑149 (MQ=255)                 ctctGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGacccaccc                                                                                                                                   >  3:364366/1‑149 (MQ=255)                     gCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGacc                                                                                                                                        >  1:310613/1‑140 (MQ=255)                     gCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGacc                                                                                                                                        <  2:310613/140‑1 (MQ=255)                     gCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTAc                                                                                                                               >  7:109822/1‑149 (MQ=255)                                 gcgcGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGt                                                                                                                    <  3:295437/148‑1 (MQ=255)                                    cgGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAGGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGc                                                                                                                <  4:23426/149‑1 (MQ=255)                                    cgGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGc                                                                                                                <  4:144088/149‑1 (MQ=255)                                      ggATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGTAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCAt                                                                                                              <  5:237765/149‑1 (MQ=255)                                       gATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATg                                                                                                             <  8:109822/149‑1 (MQ=255)                                       gATATTGCCGAGATGAAGGCGTTGTAGGGCGTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATg                                                                                                             <  4:362608/149‑1 (MQ=255)                                            tGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCttt                                                                                                        >  7:30634/1‑149 (MQ=255)                                              ccGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTc                                                                                                       >  3:73962/1‑148 (MQ=255)                                                     gAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCaa                                                                                                                            >  2:195339/1‑120 (MQ=255)                                                     gAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCaa                                                                                                                            <  1:195339/120‑1 (MQ=255)                                                        ggCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGCCCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCAt                                                                                            >  1:296916/1‑149 (MQ=255)                                                           gTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATcc                                                                                          >  4:173766/1‑148 (MQ=255)                                                           gTTGTAGGGCTTGAAGCGTCATGCATAGTGTGTGGGGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTATGCATGGCttt                                                                                                        <  3:264557/134‑1 (MQ=255)                                                                aGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACACACACTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATc                                                                                     >  4:319290/1‑148 (MQ=255)                                                                    cTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTGTGCCTACGTACCGGTCACACCCGTTCAATGGGGTTAAGACCCACCCTTACCAAGGACAGGTAGGCATTGCTTTAGACCGTTTCATCCAAATagacg                                                                                 >  1:276848/1‑148 (MQ=255)                                                                     ttGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGCCGCTTTCATACAAATcgacga                                                                                >  5:318390/1‑148 (MQ=255)                                                                              cATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTTCCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTgg                                                                       >  5:225920/1‑148 (MQ=255)                                                                              cATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTTCCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTgg                                                                       <  6:225920/148‑1 (MQ=255)                                                                               aTGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGc                                                                      >  6:314035/1‑148 (MQ=255)                                                                                       ggTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTg                                                             <  4:73962/149‑1 (MQ=255)                                                                                             gcggcTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCAt                                                                                                                                                     >  3:219022/1‑55 (MQ=255)                                                                                             gcggcTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCAt                                                                                                                                                     <  4:219022/55‑1 (MQ=255)                                                                                             gcggcTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCCGGCGAAGCGCTGGATGAg                                                       >  4:118796/1‑149 (MQ=255)                                                                                              cggcTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGc                                                      <  4:364366/149‑1 (MQ=255)                                                                                                         gTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCAATGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAAc                                           >  1:168317/1‑149 (MQ=255)                                                                                                         gTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCAATGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAAc                                           >  1:162848/1‑149 (MQ=255)                                                                                                                gAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCaaa                                                                                       >  6:260082/1‑98 (MQ=255)                                                                                                                gAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCaaa                                                                                       <  5:260082/98‑1 (MQ=255)                                                                                                                gAACACTTTGCCTACGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCg                                    <  3:319290/149‑1 (MQ=255)                                                                                                                  acacTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCTACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCgag                                  <  5:81108/149‑1 (MQ=255)                                                                                                                    acTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACTACATCGCTGGCGAATCGCTGGATGAGCAGTACACCAACTGGATCgagct                                >  2:82460/1‑147 (MQ=255)                                                                                                                      tttGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCCAGATc                               >  3:188165/1‑148 (MQ=255)                                                                                                                         gCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAAc                            <  5:314035/148‑1 (MQ=255)                                                                                                                         gCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACg                           >  8:337161/1‑149 (MQ=255)                                                                                                                          ccTCCGTACCTGTCACACCCGTTCAATGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACgg                          <  2:168317/149‑1 (MQ=255)                                                                                                                          ccTCCGTACCTGTCACACCCGTTCAATGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACgg                          <  2:162848/149‑1 (MQ=255)                                                                                                                           cTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGACCAGTACACCAACTGGATCGAGATCAACgg                          >  5:149421/1‑148 (MQ=255)                                                                                                                             ccGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACCCCAACTGGATCGAGATCAACGGCTa                       >  7:351296/1‑149 (MQ=255)                                                                                                                             ccGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCt                        >  7:315568/1‑148 (MQ=255)                                                                                                                              cGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCt                        <  7:337161/147‑1 (MQ=255)                                                                                                                                    ccccaccCCCGTGCATTGGGGTTCAGACCCACCCTTACCCACGCCAGGTACGCATCGCATTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCAAGATCAACGGCTACAACTTc                <  4:169034/146‑1 (MQ=255)                                                                                                                                      tCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCg               <  8:30634/148‑1 (MQ=255)                                                                                                                                       cacaCCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGgt             >  7:307697/1‑149 (MQ=255)                                                                                                                                        acacCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGgt             <  1:337231/148‑1 (MQ=255)                                                                                                                                        acacCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGgt             <  2:296916/148‑1 (MQ=255)                                                                                                                                        acacCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCCAGATCAACGGCTACAACTTCGgt             <  4:188165/148‑1 (MQ=255)                                                                                                                                          acCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGgtgt           >  5:118063/1‑148 (MQ=255)                                                                                                                                            ccGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGtcaa        >  1:353014/1‑149 (MQ=255)                                                                                                                                              gTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAAc                            <  7:334750/127‑1 (MQ=255)                                                                                                                                              gTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAAc                            >  8:334750/1‑127 (MQ=255)                                                                                                                                              gTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCg               <  2:563/140‑1 (MQ=255)                                                                                                                                              gTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCg               >  1:563/1‑140 (MQ=255)                                                                                                                                               ttCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGtcaatca     <  2:353014/149‑1 (MQ=255)                                                                                                                                               ttCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCCGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGtcaatca     <  3:118796/149‑1 (MQ=255)                                                                                                                                                tCATTGGGGTTCTGACCGACCCTTACCAAGTTGATGTAGGGTTGGCTTTAGACGCTTTCATCCAATTCTTCGAGTTCGCTGTCGAAGCGCTGGATGTGCAGTACTCGAACTGTTTCGATATCAACGGCTACAACTTCGGTGtcaatc      <  6:149421/147‑1 (MQ=255)                                                                                                                                                  aTTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGTCAATCaaa   <  8:315568/148‑1 (MQ=255)                                                                                                                                                   ttGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATACAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGTCAATCAAAg  >  1:85583/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGGCGATGTTCGACCACTCTGCAGACCAGGTGGCGCGGGATATTGCCGAGATGAAGGCGTTGTAGGGCTTGAAGCGTCATGCATCGGGTGTGGCGGCTCGGGCAGTGGGGGGAACACTTTGCCTCCGTACCTGTCACACCCGTTCATTGGGGTTCAGACCCACCCTTACCAAGGAGAGGTAGGCATGGCTTTCGACGCTTTCATCCAAATCGACGACATCGCTGGCGAAGCGCTGGATGAGCAGTACACCAACTGGATCGAGATCAACGGCTACAACTTCGGTGTCAATCAAAG  >  NC_002947/762837‑763130

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