Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,355,475	Δ1 bp	10.8%	intergenic (‑93/‑1)	nrdA ← / → PP_1180	ribonucleoside‑diphosphate reductase subunit alpha/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,355,475	0	A	.	10.8%	185.9 / 18.5	65	intergenic (‑93/‑1)	nrdA/PP_1180	ribonucleoside‑diphosphate reductase subunit alpha/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base . (3/4);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

TCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAA  >  NC_002947/1355327‑1355622                                                                                                                                                       |                                                                                                                                                     tCGGGGGCCTGCGGCACCTGCGCGTGCGTGTTGTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCaa                                                                                                                                                     <  2:330841/150‑1 (MQ=255)    gcggcCTGCGTCACCTGGGCCTGCGCGTACTCGCCCGTCGTTTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGCatg                                                                                                                                                    <  2:1182648/150‑3 (MQ=255)        ccTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgca‑tgc                                                                                                                                             <  1:81774/150‑1 (MQ=255)        ccTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTGTGATATTGTCCCGCCTCGATATGGTCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGCTACGCAatgca‑tgc                                                                                                                                             >  1:461808/1‑150 (MQ=255)         cTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑Tgcg                                                                                                                                            >  1:874981/1‑150 (MQ=255)          tGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑Tgcgc                                                                                                                                           <  2:294313/150‑1 (MQ=255)           gCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgca                                                                                                                                                 <  2:1270721/144‑1 (MQ=255)           gCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgca                                                                                                                                                 >  1:1270721/1‑144 (MQ=255)            cGGCACTTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGC‑ATGCATTGCGCAt                                                                                                                                         <  2:476110/150‑1 (MQ=255)            cGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCAt                                                                                                                                         <  2:1246146/150‑1 (MQ=255)                aCCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGC‑ATGCATTGCGCATcc                                                                                                                                       >  2:52619/1‑150 (MQ=255)                      gCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGC‑ATGCATTGCGCATccccgcga                                                                                                                                 <  1:52619/150‑7 (MQ=38)                        cTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCgg                                                                                                                             >  1:849285/1‑150 (MQ=255)                         tGCGGGTTCTCGCGCGGCGTGTCGGATTGCATGGCAGGGAGTCTCCATAGTCTCTGTATTGCTCAGGTCCCTTATGGGCGCGCACCGTTCCGTCCACTTGCTCGATAACCTCGCCG‑‑GCATGCGCCCGGCA‑TGCGCACTAGCCCGCGCggg                                                                                                                            >  2:410280/1‑150 (MQ=255)                                  tcGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTAGCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGAc                                                                                                                   >  1:1157957/1‑150 (MQ=255)                                       aGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCgg                                                                                                              <  1:82232/150‑1 (MQ=255)                                         tagtgtCGGGTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCTCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                            <  1:410280/148‑1 (MQ=255)                                         tCGTGTCTGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                            <  1:81631/150‑1 (MQ=255)                                           gtgtCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACa                                                                                                          >  2:16580/1‑150 (MQ=255)                                             gtCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                        <  2:461808/150‑1 (MQ=255)                                                    tGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGa                                                                                                 <  1:988870/150‑1 (MQ=255)                                                      cATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTa                                                                                                                        <  1:504882/125‑1 (MQ=255)                                                      cATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTa                                                                                                                        >  2:504882/1‑125 (MQ=255)                                                              ggtgTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                               <  1:741968/126‑1 (MQ=255)                                                              ggtgTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                               >  2:741968/1‑126 (MQ=255)                                                               gtgtCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCaa                                                                                      >  2:217268/1‑150 (MQ=255)                                                                  tctcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GCATGCGCCATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATgg                                                                                   >  2:185451/1‑150 (MQ=255)                                                                    tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                 >  2:1365077/1‑150 (MQ=255)                                                                    tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCACCTGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                 >  1:86790/1‑150 (MQ=255)                                                                      cacaGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAACGGCTTg                                                                               >  2:1363205/1‑150 (MQ=255)                                                                         aGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCCCTCGCGCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCg                                                                            <  1:797629/150‑1 (MQ=255)                                                                          gTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGc                                                                           >  2:16160/1‑150 (MQ=255)                                                                             tCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCTCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGCAATAAGGGGCAATGGCTTGCCGCTcc                                                                        <  1:403141/150‑1 (MQ=255)                                                                                     gTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGa                                                                >  2:1320271/1‑150 (MQ=255)                                                                                      ttCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                            <  2:1050659/105‑1 (MQ=255)                                                                                      ttCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                            >  1:1050659/1‑105 (MQ=255)                                                                                            cgcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCaaa                                                         <  1:127046/150‑1 (MQ=255)                                                                                             gcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAg                                                        <  1:281063/150‑1 (MQ=255)                                                                                               cTTATGGGCGCCCACCGTTCCGTCCATTTGCTAGATAACCTCGCGG‑‑GGATGCGGAATGCA‑TGCGCATCCGCCCGAGCGGGCGTACAGACTTCGGGACAGACTAAGGAATAAGGGGCAATGGATTTCCGCTCCCTGGCAGAAGTTAAAGGt                                                      >  1:930977/1‑150 (MQ=255)                                                                                                ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                            <  1:82472/95‑1 (MQ=255)                                                                                                ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                            >  2:82472/1‑95 (MQ=255)                                                                                                ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAg                                                              >  2:282721/1‑141 (MQ=255)                                                                                                ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAg                                                              <  1:282721/141‑1 (MQ=255)                                                                                                ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGtt                                                     >  1:1219823/1‑150 (MQ=255)                                                                                                ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGtt                                                     >  1:729989/1‑150 (MQ=255)                                                                                                    gggCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTgg                                                 >  2:1029466/1‑150 (MQ=255)                                                                                                       cgcCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGcc                                              <  1:1365077/150‑1 (MQ=255)                                                                                                                ttCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATaa                                                                                             >  2:562966/1‑94 (MQ=255)                                                                                                                ttCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATaa                                                                                             <  1:562966/94‑1 (MQ=255)                                                                                                                    gTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACt                                 >  1:1196836/1‑150 (MQ=255)                                                                                                                        aCTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGCCAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTg                             >  1:348717/1‑150 (MQ=255)                                                                                                                         cTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGc                            <  2:874981/150‑1 (MQ=255)                                                                                                                          ctGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGCCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCt                           >  2:974591/2‑150 (MQ=255)                                                                                                                                 aTAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGt                    <  1:217268/150‑1 (MQ=255)                                                                                                                                      cTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAg               <  2:1219823/150‑1 (MQ=255)                                                                                                                                          cgggcggATGCGC‑ATGCATTGCGCATC‑‑CCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAAGAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGTCAggcct           >  2:1024481/6‑150 (MQ=255)                                                                                                                                           gggcggATGCGC‑ATGCATTGCGCATC‑‑CCCGCGAGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctg          <  2:1343128/146‑1 (MQ=255)                                                                                                                                           gg‑‑ggATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctg          <  1:1320271/150‑1 (MQ=255)                                                                                                                                               tgATGCGC‑ATGCATTGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTg                                                  >  2:1260496/2‑108 (MQ=255)                                                                                                                                               tgATGCGC‑ATGCATTGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTg                                                  <  1:1260496/107‑1 (MQ=255)                                                                                                                                               ggATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAACGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggc        <  1:1363205/150‑1 (MQ=255)                                                                                                                                                gATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGAGTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggcc       >  1:259555/1‑150 (MQ=255)                                                                                                                                                     gcAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                               <  1:87905/89‑1 (MQ=255)                                                                                                                                                     gcAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                               >  2:87905/1‑89 (MQ=255)                                                                                                                                                     gcAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCt                                                   <  2:266521/101‑1 (MQ=255)                                                                                                                                                     gcAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCt                                                   >  1:266521/1‑101 (MQ=255)                                                                                                                                                     gcAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCaa  >  1:147546/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                     TCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCA‑TGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAA  >  NC_002947/1355327‑1355622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑