Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,239	Δ1 bp	6.6%	intergenic (‑133/+43)	PP_4606 ← / ← PP_4607	ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,229,239	0	G	.	6.6%	190.0 / 7.8	61	intergenic (‑133/+43)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base . (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGC  >  NC_002947/5229091‑5229385                                                                                                                                                     |                                                                                                                                                     ggggggTGAACTGCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACCCGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                                     <  1:324021/148‑1 (MQ=255)     ggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                                                 >  2:1193592/1‑150 (MQ=255)      gTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTCCCTGGAAACATTTTCGCTGTGGCTGCCTCTTCGCGGGTAAACCCGCCCCTa                                                                                                                                                >  2:919943/1‑150 (MQ=255)            ttCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGAGTTTCCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                         <  1:919943/150‑1 (MQ=255)                 ggACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTa                                                                                                                                    <  2:140487/150‑1 (MQ=255)                  gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAc                                                                                                                                   >  1:32470/1‑150 (MQ=255)                   aCGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAcc                                                                                                                                  >  1:1630520/1‑150 (MQ=255)                      cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGt                                                                                                                               <  1:1193592/150‑1 (MQ=255)                           cTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGTCGGCCTCGTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAaga                                                                                                                          >  2:100273/1‑150 (MQ=255)                           cTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAaga                                                                                                                          >  2:755714/1‑150 (MQ=255)                             gtgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGGTGCCACGTGCGGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTAAAGAgg                                                                                                                        >  1:1604956/1‑150 (MQ=255)                                      tGATTATGGCGAGTCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCagag                                                                                                               <  1:337593/150‑1 (MQ=255)                                       gATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              >  1:942040/1‑150 (MQ=255)                                       gATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              >  1:606355/1‑150 (MQ=255)                                           attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAAc                                                                                                          >  1:301601/1‑150 (MQ=255)                                              gCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTAc                                                                                                       <  2:606355/150‑1 (MQ=255)                                                  atcaatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑aa                                                                                                  >  2:1098073/1‑150 (MQ=255)                                                       tcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCCCCTGTGCCTGGAAATCTCCTTGCTGTTGCTGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTGCCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAgg                                                                                             <  1:851511/150‑1 (MQ=255)                                                       tcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAgg                                                                                             >  1:742461/1‑150 (MQ=255)                                                                 aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcaca                                                                                                                <  2:376109/122‑1 (MQ=255)                                                                 aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcaca                                                                                                                >  1:376109/1‑122 (MQ=255)                                                                   tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCGACGGG‑GGGTTTACCCGTGAAGAGGCTGGAACAGGCAACTACCG‑AACCAGGGTCAAACGcccc                                                                                 <  2:679973/150‑1 (MQ=255)                                                                          acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACgg                                                                          >  2:471484/1‑150 (MQ=255)                                                                              aTGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTg                                                                      >  1:1329863/1‑150 (MQ=255)                                                                               tGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGc                                                                     <  2:52606/150‑1 (MQ=255)                                                                                        cTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaaaa                                                            <  2:1329863/150‑1 (MQ=255)                                                                                         tGTTACCTCACCTGTTCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaa                                                               <  1:164818/146‑1 (MQ=255)                                                                                         tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaa                                                               >  2:164818/1‑146 (MQ=255)                                                                                         tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                                           <  2:742461/150‑1 (MQ=255)                                                                                          gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgg                                                                                                                                             >  2:1215846/1‑69 (MQ=255)                                                                                          gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgg                                                                                                                                             <  1:1215846/69‑1 (MQ=255)                                                                                          gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                          >  1:866479/1‑150 (MQ=255)                                                                                          gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                          <  1:755714/150‑1 (MQ=255)                                                                                           ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAccc                                                         >  1:709329/1‑150 (MQ=255)                                                                                              ccTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                                 <  2:52277/139‑1 (MQ=255)                                                                                              ccTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                                 >  1:52277/1‑139 (MQ=255)                                                                                                tCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCa                                                                                               <  2:601941/107‑1 (MQ=255)                                                                                                tCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCa                                                                                               >  1:601941/1‑107 (MQ=255)                                                                                                  aCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                                  >  1:8883/1‑150 (MQ=255)                                                                                                     tgtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACg                                                                           >  2:944969/1‑122 (MQ=255)                                                                                                     tgtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACg                                                                           <  1:944969/122‑1 (MQ=255)                                                                                                          cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                          <  2:709329/150‑1 (MQ=255)                                                                                                           tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAACAGGCGGGCCCAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                                         >  1:667638/1‑150 (MQ=255)                                                                                                                     tttGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcgc                               <  2:301601/150‑1 (MQ=255)                                                                                                                            tgGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCagcagc                        >  1:781054/1‑150 (MQ=255)                                                                                                                               cggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGc                     <  2:866479/150‑1 (MQ=255)                                                                                                                                 gcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTc                   <  2:942040/150‑1 (MQ=255)                                                                                                                                 gcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTc                   >  1:132428/1‑150 (MQ=255)                                                                                                                                  ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  <  2:1417663/150‑1 (MQ=255)                                                                                                                                  ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  >  2:687779/1‑150 (MQ=255)                                                                                                                                  ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  <  2:32470/150‑1 (MQ=255)                                                                                                                                  ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  <  2:1630520/150‑1 (MQ=255)                                                                                                                                   ctctTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCgg                 >  2:294708/1‑150 (MQ=255)                                                                                                                                   ctctTCACGGGTAAACCC‑CCCGTA‑GGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                                  >  2:1069864/1‑115 (MQ=38)                                                                                                                                   ctctTCACGGGTAAACCC‑CCCGTA‑GGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                                  <  1:1069864/115‑1 (MQ=38)                                                                                                                                    tctTCGCGAGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                >  1:711328/1‑150 (MQ=255)                                                                                                                                      ttCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCa              >  2:312504/1‑150 (MQ=255)                                                                                                                                             gTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGaa       >  1:1455062/1‑150 (MQ=38)                                                                                                                                              tAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAg      >  1:643807/1‑150 (MQ=255)                                                                                                                                                  cccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc  <  2:8883/150‑1 (MQ=255)                                                                                                                                                  cccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAAC‑ACCGCCTCCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc  <  1:100273/150‑1 (MQ=255)                                                                                                                                                     gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGc                                                                                    <  2:562563/65‑1 (MQ=255)                                                                                                                                                     gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGc                                                                                    >  1:562563/1‑65 (MQ=255)                                                                                                                                                     gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                >  1:740455/1‑133 (MQ=255)                                                                                                                                                     gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                <  2:740455/133‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGC  >  NC_002947/5229091‑5229385

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