Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	715,012:1	+T	13.6%	intergenic (+75/+2)	ispH → / ← PP_0607	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	715,012	1	.	T	13.6%	216.9 / 34.0	66	intergenic (+75/+2)	ispH/PP_0607	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base . (29/28);  new base T (6/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGG  >  NC_002947/714871‑715156                                                                                                                                               |                                                                                                                                                 ggACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTagcag                                                                                                                                          >  1:316507/1‑150 (MQ=255)           gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGc                                                                                                                                >  1:68840/1‑150 (MQ=255)                      cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                     >  2:1199295/1‑150 (MQ=255)                       aCCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCg                                                                                                                    >  2:1298094/1‑150 (MQ=255)                         cTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTa                                                                                                                  >  1:842018/1‑150 (MQ=255)                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGGCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                 <  1:1670799/150‑1 (MQ=255)                           tttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAg                                                                                                                <  2:68840/150‑1 (MQ=255)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCTACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                       <  1:301282/150‑1 (MQ=255)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                              <  2:951214/127‑3 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                              <  2:951196/127‑3 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                              >  1:951196/1‑125 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcag                                                                                                                              >  1:951214/1‑125 (MQ=39)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTa                                                                                                                                              >  1:1635983/1‑111 (MQ=255)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTa                                                                                                                                              <  2:1635983/111‑1 (MQ=255)                                    ccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                       <  2:471918/150‑1 (MQ=255)                                     cAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                                                                      >  2:1031303/1‑150 (MQ=255)                                          agCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGcc                                                                                                 <  1:1298094/150‑1 (MQ=255)                                                 ggTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGAt                                                                                          <  2:316507/150‑1 (MQ=255)                                                     cgctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCa                                                                                      <  2:842018/150‑1 (MQ=255)                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAg                                                                                                                            >  2:1298189/1‑111 (MQ=255)                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAg                                                                                                                            <  1:1298189/111‑1 (MQ=255)                                                              aTCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCTTTGCGCGAAACAAGACCGCTCCTGCAATGACAGC‑GCTAGCAGCGGTCATTGCTCATGTCCTCCGTAAGCAACCTGACTCTCCCGCTTGATGCCAGCACCACCt                                                                             <  2:220822/150‑1 (MQ=255)                                                              aTCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTc                                                                                                                              <  2:647055/101‑1 (MQ=255)                                                              aTCTGACGGGCAGTACCGCCGCCAATGCAGGGGATGCTGTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTc                                                                                                                              >  1:647055/1‑101 (MQ=255)                                                               tCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTg                                                                            >  2:1516790/1‑150 (MQ=255)                                                                     gggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                      >  1:1084229/1‑150 (MQ=255)                                                                      ggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                     <  1:488900/150‑1 (MQ=255)                                                                         aGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTg                                                                  >  1:705840/1‑150 (MQ=255)                                                                         aGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTg                                                                  <  1:985926/150‑1 (MQ=255)                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                 >  2:978382/1‑150 (MQ=255)                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                 >  1:240938/1‑150 (MQ=255)                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATACCAGCACCACCTGATATCGGCTGc                                                                 >  1:1734360/1‑150 (MQ=255)                                                                           tACCGCCGCAATTGAAGGGGCTGCTTGGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCg                                                                >  2:1392161/1‑150 (MQ=255)                                                                            cccgacgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCTACACAAGACCGCTCCTGCAATGACCGC‑GCTAGGAGCAGTCATTTCTCAGGTCCTCTGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCGCCTGATATCGGCTGCGa                                                               <  2:1740389/149‑1 (MQ=255)                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCg                                                                                                >  2:1100887/1‑116 (MQ=255)                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCg                                                                                                <  1:1100887/116‑1 (MQ=255)                                                                                 cgccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagcagc                                                          <  2:625355/150‑1 (MQ=255)                                                                                      aTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGgcgc                                                     <  2:705840/150‑1 (MQ=255)                                                                                        gCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGt                                                   <  1:1516790/150‑1 (MQ=255)                                                                                         cAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTCCACCGTTCCTCCGTAAGCGCCCGGAGTCGGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGtt                                                  <  1:236436/150‑1 (MQ=255)                                                                                          aGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                 >  1:1434140/1‑134 (MQ=255)                                                                                          aGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                 <  2:1434140/134‑1 (MQ=255)                                                                                              gctgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCGAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGcac                                             >  2:1726961/1‑150 (MQ=255)                                                                                                tgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTGAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATGGGCTGCGAGCAGCGGCGCGTTAGCACAt                                           >  1:1788923/1‑150 (MQ=255)                                                                                                tgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACAt                                           >  1:1402980/1‑150 (MQ=255)                                                                                                     tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCa                                      <  1:1199295/150‑1 (MQ=255)                                                                                                           ccTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTg                                >  2:646244/1‑150 (MQ=255)                                                                                                              tGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTgcgc                             <  2:1402980/150‑1 (MQ=255)                                                                                                                      acaAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAg                     >  2:617095/1‑150 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCGGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                             >  1:1851493/1‑73 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCtt                                                                                             <  2:1851493/73‑1 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGtt                <  1:1536393/150‑1 (MQ=255)                                                                                                                           aCCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGtt                >  1:189375/1‑150 (MQ=255)                                                                                                                               cTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc            >  1:376821/1‑150 (MQ=255)                                                                                                                                tCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCg           <  1:1031303/150‑1 (MQ=255)                                                                                                                                 tgagCAATGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGa          >  1:1565167/4‑150 (MQ=255)                                                                                                                                   tGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGa                                                               >  2:1897516/1‑95 (MQ=255)                                                                                                                                   tGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGa                                                               <  1:1897516/95‑1 (MQ=255)                                                                                                                                   tGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACg        >  1:498037/1‑150 (MQ=255)                                                                                                                                     cAATGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACGACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTc      >  2:771809/1‑150 (MQ=255)                                                                                                                                      aaTGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                      <  1:1448624/85‑1 (MQ=38)                                                                                                                                      aaTGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                      >  2:1448624/1‑85 (MQ=38)                                                                                                                                      aaTGACCGC‑GCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                     >  1:1501705/1‑38 (MQ=255)                                                                                                                                      aaTGACCGC‑GCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                     <  2:1501705/38‑1 (MQ=255)                                                                                                                                       aTGACCGCTGCTAGC‑GCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACTTTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCAt    >  2:118813/1‑150 (MQ=255)                                                                                                                                         gACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCAAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATgg  <  1:646244/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 GGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGG  >  NC_002947/714871‑715156

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